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dvanced  durability  analysis  ^design  tools’Vhave  been  developed  for  metallic  aircraft 
structures.  These  tools  can  be  used  to  evaluate  durability  design  requirements  for  func- 


structures.  These  tools  can  be  used  to  evaluate  durability  design  requirements  tor  tunc- 
tional  impairments  due  to t  excessive  cracking  and  (-27  fuel  leakage/ligament  breakage. 

The  methodology  accounts  for  the  initial  fatigue  quality  variation  of  structural  details, 
the  crack  growth  accumulation  for  a  population  of  structural  details  under  specified  design 
conditions  and  structural  properties.  Step-by-step  procedures  are  provided.  This  volume  is 
limited  to  the  analytical  methods,  technical  aspects,  concepts  and  philosophy  for  the  dur¬ 
ability  analysis  of  metallic  aircraft  structures. 

The  methodology  reflects  a  probabilistic  approach,  a  fracture  mechanics  philosophy  and  both 
deterministic  and  stochastic  crack  growth  methods.  It  can  be  used  to  predict  the  probabil¬ 
ity  of  crack  exceedance  at  any  service  time  and/or  the  cumulative  distribution  of  the  time^ 
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18.  (continued)  probability  of  crack  exceedance,  cumulative  distribution  of  TTCI. 

19.  TconfinuedT^o-crack  initiation  cfTTeif>  at  any  crack  size.  The  methodology  applies  to  the 
small  crack  size  range  associated  with  exp^ssiyie  cracking  (e.g.,  <0.05®)  and  to  large 
through-the-thickness  cracks  (e.g.,  0. 50&-0. 75^  associated  with  fuel  leakage/ligamenc 
breakage . 

No  matter  what  form,  location  or  combination  the  as-manufactured  flaws  may  have  in  fastener 
holes  (e.g.,  scratches,  burrs,  microscopic  imperfections,  etc.)  or  whatever  the  source  of 
fatigue  cracking  may  be,  a  practical  method  of  representing  the  reality  of  the  as- 
manufactured  condition  is  needed  for  durability  analysis.  This  is  taken  care  of  by  the 
equivalent  initial  flaw  concept,  - 

Initial  fatigue  quality  of  a  structural  detail  (e.g.,  fastener  holes,  cutouts,  fillets,  lugs, 
etc.)  is  represented  by  an  equivalent  initial  flaw  size  distribution.  An  equivalent  initial 
flaw  (EIFS)  is  an  artificial  crack  size  which  results  in  an  actual  crack  size  at  an  actual 
point  in  time  when  the  initial  flaw  is  grown  forward.  It  is  determined  by  beck-extrapolating 
fractogr&phic  results.  It  has  the  following  characteristics: 

(1)  An  EIFS  is  an  artificial  crack  assumed  to  represent  the  initial  fatigue  quality  of  a 
structural  detail  in  the  as-manufactured  condition  whatever  the  source  of  fatigue  cracking 
may  be,  (2)  it  has  no  direct  relationship  to  actual  initial  flaws  in  fastener  holes  such  as 
scratches,  burrs,  mlerodafects,  etc.,  and  it  cannot  be  verified  by  NDI,  (3)  it  has  s  universal 
crack  shape  in  which  the  crack  size  is  measured  in  the  direction  of  crack  propagation,  (4) 

EIFSs  are  in  a  fracture  mechanics  format  but  they  are  not  based  on  linear  elastic  fracture 
mechanics  (LEFM)  principles  and  are  not  subject  to  such  laws  or  limitations  such  as  the  "short 
crack  effect,"  (5)  it  depends  on  the  fractographic  data  used,  the  fractographic  crack  size 
range  used  for  the  back-extrapolation  and  the  crack  growth  rate  model  used, (6)  it  must  be 
grown  forvsrd  in  a  manner  consistent  with  the  basis  for  the  EIFS,  (7)  EIFSs  ere  not  unique  - 
a  different  set  is  obtained  for  each  crack  growth  law  used  :or  the  back-extrapolation,  and 
(8)  EIFSs  are  not  necessarily  comparable  with  minimum  crack  sizes  obtained  from  fracto¬ 
graphic  results  from  tear-down  inspections. 

Recommendations  for  durability  analysis  are  as  follows:  (1)  define  the  equivalent  initial 
flaw  size  distribution  (EIFSO)  using  fractographic  data  in  the  small  crack  size  region 
(e.g.,  0.01"-0.05") ,  (2)  use  fractographic  data  pooling  procedure  and  statistical  scaling 
technique  to  estimate  the  EIFSD  parameters  in  a  "global  sense"  for  a  "single  hole  popula¬ 
tion"  basis  (3)  use  the  deterministic  crack  growth  approach  (DCGA)  in  the  small  crack  size 
region  and  (4)  use  the  two-segment  deterministic-stochastic  crack  growth  approach  (DCGA- 
SCGA)  for  applications  in  the  large  crack  size  region  (e.g.,  0.50"-0.75";  the  two-segment 
deterministic  crack  growth  approach  (DCGA-DCGA)  is  also  reasonable  but  it  is  slightly  less 
conservative  than  the  DCGA-SCGA. 

Procedures  have  been  developed  for  defining  initial  fatigue  quality.  These  procedures  could 
be  used  to  standardize  the  way  initial  flaw  sizes  are  determined  from  fractographic  data. 

A  better  understanding  of  initial  flaw  sizas  (i.a.,  what  they  are  and  limitations)  has  been 
developed.  For  consistent  durability  analysis  pradictions,  equivalent  initial  flaws  must 
be  used  in  the  same  context  for  which  they  were  defined.  This  means  that  equivalent  initial 
flaws  must  be  grown  forward  in  the  same  manner  the  EIFSs  were  established  by  beck-extrapolat¬ 
ing  fractographic  results. 
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FOREWORD 


This  rsport  was  prsparsd  by  Gsnsral  Dynasties,  Fort  Worth 
Division,  undar  ths  "Advancsd  Durability  Analysis"  program 
(Air  Fores  Contract  F33615-84-C-3208)  for  ths  Air  Fores 
Wright  Aeronautical  Laboratories  (AFWAL/FIBEC) .  Margery  E. 

Art ley  was  ths  Air  Fores  Project  Engineer;  Dr.  John  W. 

Lincoln  of  ASD/ENFS  and  James  L.  Rudd  of  AFWAL/FIBEC  were 
technical  advisors.  Dr.  S.  D.  Manning  of  ths  Gsnsral 
Dynamics'  Structures  Technology  Staff  was  ths  program  manager 
and  co-principal  investigator  along  with  Dr.  J.  N.  Yang  of 
United  Analysis  Incorporated  (Vienna,  VA). 

This  report  (Vol.  I)  documents  the  analytical  methods 
and  theoretical  aspects  developed  under  the  "Advanced 
Durability  Analysis"  program.  Other  volumes  for  this  program 
are  as  follows: 

o  Volume  II  -  Analytical  Predictions,  Test  Results, 

and  Analytical/Experimental  Correlations 

o  Volume  III  -  Fractographic  Test  Data 

o  volume  IV  -  Executive  Summary 

o  Volume  V  -  Durability  Analysis  Software  User's  Guide 
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SECTION  I 


INTRODUCTION 
1.1  OBJECTIVE  AND  SCOPE 

This  program  had  thraa  objaotivaa:  (1)  racommand  im- 
provamanta  to  the  currant  Air  Forca  durability  dasign  re¬ 
quirements  (MIL-A8866B) [1] ,  (2)  davalop  a  probabiliatic  dur¬ 
ability  analyaia  mathod  capable  of  predicting  tha  durability 
of  advanced  metallic  aircraft  atructuraa  for  such  functional 
impairment  dafinitiona  aa  axcaaaiva  cracking,  fuel  leakage 
and  ligament  breakage,  and  (3)  update  tha  currant  Air  Force 
Durability  Daaign  Handbook  (AFWAL-TR-83-3027) [2 ] .  Fatigue 
cracking  is  the  form  of  degradation  to  be  considered. 

This  report  documents  the  advanced  durability  analysis 
methodology  for  metallic  aircraft  structures  developed  under 
the  Phase  I  effort.  Phase  I  included  three  tasks  as  follows: 
(1)  Task  I  -  refine,  optimize  and  increase  confidence  in  the 
method  used  to  determine  the  equivalent  initial  flaw  size 
distribution  (EIFSD)  for  clearance-fit  fastener  holes,  (2) 
Task  II  -  refine  procedure  and  guidelines  for  predicting 
crack  exceedance  in  the  small  crack  size  region  (e.g., 
<0.10”),  and  (3)  Task  III  -  extend  the  prototype  durability 
analysis  methodology  to  large  crack  sizes  (e.g.,  0.5"  -  0.75" 
through-the-thickness)  associated  with  functional  impairment 
due  to  fuel  leakage  and  ligament  breakage. 

This  volume  is  limited  to  the  analytical  methods,  theo¬ 
retical  aspects,  concepts  and  philosophy  for  advanced  durab¬ 
ility  analysis.  The  methodology  presented  in  this  report  is 
evaluated  and  justified  for  durability  analyses  in  Volume  II 
[3].  Three  other  volumes  for  the  program  are  as  follows: 
(1)  Volume  III  -  Fractographic  Test  Data  (4],  (2)  Volume  IV  - 
Executive  Summary  [5],  and  (3)  Volume  V  -  Durability  Analysis 
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Software  User's  Guida  [6].  Tha  Durability  Design  Handbook 
[7]  haa  also  baan  updated  for  tha  advancamanta  of  thia  pro¬ 
gram. 

1 . 2  BACKGROUND 

Tha  original  Air  Force  durability  daaign  raquiramanta 
[1]  ara  both  briaf  and  vagua.  Naw  raquiramanta  hava  racantly 
baan  propoaad  [8],  including  aaaumptiona  for  tha  initial  fa- 
tigua  quality  (IFQ)  of  a  faatanar  hola.  A  0.01"  initial 
quarter-circular  corner  crack,  asaumed  to  exist  in  a  faatanar 
hola,  ia  intandad  to  rapraaant  tha  IFQ  of  tha  atructura. 
However,  tha  aama  initial  crack  aiza  ia  uaad  to  rapraaant  tha 
IFQ  irraapactiva  of  tha  material  aalactad,  tha  hola  drilling 
procedurea  used,  and  tha  manufacturing  and  aaaambly  techni¬ 
ques  used.  A  more  realistic  representation  of  tha  IFQ  of 
structural  details  is  needed,  including  standard  methods  for 
defining  IFQ. 

A  durability  analysis  methodology  [2,9-17]  has  bean  de¬ 
veloped  for  satisfying  the  durability  design  requirements  of 
MIL-A-8866B.  This  methodology  accounts  for  IFQ  variation, 
fatigue  crack  growth  accumulation  in  a  population  of  struc¬ 
tural  details,  structural  properties  and  durability  analysis 
conditions  (e.g.,  load  spectra,  stress  level,  %  bolt  load 
transfer,  etc.).  A  statistical  distribution  of  equivalent 
initial  flaw  size  (EIFS)  is  used  to  represent  the  IFQ  of 
structural  details  (t  g.,  fastener  holes,  fillets,  cutouts, 
lugs,  etc.).  The  methodology  reflects  a  fracture  mechanics 
philosophy,  a  probabilistic  format  and  deterministic  crack 
growth.  The  equivalent  initial  flaw  size  distribution 
(EIFSD)  is  grown  forward,  using  the  applicable  durability  de¬ 
sign  conditions,  to  a  selected  service  time  using  a  determin¬ 
istic  crack  growth  approach.  Additional  references  on  the 
methodology  are  available  [18-26]. 
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The  prototype  durability  analysis  mathodology  uses  an 
"engineering  approach"  to  represent  the  XFQ  of  structural  de¬ 
tails.  XFQ  is  represented  by  an  equivalent  initial  flaw  size 
distribution  (EXFSD)  function.  Engineering  principles  rather 
than  mechanistic-based  theories  for  microstructural  behavior 
are  used.  This  approach  is  justified  for  durability  design 
applications. 

An  empirical-based  crack  growth  rate  model  and  fracto- 
graphic  results  are  used  to  define  the  EXFSD.  Equivalent  in¬ 
itial  flaws  are  determined  by  back-extrapolating  fractograph- 
ic  results  to  time  zero.  Such  flaws  are  artificial  rather 
than  actual  initial  flaws  in  a  fastener  hole,  such  as  a 
scratch,  a  burr,  or  a  void.  Since  EXFSs  cannot  be  verified 
by  NDI ,  indirect  methods  are  used  to  justify  their  use  for 
durability  analyses. 

The  methodology  uses  an  EXFSD  function  to  represent  the 
IFQ.  However,  the  parameters  for  describing  these  distribu¬ 
tions  have  not  been  optimized.  The  sensitivity  of  the  param¬ 
eters  with  respect  to  the  material,  geometry,  crack  size 
range,  maximum  stress  level,  amount  of  load  transfer,  and 
load  spectra  has  not  bean  established.  Such  an  optimization 
and  sensitivity  study  is  needed  to  adequately  represent  the 
initial  fatigue  quality  of  aircraft  structure. 

Once  the  EXFSD  has  been  determined,  a  service  crack 
growth  master  curve  is  used  to  predict  the  probability  of 
crack  exceedance,  p(i,T),  for  structural  details  at  future 
points  in  time.  However,  the  methods  for  predicting  p(i,T) 
have  not  been  optimized.  Such  an  optimization  study  is  need¬ 
ed  to  accurately  predict  the  probability  that  a  crack  will 
exceed  a  specific  size  at  the  point  in  time  of  interest. 

Fractographic  data  are  needed  to  determine  suitable 
EIFSDs  for  durability  analysis.  However,  design  guidelines 
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do  not  currently  exist  for  opacifying  tha  apaciman  geomet¬ 
ries,  number  of  holaa  par  specimen,  number  of  maximum  atraaa 
levels,  numbar  of  load  tranafar  1 avals ,  numbar  of  load  spec- 
tra ,  and  crack  aiza  rangaa  which  ara  naadad  to  adequately  re- 
praaant  tha  initial  fatigue  quality  of  aircraft  structure. 
Similarly,  guidalinaa  do  not  exist  for  making  accurate  p(i,7) 
pradlctiona. 

Tha  durability  analysis  methodology  was  based  on  a  dur¬ 
ability  definition  involving  tha  rework  of  fastener  holaa 
(i.e.,  tha  reaming  out  of  fastener  holes  to  tha  next  nominal 
hole  size) .  Such  a  definition  involves  small  crack  sizes 
(i.e.,  lass  than  0.1").  however,  durability  can  also  include 
definitions  of  functional  impairment  such  as  fuel  leakage  and 
ligament  breakage.  Such  definitions  normally  involve  crack 
sizes  greater  than  0.1".  Hence,  the  durability  analysis 
methodology  needs  to  be  extended  to  include  large  crack 
sizes.  Such  an  extension  for  large  crack  sizes  indicates 
that  some  of  the  small-crack-size  assumptions  used  for  the 
durability  analysis  methodology  are  invalid  and  must  be  re¬ 
placed  with  more  appropriate  formulations. 
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SECTION  II 


DURABILITY  ANALYSIS  METHODS  FOR  SMALL  CRACK 

SIZE  RANGE 

Advanced  durability  analysis  methodology  for  metallic 
airframes  is  described  and  discussed  in  this  section  for  ev¬ 
aluating  durability  requirements  for  functional  impairment 
due  to  excessive  cracking.  General  principles  and  procedures 
are  presented  for  application  in  the  small  crack  size  range. 
The  durability  analysis  extension  covering  the  large  crack 
size  range  (e.g.,  crack  sizes  associated  with  functional  im¬ 
pairment  due  to  fuel  leaks/ ligament  breakage)  is  presented  in 
Section  III. 

A  durability  analysis  methodology  for  the  small  crack 
size  region  (e.g.,  <0.10")  was  developed  under  the  previous 
program  [2,9-17],  Under  the  current  program  several  improve¬ 
ments  and  advancements  have  been  made  in  the  durability  anal¬ 
ysis  methodology  for  application  in  the  small  crack  size 
range.  In  particular,  improved  methods,  procedures  and 
guidelines  have  been  developed  for:  (1)  defining  a  suitable 
equivalent  initial  flaw  size  distribution  ,  and  (2)  making 
durability  analysis  predictions  for  the  probability  of  crack 
exceedance  and  cumulative  distribution  of  TTCI.  Significant 
improvements  and  refinements  have  been  made  in  the  methods 
for  defining  the  initial  fatigue  quality  of  structural  de¬ 
tails.  Methods  and  software  for  an  IBM  compatible  PC  [6] 
have  been  developed  for  estimating  and  optimizing  the  EIFSD 
parameters.  Technical  details  of  the  advancements  made  are 
described  and  discussed  in  this  Volume  (I) .  The  improved 
methods  and  procedures  are  evaluated  and  justified  in  Volume 
II  [3]. 

In  this  section  the  advanced  durability  analysis  method¬ 
ology  for  the  small  crack  size  range  is  summarized  and  de- 
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tails  ara  given  in  Appendices  A-G.  The  advanced  durability 
analysis  methodology  presented  in  this  section  is  formulated 
based  on  the  Weibull  compatible  EIFSD  function  and  the  deter¬ 
ministic  crack  growth  approach  (DCGA) .  Formulations  for  the 
following  EIFSD  functions  are  presented  in  Appendix  D:  two- 
parameter  Weibull,  lognormal  and  lognormal  compatible. 

Five  different  methods  were  investigated  under  this  pro¬ 
gram  for  estimating  the  EIFSD  parameters:  (1)  combined  least 
square  sums  approach  (CLSSA) ,  (2)  method  of  moments,  (3)  ho¬ 
mogeneous  EIFS  approach  (HEIFS) ,  (4)  upper  tail  fit  and  (5) 
nonlinear  approach.  Only  the  CLSSA  is  presented  in  this  sec¬ 
tion  and  the  other  approaches  are  described  and  discussed  in 
Appendix  F. 


2.1  DURABILITY  ANALYSIS  APPROACH 

The  durability  analysis  approach  for  the  small  crack 
size  range  (e.g.,  <0.10")  is  conceptually  described  in  Fig. 
1.  This  approach  reflects  a  probabilistic  format  with  a  de¬ 
terministic  crack  growth  law.  The  initial  fatigue  quality  of 
structural  details  is  represented  by  an  equivalent  initial 
flaw  size  distribution  (EIFSD) .  The  durability  analysis  ap¬ 
proach  accounts  for  initial  fatigue  quality,  crack  growth  ac¬ 
cumulation  in  a  population  of  structural  details,  load  spec¬ 
tra  and  structural  properties.  Essential  elements,  equations 
and  features  of  the  approach  are  described  in  the  following. 

The  Weibull  compatible  distribution  function  proposed 
by  Yang  and  Manning  [13,18]  has  been  found  to  be  reasonable 
for  representing  the  EIFS  cumulative  distribution  [2,13,18, 
24-30] . 
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Figure  1.  Deterministic  Crack  Growth  Approach  (DCGA). 
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in  which  ra(o)M  "  pCa(°)lxl»  a(0)  ■  EIFS  *  crack  size  at 
tine  t  *  0,  X  EIFS  upper  bound  linit  and  *•  and ^  are  empir¬ 
ical  paraneters. 

A  simple  deterministic  crack  growth  rate  model  for  the 
small  crack  size  region  designed  and  proposed  by  Yang  and 
Manning  [13,18],  is  used 

da(t)/dt  -  Q[a(t)]b  (2) 

where:  da(t)/dt  -  crack  growth  rate,  a(t)  -  crack  size  at 
any  time  t;  Q  and  b  are  empirical  crack  growth  rate  parame¬ 
ters.  Using  the  special  case  b  -  1  and  integrating  Eg.  2, 
one  obtains 

a(t)  -  a (0) exp (Qt)  (3) 

When  the  equivalent  initial  flaw  size  takes  a  sample  value 
x,  i.e.,  a(0)  ■  x,  it  will  grow  to  the  corresponding  sample 
value  Xj^  at  time  t,  i.e.,  a(t)  -  x.  It  follows  from  Eg.  3 
that 


Xj^  -  x  exp  (Qt) 


(4) 


Equation  4  is  referred  to  as  the  "EIFS  Master  Curve." 


The  cumulative  distribution  of  crack  size  a(t)  at  any 
time  t,  denoted  by  (*)  ■  pCa(t)<x],  can  be  obtained  from 
the  Weibull  compatible  distribution  of  a(0)  given  by  Eg.  1 
through  the  transformation  of  Eg.  3  as  follows  (Fig.  2) . 


p.(t) <x>  ■  •«> 


In  ( x  lx  )  +  Qt 
u _ 7 

4> 


1  :  *  >° 


and  t>0 


(5) 


The  probability  that  the  crack  size  a (T)  at  any  service 
time  will  exceed  any  given  crack  size  x^,  denoted  by  p(i,T)  » 
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P[a( T  )  >  x^],  is  referred  to  as  the  crack  exceedance  prob 

ability  (Fig.  1). 


p(i,T)  -  l-Fa(r)(x1)  -  1  -  exp  - 


ln(xjxx)+  Qfpj 


(6) 


If  FT(t)  denotes  the  cumulative  distribution  of  time  to 
initiate  a  given  crack  size  x^,  referred  to  as  time  to  crack 
initiation  (T  or  TTCI) ,  then  it  follows  that 


FrU)  -  P[T<t]  =1  -Fau|(x)  -  p(i,t) 


(7) 


*  1  -  exp 


In  ( x  /x.)  +  Qt 
u  i 

♦ 


Hence,  the  prediction  of  the  cumulative  distribution  of 
TTCI,  FT(t) ,  can  be  made  using  Eq.  7.  The  expression  for 
FT(t),  given  in  Eq.  7,  can  also  be  obtained  using  Eq.  1 
through  a  transformation  of  Eq.  4  as  conceptually  described 
in  Fig.  3. 


Equations  5-7  are  based  on  the  Weibull  compatible  EIFSD 
function  given  in  Eq.  1.  Similar  expressions,  based  on  the 
two-parameter  Weibull,  lognormal  and  lognormal  compatible 
EIFSD  functions,  are  given  in  Appendix  D. 


2.2  INITIAL  FATIGUE  QUALITY  REPRESENTATION 


The  initial  fatigue  quality  (IFQ)  defines  the  initially 
manufactured  state  of  a  structural  detail  or  details  with 
respect  to  initial  flaws  in  a  part,  component,  or  airframe 
prior  to  service.  An  equivalent  initial  flaw  size  distribu¬ 
tion  (EIFSD)  is  used  to  represent  the  IFQ  of  replicate  de¬ 
tails  (e.g.,  fastener  holes,  cutouts,  lugs,  fillets,  etc.). 

An  equivalent  initial  flaw  (EIFS)  is  an  artificial  crack 
size  which  results  in  an  actual  crack  size  at  an  actual  point 
in  time  when  the  initial  flaw  is  grown  forward.  It  is  deter¬ 
mined  by  back-extrapolating  fractographic  results.  It  has 
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Figure  3.  Transformation  of  the  Cumulative  Distribution  of  EIFS, 

F  .  . (x) ,  to  Obtain  the  Cumulative  Distribution  of  TTCI 


the  following  characteristics:  (1)  an  EIFS  is  an  artificial 
crack  assumed  to  represent  the  initial  fatigue  quality  of  a 
structural  detail  in  the  as-manufactured  condition  whatever 
the  source  of  fatigue  cracking  may  be,  (2)  it  has  no  direct 
relationship  to  actual  initial  flaws  in  fastener  holes  such 
as  scratches,  burrs,  microdefects,  etc.,  and  it  cannot  be 
verified  by  NDI ,  (3)  it  has  a  universal  crack  shape  in  which 
the  crack  size  is  measured  in  the  direction  of  crack  propaga¬ 
tion,  (4)  EIFSs  are  in  a  fracture  mechanics  format  but  they 
are  not  based  on  linear  elastic  fracture  mechanics  (LEFM) 
principles  and  are  not  subject  to  such  laws  or  limitation- 
such  as  the  "short  crack  effect,"  (5)  it  depends  on  the 
fractographic  data  used,  the  fractographic  crack  size  range 
used  for  the  back-extrapolation  and  the  crack  growth  rate 
model  used,  (6)  it  must  be  grown  forward  in  a  manner  consis¬ 
tent  with  the  basis  for  the  EIFS,  (7)  EIFSs  are  not  unique  - 
a  different  set  is  obtained  for  each  crack  growth  law  used 
for  the  back-extrapolation , and  (8)  EIFSs  are  not  necessarily 
comparable  with  minimum  crack  sizes  obtained  from  fracte- 
graphic  results  from  tear-down  inspections. 

An  EIFSD  is  an  "engineering  approach"  for  representing 
IFQ.  It  is  not  based  on  mechanistic  theories  for  microstruc- 
tural  behavior.  This  approach  for  representing  the  IFQ,  when 
properly  used,  can  provide  reasonable  durability  analysis  re¬ 
sults  for  metallic  airframes. 

Since  EIFSs  are  determined  by  back-extrapolating  fracto¬ 
graphic  results,  they  are  not  strictly  "generic,"  because  the 
fractographic  results  are  obtained  from  fatigue  cracks  gener¬ 
ated  under  specified  test  conditions  (e.g.,  load  spectrum, 
stress  level,  %  bolt  load  transfer,  fastener  hole  type/fit., 
etc.).  Therefore,  the  EIFSs  depend  on  the  test  conditions  re¬ 
flected.  The  real  issue  of  the  EIFS  approach  is  not  whether 
the  EIFSs  or  EIFSDs  are  "generic"  or  not  but  rather  can  reas¬ 
onable  durability  analysis  predictions  for  p(i,T)  and/or 
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FT(t)  be  obtained  for  a  given  set  of  conditions  using  EXFSs 
based  on  another  set  of  conditions  (e.g.,  different  load  spec 
truxn,  stress  level,  %  bolt  load  transfer,  etc.)?  Under  this 
program  we  have  shown  that  fractographic  data  pooling  and 
statistical  scaling  techniques  can  be  effectively  used  to  de¬ 
velop  meaningful  ElFSDs  for  practical  durability  analyses. 
Hence,  the  EIFSD  can  be  established  for  a  given  set  of  condi¬ 
tions  and  the  same  EIFSD  can  be  used  to  make  reasonable  pre¬ 
dictions  for  another  set  of  conditions.  Data  pooling  and 
statistical  scaling  techniques  are  described  and  discussed  in 
Sections  2.3.1  and  2.3.2,  respectively. 

Durability  analysis  concepts,  issues  and  philosophy  are 
discussed  in  detail  in  Appendix  A.  This  includes  the  current 
understanding  of  EIFS  concepts,  principles  and  practices. 

2.3  ESTIMATION  OF  EIFSD  PARAMETERS 

The  following  Weibull  compatible  EIFSD  parameters  in  Eq. 
1  must  be  estimated  and  optimized  to  define  IFQ:  xu ,  and  <p . 
In  this  section  methods  are  presented  for  estimating  these 
parameters,  including  procedures  for  pooling  fractographic 
results,  statistical  scaling,  optimization  and  checking  suit¬ 
ability  of  the  candidate  EIFSD  for  durability  analysis. 

The  methodology  for  determining  a  suitable  EIFSD  for 
durability  analysis  includes  several  elements:  (1)  selection 
and  use  of  fractographic  data,  (2)  fractographic  crack  size 
range,  AL-AU,  and  reference  crack  size,  xx  for  time-to-crack- 
initiation  (TTCI) ,  (3)  deterministic  crack  growth  rate  model, 
(4)  crack  size-time  relationship  or  "EIFS  Master  Curve",  (5) 
EIFSD  function  (e.g.,  Weibull  compatible  or  lognormal  compat¬ 
ible)  ,  (6)  determination  of  EIFSs  and/or  TTCIs  ,  (7)  statis¬ 
tical  scaling  procedure  for  obtaining  homogenous  data,  (8) 
estimation  of  EIFSD  parameters  based  on  data  pooling  proce¬ 
dure  (i.e.,  combined  least  square  sums  approach  (CLSSA)  or 
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method  of  moments  approach),  (9)  iterative  procedure  for  op¬ 
timizing  the  EZFSD  parameters  and  (10)  justification  of  can¬ 
didate  EIFSD  for  durability  analysis. 

Detail  procedures  and  methods  for  estimating  and  opti¬ 
mizing  the  EIFSD  parameter  are  given  in  the  following  subsec¬ 
tions  and  in  Appendices  B,  C,  E,  F  and  G.  The  general  proce¬ 
dure  for  defining  IFQ  is  summarized  below  and  key  elements 
are  described  in  Fig.  4. 

1.  Select  suitable  EIFSD  function  for  representing  the 
initial  fatigue  quality  (e.g.,  Weibull  compatible  or  Lognor¬ 
mal  compatible) . 

2.  Select  fractographic  data  set(s)  to  be  used  to  de¬ 
termine  the  EIFSD.  The  data  sets  should  be  for  the  same 
material,  same  type  load  spectrum  (e.g.,  fighter,  bomber  or 
transport)  and  type  fastener/  hole/fit  (i.e.,  straight  bore 
or  countersunk) . 

3.  Select  fractographic  crack  size  range,  AL-AU,  to  be 
used  to  determine  the  crack  size-time  relationship  from  a 
suitable  deterministic  crack  growth  rate  model.  Also  select 
a  reference  crack  size,  x1#  for  the  TTCI's  for  each  fracto¬ 
graphic  data  set. 

4.  Estimate  the  crack  growth  parameter,  Q^,  or  "pooled 
Q" ,  value  for  each  fractographic  data  set  using  the  fracto¬ 
graphic  data  in  the  AL-AU  crack  size  range  (see  Appendix  E) . 
A  DCGA-based  crack  size-time  relationship  or  EIFS  master 
curve  is  defined  in  terms  of  Q^.  This  relationship  defines 
the  equivalent  initial  flaw  size  (EIFS)  corresponding  to  a 
given  crack  size  x1  at  time  T  or  vice  versa. 

5.  Compute  TTCI's  for  each  data  set(s)  using  the  re¬ 
ference  crack  size,  x^. 
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— -ltial  Steps  in  Procedure  Leading  to  Estimation  of 
EIFSD  Parameters. 


6.  Select  either  the  "TTCI  fit"  or  the  "EIFS  fit"  ap¬ 
proach  to  fit  the  EIFSD  parameters .  Refer  to  Appendix  F. 

7.  If  the  EIFSe  are  to  be  used  to  fit  the  EIFSD  para¬ 
meters,  compute  the  EIFSs  for  each  fractographic  data  set 
using  the  TTCIs  from  step  5  and  the  EXFS  master  curve  from 
step  4.  Each  EIFS  in  a  given  data  set  is  determined  using 
the  same  crack  size-time  relationship  or  EIFS  master  curve. 
Hence,  each  EIFS  in  the  data  set  has  a  corresponding  TTCI  and 
vice  versa. 


8.  Select  the  data  pooling-parameter  fitting  method  to 
be  used  to  estimate  the  EIFSD  parameters:  (1)  combined  least 
square  sums  approach  (CLSSA)  or  (2)  method  of  moments  ap¬ 
proach.  Both  approaches  include  a  scaling  procedure  for  mak¬ 
ing  the  fractographic  results  homogeneous  for  different  data 
sets,  with  scaling,  results  for  different  fractographic  data 
sets  can  be  "mixed  and  matched"  to  estimate  the  EIFSD  parame¬ 
ters  in  a  global  sense. 

9.  Optimize  the  EIFSD  parameters  using  an  iterative 
procedure.  Details  are  given  in  Section  2.3.4. 

10.  Justify  using  the  EIFSD  for  durability  analysis. 
Correlate  predictions  for  the  cumulative  distribution  of  TTCI 
and/or  the  cumulative  distribution  of  crack  size  with  fracto¬ 
graphic  results  for  selected  data  sets.  The  goodness-of-fit, 
based  on  the  proposed  EIFSD,  should  be  checked  using,  as  a 
minimum,  data  sets  used  to  estimate  the  EIFSD.  We  recommend 
that  other  fractographic  data  sets  also  be  used  to  check  the 
"goodness-of-fit. " 

Computer  software  is  available  in  Volume  V  [6]  to  esti¬ 
mate,  optimize  and  justify  the  EIFSD  parameters  for  durabil¬ 
ity  analysis,  including  a  goodness-of-fit  plotting  capabil¬ 
ity. 
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2.3.1  Data  Pooling 


Data  pooling  la  a  practical  way  to  define  the  EIFSD  pa¬ 
rameters  for  diffarant  data  aata  (e.g.,  a am a  typa  of  load 
spectrum  but  diffarant  atraaa  lavala  and  %  bolt  lead  trans- 
far) .  Pooling  affactivaly  incraaaaa  tha  samp la  size  and  con- 
fidanca  in  tha  EIFSD  parameters.  Likewise,  thia  ia  a  reason- 
ablo  approach  for  justifying  an  EIFSD  for  more  ganaral  appli¬ 
cations.  Thia  is  an  important  parapactiva  for  practical  dur¬ 
ability  analysis. 

Tha  fractographic  data  pooling  concapt  io  conceptually 
illustrated  in  Fig.  5  using  two  fractographic  data  seta.  Tha 
applicable  TTCI  values  for  each  data  aet  are  represented  by  a 
TTCID  as  shown  in  Fig.  5.  Tha  basic  idea  of  data  pooling  is 
that  tha  crack  growth  results  for  each  data  set  have  a  common 
EIFSD.  For  example,  if  the  TTCIs  for  each  data  set  are  re¬ 
gressed  backwards  to  time  zero  using  the  applicable  EIFS  mas¬ 
ter  curve  for  each  data  set,  the  resulting  EIFSs  have  the 
same  EIFSD.  Here's  another  way  tc  look  at  this.  The  TTCID 
for  different  fractographic  data  sets  can  be  determined  using 
the  same  EIFSD. 

The  EIFSD  parameters  can  be  estimated  using  either  an 
"EIFS  fit"  or  a  "TTCI  fit."  Either  approach  will  give  the 
same  EIFSD  parameters.  There  are  subtle  differences  between 
the  two  aproaches. 

With  the  "EIFS  fit,"  TTCIs  for  a  given  reference  crack 
si z®,  x1#  are  back-extrapolated  to  time  zero  using  the  appli¬ 
cable  EIFS  master  curve  (e.g.,  Eg.  4)  to  obtain  the  EIFSs  for 
each  data  set.  The  resulting  EIFSs  and  a  "fitting  procedure" 
are  then  used  to  estimate  the  EIFSD  parameters  in  a  global 
sense  for  the  particular  EIFSD  function  used  (see  Fig.  6(a)), 
i.e.,  Eg.  1.  Different  "fitting  procedures"  have  been  inves- 
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tigatsd  but  only  ths  combined  least  square  turns  approach  will 
ba  oonsldarad  In  Saotlon  ZI  (raf.  Saotlon  2.3.3).  othar 
"fitting  procaduras"  ara  daacribad  and  diaoussad  in  Appendix 
F. 


A  ganaral  axpraaaion  for  tha  cumulative  distribution  of 
TTCX ,  FT(t),  is  givan  in  Eq,  ?.  Thia  axpraaaion  is  baaad  on 
a  transformation  of  tha  Waibull  compatible  EIFSD  using  tha 
EXFS  master  curve  givan  in  Eq.  4  (sea  Fig.  3) .  Equation  7 
contains  all  tha  EXFSD  parameters  in  Eq.  l  (i.e.,  xu,<*and£) 
plus  tha  crack  growth  rata  parameter  Q.  Therefore,  with  tha 
"TTCX  fit"  tha  EXFSD  parameters  can  ba  estimated  using  TTCX 
values  for  each  data  sat  considered  in  tha  data  pooling  pro¬ 
cedure  (Fig.  6(b)).  Q  can  be  obtained  from  fractographic  re¬ 
sults  and/or  using  a  suitable  analytical  crack  growth  program 
[e.g.,  38,  39].  Detailed  procedures  for  computing  Q  are  given 
in  Appendix  E. 

With  the  "TTCI  fit"  the  EIFSD  parameters  can  be  estimat¬ 
ed  without  computing  EXFSs.  Since  TTCIa  are  actual  physical 
quantities  and  EIFSs  ara  not,  some  may  prefer  to  quantify  the 
EIFSD  parameters  using  TTCX  values.  The  same  Q  value  is  used 
for  both  the  "EIFS  fit"  and  the  "TTCI  fit."  In  any  case,  it 
can  be  shown  that  the  same  EXFSD  parameters  are  obtained  from 
a  given  xu  value  using  either  approach. 

2.3.2  Statistical  Scaling 

The  IFQ  or  EIFSD  for  fastener  holes  is  defined  for  an 
"equivalent  single  hole  population."  Therefore,  the  fatigue 
cracking  resistance  of  each  fastener  hole  in  each  test  speci¬ 
men  is  accounted  for  in  the  definition  of  the  EXFSD.  Test 
specimens  for  acquiring  fatigue  crack  growth  data  may  have 
one  or  more  fastener  holes  per  specimen.  Some  specimens  may 
or  may  not  be  fatigue  tested  to  failure.  Also,  every  fasten¬ 
er  hole  in  each  replicate  test  specimen  may  not  contain  a 
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measurable  fatigue  crack  or  else  tha  crack  ia  too  small  or 
complex  (a.g.,  multiple  crack  origins  and  branching)  for 
fractographic  analysis.  A  statistical  scaling  tachniqua  has 
baan  davalopad  (saa  Appandix  B  for  details)  for  determining 
tha  EZFSD  for  an  '’equivalent  single  hole  population"  based  on 
tha  largest  fatigue  crack  par  specimen.  Hence,  it  is  neces¬ 
sary  to  read  only  the  fractographic  results  for  the  largest 
crack  par  specimen.  Essential  elements  are  conceptually  des¬ 
cribed  in  Fig.  B.l.  This  technique  is  very  general  and  is 
independent  of  the  distribution  functions  used.  It  accounts 
for  the  number  of  fastener  holes  per  test  specimen  in  a  given 
fractographic  data  set.  It  minimizes  the  fractographic  read¬ 
ing  requirements,  permits  a  maximum  utilization  of  the  avail¬ 
able  fractographic  data  and  allows  for  "mixing  and  matching" 
of  fractographic  data  for  the  largest  crack  in  specimens  with 
a  different  number  of  holes. 


Details  of  the  statistical  scaling  technique  developed 
are  given  in  Appendix  B.  Essential  features  and  key  equa¬ 
tions  are  summarized  in  the  following. 


Let  the  cumulative  distribution  of  EIFS  for  a  single 
hole  population  be  denoted  by  ya(o) and  that  of  the  maxi¬ 
mum  EIFS  based  on  the  largest  fatigue  crack  per  specimen  with 
£  fastener  holes,  be  denoted  by  Fa  ^ (x) .  Assuming  that  fa¬ 
tigue  cracking  in  each  fastener  hot;  of  a  specimen  is  statis¬ 


tically 
t0  Fa(0) 


independent  of  the  other  holes,  Fa^0) (x) 
(x)  through  the  following,  * 


is  related 


FS(0){X) 


(8) 


where  X  «  number  of  fastener  holes  per  specimen,  similar  ex¬ 
pression  for  the  cumulative  distribution  of  TTCI  are  given  in 
Eqs.  9  and  10. 
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Fjt)  *  1  -  1  -  F  U) 


ut 


F  «)  8  l_ 
1 1 


1  -FT(t) 


(9) 

(10) 


where  FT(t)  ■  cumulative  distribution  of  TTCI  for  an  equiva* 
lent  single  hole  population,  and  F-  (t)  -  cumulative  distri* 
bution  of  the  minimum  TTCI  in  Wholes. 


In  a  similar  manner,  let  Fft^j  (x)  denote  the  cumulative 

distribution  of  crack  size  at  any  time  t  for  an  "equivalent 

* 

single  hole  population"  for  Fa(0)(x).  Let  Fft  (x)  denote 
the  cumulative  distribution  of  crack  size  at  an§  t  based  on 


the  EIFS  cumulative  distribution,  F 
is  related  to  F,  ... (x)  as  follows. 


atf(0) 


(x) .  Then,  Fft(t)(x) 


Fc(U>ix)  " 


(ID 


Let  FT(t)  denote  the  cumulative  distribution  of  TTCI  for 
an  "equivalent  single  hole  population"  and  FT  (t)  denote  the 
cumulative  distribution  of  TTCI  based  on  &  smallest  TTCI 
per  specimen  with  JL  fastener  holes.  Then,  the  following  re¬ 
lationship  between  F-(t)  and  F_  (t)  holds. 

1  X  0 


Fr  (<> 
( 


=  1-  [  1  -FT(t) 


(12) 


The  simple  scaling  technique  described  in  this  section 
has  been  incorporated  into  the  procedure  for  estimating  the 
EIFSD  parameters  for  one  or  more  fractographic  data  sets. 
Computer  software  for  the  IBM  compatible  PC  is  available  for 
estimating  the  EIFSD  parameters,  for  checking  goodness-of-fit 
and  for  plotting  results.  The  statistical  scaling  technique 
is  evaluated  in  Volume  II  [ 3 ] . 
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2.3.3  Combined  Least  Square  Sums  Approach  (CLSSA) 


The  data  pooling  and  statistical  scaling  concepts  des¬ 
cribed  in  Sections  2.3.1  and  2.3.2,  respectively,  are  impor¬ 
tant  considerations  for  estimating  the  EIFSD  parameters  for 
fastener  holes  or  other  types  of  structural  details.  Data 
pooling  is  important  for  two  reasons:  (1)  the  sample  size 
is  increased  and  so  is  the  confidence  in  the  estimated  EIFSD 
paameters  and  (2)  the  derived  EIFSD  can  be  justified  for  a 
broader  range  of  variables  (e.g.,  stress  level,  %  bolt  load 
transfer,  etc.)  than  that  for  a  single  fractographic  data 
set. 


The  desired  EIFSD  should  account  for  the  fatigue  crack¬ 
ing  resistance  of  each  fastener  hole  in  each  test  specimen 
for  each  fractographic  data  set  to  obtain  a  "single  hole  pop¬ 
ulation  basis."  Test  specimens  for  one  fractographic  data 
set  may  have  a  different  number  of  fastener  holes  than  an¬ 
other  data  set.  A  method  is  needed  for  "mixing  and  matching" 
or  pooling  the  fractographic  results  for  different  data  sets 
so  that  the  EIFSD  parameter  (e.g.,  xu,  oc  and  in  Eq.  1)  can 
be  estimated  in  a  global  sense.  A  statistical  scaling  tech¬ 
nique  is  described  in  Section  2.3.2  and  Appendix  B  for  scal¬ 
ing  the  fractographic  results  of  test  specimens  with  one  or 
more  fastener  holes  to  an  "equivalent  single  hole  population" 
basis. 

In  this  section,  a  method,  referred  to  as  the  "combined 
least  square  sums  approach"  (CLSSA) ,  is  presented  for  esti¬ 
mating  EIFSD  parameters  for  one  or  more  fractographic  data 
sets.  This  data  pooling  method  accounts  for  statistical 
scaling.  In  this  section  the  CLSSA  is  described  and  discuss¬ 
ed  for  the  Weibull  compatible  EIFSD  function.  Other  EIFSD 
functions  and  methods  are  described  in  Appendix  F. 
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The  EIFSD  parameters  for  tha  Waibull  compatible  EIFSD 
function  (i.e.,  xu,  otand^)  can  ba  estimated  using  aithar 
TTCI  or  EIFS  input  data  for  ona  or  mora  fractographic  data 
sats.  Either  approach  will  give  tha  same  EIFSD  parameter 
values.  The  TTCIs  or  EIFSs  (whichever  is  used)  for  each 
fractographic  data  set  are  ranked  and  least  square  sums  are 
obtained  for  each  data  set  using  the  applicable  EIFSD  func- 
tion  and  data  used  to  fit  the  EIFSD  parameters.  Each  least 
square  sum  accounts  for  the  scaling  factor  for  each  data  set. 
The  combined  least  square  sums  for  all  data  sets  are  obtained 
by  summing  the  applicable  value  for  each  data  set.  Then,  the 
EIFSD  parameters  are  estimated  using  the  combined  least 
square  sums  and  a  global  fit. 

The  procedure  for  implementing  the  CLSSA  is  given  below: 

1.  Select  fractographic  data  sets  to  be  used  in  the 
pooling  procedure  to  estimate  the  EIFSD  parameters  in  a  glo¬ 
bal  sense. 

2.  Acquire  the  TTCIs  or  EIFSs  for  each  data  set  to  be 
used  to  estimate  the  EIFSD  parameters. 

3.  Define  scaling  factor,  ^  ,  for  each  data  set  to  be 

used. 

4.  Select  the  EIFSD  function  to  be  used  and  then  trans¬ 
form  it  into  a  linear  least  square  fit  form.  If  EIFSs  are 
used  for  the  fit,  use  the  expression  for  the  cumulative  dis¬ 
tribution  of  EIFS  for  the  fit.  On  the  other  hand,  if  TTCIs 
are  used  for  the  fit,  use  the  expression  for  the  cumulative 
distribution  of  TTCI  based  on  the  transformation  of  the  EIFSD 
and  the  crack  growth  model  used  (e.g.,  Eq.  2). 

5.  Assume  a  value  for  the  EIFS  upper  bound  within  the 
allowable  limits  for  fastener  holes:  largest  EIFS  in  any 
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data  sat  <xu<  0.05".  Compute  tha  combined  least  square  sums 
for  all  data  sets  combined.  Use  these  results  to  estimate 
the  EIFSD  parameters. 

6.  The  EIFSD  parameters  can  be  optimized  within  the  xu 
constraints.  A  reasonable  way  to  do  this  is  to  minimize  the 
total  standard  error  for  all  data  sets  used  to  determine  the 
EIFSD  parameters.  The  optimum  xu  can  be  determined  by  trial 
and  error  until  the  total  standard  error  has  been  minimized. 


The  CLSSA  is  described  in  detail  in  Section  F.2.  Essen¬ 
tial  elements  and  equations,  based  on  the  Weibull  compatible 
EIFSD,  are  summarized  in  the  following  for  both  the  "TTCI 
fit"  and  "EIFS  fit"  approaches. 


It  has  been  demonstrated  in  Volume  II  [3]  that  the  "TTCI 
fit"  and  the  "EIFS  fit"  will  give  the  same  EIFSD  parameter 
values.  From  a  philosophical  standpoint,  the  TTCI  fit  is  re¬ 
commended  because  TTCIs  can  be  directly  verified  from  the 
fractographic  data,  and  EIFSs  do  not  have  to  be  computed; 
whereas,  the  EIFSs  are  artificial  quantities  that  cannot  be 
verified  by  NDI .  If  the  EIFS  fit  is  used  the  resulting  EIFSD 
has  to  be  grown  forward  anyway  to  test  the  goodness-of-f it  of 
the  Fa£tj(x)  and  FT ( t)  predictions.  It's  all  a  matter  of 
choice  -  either  method  is  acceptable.  Software  is  available 
in  Volume  V  [6]  for  an  IBM  compatible  PC  to  implement  either 
method  for  the  Weibull  compatible  EIFSD  function. 


2. 3. 3.1  TTCI  Fit. 


The  time-to-crack  initiation  (TTCI)  is  defined  as  the 
time  required  to  initiate  a  specified  crack  size,  x^,  in  a 
structural  detail  in  the  as-manufactured  condition.  For  ex¬ 
ample,  in  Fig.  7  fractographic  results  are  depicted  for  three 
fatigue  cracks.  The  time  to  initiate  a  crack  size,  x^,  for 
each  of  these  tracks  is  t^,  t2  and  t3,  respectively,  as  shown 

in  Fig.  7. 
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Figure  7.  Tlme-To-Crack-Inltlatlon  (TTCI)  Concept. 
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TTCla  can  ba  determined  froa  tha  fractographic  raaulta 
in  ona  of  tha  following  ways:  (1)  diractly  froa  tha  fracto- 
graphic  raaulta  when  tha  maaaurad  crack  aiza  ia  axactly  equal 
to  (2)  by  interpolation  (e.g.,  thraa-point  Lagrangian) , 
and  (3)  by  extrapolation  (saa  Section  G.3  for  dataila) . 


The  Weibull  compatible  E1FSD  parameters  oC  and  ^  can  ba 
determined  for  a  given  xu  uaing  Eq.  13  and  14,  respectively. 
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The  terms  in  Eqs.  13  and  14  are  defined  as  shown  in  Eq.  15. 
The  expression  for  the  total  standard  error  is  given  in  Eq. 
16. 
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In  Eqs.  13  through  16 ,  -  numb or  of  fastener  holaa  for  tha 

raplicata  specimen  in  data  aat  i,  M  *  number  of  fractographlc 
data  sats  usad  in  tha  pooling  procedure,  -  numbar  of  TTCI 
valuaa  for  tha  ith  data  sets,  xy  «  EIFS  uppar  bound  limit,  x 
-  aQ  -  rafaranca  crack  siza  for  TTCls,  «  poolad  Q  valua 
for  tha  ith  data  aat  (saa  Eq.  2),  t^  -  jth  TTCI  valua  for 
tha  ith  data  sat,  F_  (t,,)  -  jth  valua  for  tha  cumulativa 

*  8*i  *  J 

distribution  of  TTCI  minimums  (i.a.,  based  on  tha  largast 
fatigua  crack  par  spaciman  in  aach  data  sat)  for  tha  ith  data 

sat.  Rafar  to  Appandix  F  for  datails. 

Tha  Waibull  compatible  parameters  o t  and  $  for  a  given  xu 
can  be  determined  using  tha  same  aquations  usad  for  tha  TTCI 
fit  approach  (i.a.,  Eqs.  13  and  14).  For  tha  EIFS  fit  the 
expressions  shown  in  Eq.  17  should  be  used  in  Eq.  13  and  14. 
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The  expression  for  tha  total  standard  error,  TSE,  is 
given  in  Eq.  18 . 
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In  Eqs.  17  and  18  all  terms  are  defined  tha  same  way  as 
those  for  Eq.  15  except  the  following:  »  number  of  EIFS 

values  for  the  ith  data  set,  x..»  jth  EIFS  value  for  the  ith 

■i 

data  set,  and  F  “jth  valua  for  tha  cumulative 

distribution  of  EIFS  maximums  for  tha  ith  data  sat. 


2. 3. 3. 2  EIFS  Fit,  Methods  for  determining  EIFSe  ueing 
the  DCGA  are  described  in  Appendix  E.  The  recommended  proce¬ 
dures  for  computing  DCGA-based  EIFSs  are  summarized  below: 

1.  Select  the  fractographic  data  to  be  used. 

2.  Select  the  fractographic  crack  size  range,  AL-AU, 
to  be  used  to  compute  the  "pooled  Q"  value  (e.g.,  0.01"-0.05" 
is  considered  reasonable  for  defining  IFQ  for  fastener 
holes)  . 


3.  Select  a  reference  crack  size,  x^,  within  the  fol¬ 
lowing  constraint:  AL  <xx<AU.  Then,  determine  the  TTCIs  for 
the  largest  fatigue  crack  for  each  specimen  in  the  data  set. 

4.  Compute  the  "pooled  Q"  value  for  the  given  fracto¬ 
graphic  data  set  using  Eg.  C-13.  Q  is  based  on  fractographic 
results  for  the  selected  AL-AU  range.  The  pooled  Q  value  is 
used  in  Eq.  4  to  define  the  crack  size-time  relationship  for 
the  data  set.  This  relationship  is  used  as  an  "EIFS  master 
curve . " 


5.  Use  Eg.  4  to  determine  the  EIFS  corresponding  to 
each  TTCI  in  the  data  set  as  conceptually  described  in  Fig. 
8. 


The  EIFSs  for  one  or  more  data  sets  can  be  used  to  fit 
the  EIFSD  parameters  for  the  Weibull  compatible  EIFSD  func¬ 
tion.  Methods  and  equations  are  given  in  detail  in  Appendix 
F. 


2.3.4  Optimization  of  Parameters  and  Goodness-of-Fit 

The  purpose  of  this  section  is  to:  (1)  describe  an  it¬ 
erative  procedure  for  optimizing  the  EIFSD  parameters  for 
structural  details  and  (2)  describe  a  way  to  check  the  good- 
ness-of-fit  of  the  candidate  EIFSD  and  justify  it  for  durab¬ 
ility  analysis. 
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Not* :  Datarmina  EXFSa  corraaponding  to  TTCXa  for  givan 
rafaranc*  crack  aiza.  Fit  tha  EIFSD  parametara 
using  EXFSa. 


xi  ■  *x  •*p(-QTi  ) 


EIFSD 


TIME 


Figure  8.  EIFS  Approach  for  Estimating  EIFSD  Parameters. 


A  general  procedure  for  optimizing  th®  EZFSD  parameters 
and  checking  goodn®®a-of-fit  for  compatible  typ®  EIFSD  func¬ 
tion®  is  shown  in  Fig.  9.  Th®  following  it«rativ®  approach 
is  consid®r®d  raasonabl®  for  optimizing  th®  EIFSD  paramatar®. 
An  upp®r  bound  EIFS  limit,  xu,  is  assumad  with  th®  following 
limits  for  fastanar  holes i  largast  EIFS  in  any  data  sat 
<xu<0.05".  If  tha  "TTCI  fit"  is  usad,  the  smallest  TTCI  for 
each  data  sat  can  be  converted  to  an  EIFS  using  tha  EIFS-time 
relationship  of  Eq.  4.  The  parameters  U  and  ^  in  tha  Weibull 
compatible  EIFSD  function,  Eq.  1,  can  be  estimated  using  the 
CLSSA  described  in  Section  2.3.3.  Then,  the  total  sum  squar¬ 
ed  error  can  be  computed  using  either  the  TTCls  or  EIFSs  for 

those  data  sets  used  to  estimate  ot  and  ^  .  The  optimum  oc  and 
4>  for  a  given  xu  is  obtained  when  the  total  standard  error, 
TSE,  for  all  data  sets  used  in  the  pooling  procedure  has  been 
minimized  within  the  allowable  limits  for  xu  (i.e.,  largest 
EIFS  in  any  data  set  <xu<0.05"  for  fastener  holes) .  An  iter¬ 
ative  procedure  is  used  to  optimize  xu,  it  and  <P  .  Software  is 
available  in  Volume  V  [6]  to  efficiently  optimize  the  EIFSD 
parameters  using  the  procedure  described.  After  estimating 
the  EIFSD  parameters,  the  EIFSD  needs  to  be  tested  to  deter¬ 
mine  if  it  is  suitable  for  representing  the  initial  fatigue 
quality  of  the  structural  details  in  question. 

The  procedure  used  to  estimate  the  EISD  parameters  ac¬ 
counts  for  the  data  fit.  For  example,  suppose  the  CLSSA  and 
a  "TTCI  fit"  is  used  to  estimate  the  EIFSD  parameters.  In 
this  case,  the  EIFSD  parameters  are  optimized  in  a  least 
square  sense  using  the  cumulative  distribution  of  TTCI.  Fur¬ 
ther  testing  of  the  resulting  EIFSD  is  recommended,  including 
goodness-of-f it  plots. 

The  EIFSD  is  determined  using  the  fractographic  results 
for  one  or  more  data  sets.  Each  data  set  is  assumed  to  have 
the  same  type  spectrum  (e.g.,  fighter,  bomber  or  transport), 
fastener/hole/fit  and  material.  Stress  level,  %  bolt  load 
transfer,  and  number  of  fastener  holes  can  vary  for  each 
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Figure  9.  General  Procedure  for  Optimizing  EIFSD  Paramati 
and  Checking  Goodnesa-of-Flt  for  Compatible  T5 
EIFSD  Function. 
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fractographic  data  aat  to  ba  uaad  to  estimate  tha  EZFSO  pa- 
raiaatara  in  a  global  aanaa.  In  most  cases,  tha  durability 
daaign  conditions  will  not  likaly  ba  tha  sama  as  thosa  re- 
flactad  in  any  ona  of  tha  fractographic  data  sats  usad  to 
astimata  tha  EIFSD  paramatars.  Tha  EIFSO  is  justifiad  for 
other  design  conditions  by  showing  that  "reasonable"  durab¬ 
ility  analysis  results  can  ba  obtained  for  thosa  fractograph¬ 
ic  data  sats  usad  to  define  tha  IFQ  and  for  selected  data 
sets  not  usad  to  define  tha  IFQ.  In  such  case,  it  is  assumed 
that  tha  candidate  EIFSD  will  ba  applied  to  tha  same  or  simi¬ 
lar  material,  sama  type  load  spectrum  and  same  type  fastener/ 
hole/fit.  For  durability  analysis  applications,  stress 
level,  %  bolt  load  transfer  and  number  of  fastener  holes  can 
vary. 


The  EIFSD  can  be  justified  using  the  fractographic  data 
used  to  estimate  the  EIFSD  parameters  and/or  other  fracto¬ 
graphic  data  sets.  For  example,  tha  "pooled  Q"  value  for 
each  fractographic  data  set  to  be  considered  and  the  candi¬ 
date  EIFSD  can  be  used  to  predict  the  cumulative  distribution 


of  crack  size  at  any  time  t, 
tribution  of  TTCI  at  any  crack  size, 


ra(t) <x) ' 


or  the  cumulative  dis- 


FT(t) 


The  predict- 

k 

ed  cumulative  distributions  can  ba  plotted  and  the  results 
compared  with  the  ranked  fractographic  observations  for 
either  crack  size  or  TTCI.  Such  plots  are  conceptually  de¬ 
picted  for  the  probability  of  crack  exceedance  and  cumulative 
distribution  of  TTCI  in  Figs.  10  and  11,  respectively. 


The  goodness-of-f it  in  the  tail  of  the  applicable  cumu¬ 
lative  distribution  is  of  most  interest.  For  example,  the 
smaller  TTCIs  for  a  given  crack  size,  x1#  and  the  largest 
crack  sizes  for  a  given  time,  t,  are  of  most  interest  for 
judging  goodness-of-f it. 


Different  criteria  could  be  used  to  judge  goodness-of- 
fit.  For  example,  a  visual  inspection  of  the  plotted  results 
or  the  total  standard  error  might  be  used. 
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2.4  DURABILITY  ANALYSIS  PREDICTIONS 


Durability  anaysis  pradictiona  for  tha  probability  of 
crack  exceedance,  p(i,T),  and  tha  cumulative  distribution  of 
TTCI ,  PT(t),  can  ba  determined  using  an  appropriata  EIFSD 
function,  p  <x>,  and  a  suitabla  sarvica  crack  growth  mas- 
tar  curva  (SCGMC) .  Guidalinas  and  procaduras  ara  prasantad 
in  this  saction  for  datarmining  tha  SCGMC  and  for  making  dur¬ 
ability  analysis  pradictions  for  p(i,T)  and  FT(t)  in  tha 
small  crack  siza  ranga  (a.g.,  <0.10"). 

2.4.1  Sarvica  Crack  Growth  Mastar  Curva 

A  Sarvica  Crack  Growth  Mastar  Curva  (SCGMC)  dascribas 
tha  datarministic  crack  siza-tima  ralationship  for  a  givan 
stress  region.  It  is  used  to  grow  tha  EIFSD  forward  to  pre¬ 
dict  the  following:  (1)  crack  size  at  any  time,  t;  (2)  the 
probability  of  exceeding  a  given  crack  size,  x.,  at  any  time, 
t;  and  (3)  the  cumulative  distribution  of  TTCI  at  any  refer¬ 
ence  crack  size,  aQ.  The  SCGMC  depends  on  tha  conditions  to 
be  reflected  in  the  durability  analysis  (e.g.,  loading  spec¬ 
trum,  stress  level,  %  bolt  load  transfer,  etc.)  and  the  basis 
for  the  EIFS  or  EIFSD  to  be  used. 

A  crack  initiation  process  takes  place  in  the  small 
crack  size  range.  Eventually,  as  a  dominant  fatigue  crack 
emerges,  a  crack  propagation  process  takes  over.  The  actual 
crack  size  where  tha  process  changes  from  initiation  to  pro¬ 
pagation  is  debatable.  Suppose  we  assume  that  crack  sizes 
below  0.05"  are  governed  by  tha  initiation  process  and  larger 
cracks  are  governed  by  the  crack  propagation  process.  In 
this  context,  tha  empirical  crack  size-time  relationship  for 
the  small  crack  size  range,  referred  to  as  the  "service  crack 
growth  master  curve"  (SCGMC) ,  then  represents  the  initiation 
process.  This  is  an  "engineering  approach"  to  a  complex 
problem.  Yet,  such  an  approach  is  justified  for  durability 
analysis  until  such  time  that  suitable  mechanistic-based 
models  are  developed  for  practical  applications. 
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In  most  durability  design  situations,  a  suitabls  LEFM 
analytical  crack  growth  program  [e.g., 38,  39]  is  usad  to  de- 
vslop  ths  scgmc  for  ths  dssirsd  analysis  conditions  because 
applicabla  fractographic  raaults  will  not  liksly  bs  avail- 
abls.  Ths  following  gsnsral  procsdurs  for  dsvsloping  a  scgmc 
is  rscommsndsd  for  durability  analysis  applications  in  the 
small  crack  sizs  region. 

1.  Define  the  basis  for  the  EIFSD  to  be  used  in  the 
durability  analysis.  For  example:  (1)  what  fractographic 
crack  size  range,  AL-AU,  was  used?;  (2)  basis  for  EIFS  master 
curve  (e.g.,  given  empirical  crack  growth  law,  LEFM  program 
applied  without  crack  size  restrictions  or  fracture  mechan¬ 
ics-based  program  that  included  the  "short  crack  effect") ; 
and  (3)  what  method  was  used  to  define  the  EIFS  master  curve, 
including  criterion  for  goodness-of-fit  and  crack  shape? 

NOTE:  An  empirical  EIFS-TTCI  relationship  (e.g., 

Eq.  4)  is  recommended  for  general  applica¬ 
tions  to  assure  consistent  EIFSDs  will  be 
obtained  by  different  aerospace  contrac¬ 
tors  for  the  same  fractographic  data  base. 

2.  Use  a  suitable  analytical  crack  growth  program  to 
"curve  fit"  or  tune  the  EIFS  master  curve  or  curves  (where 
the  fractographic  data  pooling  procedure  is  used)  in  the 
fractographic  crack  size  range,  AL-AU,  used  to  define  the 
EIFSD.  In  principle,  it  doesn't  matter  if  the  analytical 
crack  growth  program  is  limited  to  LEFM  applications  or  not 
because  crack  sizes  will  be  limited  to  sizes  where  LEFM  prin¬ 
ciples  apply  anyway.  The  "curve  fit"  to  the  EIFS  master 
curve  in  the  selected  AL-AU  range  is  accomplished  as  follows: 
(1)  plot  the  EIFS  master  curve  to  cover  the  applicable  AL-AU 
range;  (2)  select  a  crack  growth  model  (e.g.,  Walker-AK,  For¬ 
man,  etc.);  (3)  select  da/dN  versus  AK  data  and  calibrate  the 
crack  growth  model  parameters  for  given  material;  (4)  select 


a  load-retardation  nodal  (a.g.,  nodlflad  Willenborg,  Wheeler, 
ate.) ;  and  (5)  by  trial  and  error,  datarmina  tha  ramaining 
nodal  paranatara  raquirad  to  obtain  a  raaaonabla  curva  fit 
using  tha  analytical  crack  growth  progran.  Tha  goodness-of- 
fit  is  dataminad  subjactivaly.  Tha  abova  procadura  is  con- 
captually  dascribad  in  Fig.  12(a),  Tha  "curva  fit"  is  accom- 
plishad  using  tha  applicabla  conditions  raflactad  in  tha  EIFS 
nastar  curva  (i.a.,  load  spactrun,  strass  level,  %  bolt  load 
transfar,  hola  typa/dianatar) . 

3.  Tha  tunad  analytical  crack  growth  progran  is  then 
used  to  predict  tha  crack  growth  over  tha  applicabla  AL-AU 
range  using  tha  applicabla  durability  analysis  conditions. 
For  axanpla,  a  specific  spactrun,  strass  laval  and  %  bolt 
load  transfar  and  assunad  crack  shapa/geonetry  are  used  to 
predict  tha  crack  size  at  a  given  tine.  This  step  is  concep¬ 
tually  illustrated  in  Fig.  12b  for  three  different  stress 
levels  (  0^  >0*2  >0^).  The  crack  size-time  predictions  in 
the  AL-AU  range  are  indicated  in  Fig.  12.  The  analytical 
crack  growth  progran  is  used  to  naka  crack  size-time  predic¬ 
tions  for  crack  size  where  LEFM  principles  apply.  Procedures 
and  assumptions  for  tha  analytical  crack  growth  analysis  are 
assumed  to  be  tha  same  as  those  used  for  a  typical  damage 
tolerance  analysis. 

4.  Crack  size-time  predictions  (i.a.,  a(t)  versus  t) 
in  the  AL-AU  range  for  a  given  strass  region,  depicted  in 
Fig.  12(b),  can  be  used  to  determine  tha  SCGMC  as  follows; 
An  expression  for  tha  SCGMC  can  be  obtained  from  tha  "EIFS 
master  curva,"  given  in  Eq.  4,  by  substituting  y1 ^(T)  “  x  - 
EIFS  and  T  ■  t  to  obtain  Eq.  19.  Y^fT)  is  an  EIFS  in  tha 
EIFSD  corresponding  to  a  crack  size  xx  at  time  7"  in  tha  ith 
stress  region. 

¥xi(T)  -  exp  ( -QjT )  (19) 
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(Analytical  Crack 

Growth  “Curve  Fit"  CRACK  I 


and  Crack  Growth  Model. 


In  Eq.  19,  the  crack  growth  rata  parameter  Q^an  empirical 
constant,  can  be  determined  using  the  crack  size-time 
predictions  for  the  AL-AU  range  as  shown  in  Fig.  12(b). 
Methods  for  estimating  "Qj|  are  given  in  Appendix  C.  As  shown 
in  Fig.  12(b),  LEFM  principles  are  used  only  for  the  crack 
size  range  where  such  principles  apply. 

2.4.2  Probability  of  Crack  Exceedance  and  Correlative 

Distribution  of  TTCI 

Once  the  SCGMC  has  been  determined,  consistent  with  the 
basis  for  the  EIFSD  used,  the  EIFSD  is  grown  forward  using 
the  applicable  SCGMC  (based  on  DCGA)  to  predict  the  probabil¬ 
ity  of  crack  exceedance,  p  (i,T),  at  a  given  time  and  the 
cumulative  distribution  of  TTCI,  FT(t),  for  a  given  refer¬ 
ence  crack  size.  Such  predictions  are  conceptually  described 
in  Fig.  13.  The  SCGMC  defines  the  initial  flaw  size,  y^CT)  , 
which  grows  to  a  crack  size  x1  at  time  T  under  the  desired 
design  conditions.  Given  yli(  T  )  and  the  EIFSD  function 
(e.g.,  Weibull  compatible,  lognormal  compatible,  two-paramet¬ 
er  Weibull,  lognormal,  etc.)  predictions  for  p(i,T)  and 
FT(t)  can  then  be  made  as  described  in  Fig.  13.  The  position 
of  y^CT)  in  the  EIFSD  determines  the  resulting  p(i,7")  and/ 
or  FT(t) . 

The  concept  for  determining  p(i,T)  and  FT(t)  is  very 
general  and  can  be  applied  to  different  EIFSD  functions  and 
SCGMC.  For  the  Weibull  compatible  EIFSD  function,  p(i,  T*  ) 
and  FT(t)  can  be  computed  using  Eqs.  6  and  7,  respectively. 
In  Eq.  6,  x1  is  a  variable  and  T  is  a  fixed  value  for  p(i,T) 
predictions.  In  a  similar  manner,  in  Eq.  7,  x1  is  a  fixed 
value  (i.e.,  reference  crack  size  for  TTCIs)  and  t  is  a  var¬ 
iable  for  FT(t)  predictions. 
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P(i,T)  - 


•robability  of  Crack  Exceedance,  p(i;T)  (*>)  Cumulative  Distribution 


2.5  DURABILITY  ANALYSIS  SOFTWARE 


Computer  software  has  baan  developed  for  an  IBM  PC  for 
implementing  tha  advancad  durability  analysis  methodology 
described  in  this  section.  This  includes  software  for:  (l) 
storing  fractographic  data  on  floppy  disks,  (2)  screening  and 
plotting  the  fractographic  data  for  a  selected  data  set,  (3) 
optimizing  the  EIFSD  parameters  for  one  or  more  fractographic 
data  sets,  (4)  plotting  and  checking  the  goodness-of-fit  of 
p(i,7l»  ra(t)*x*  and  predictions  for  a  given  EIFSD. 
Software  is  also  available  for  implementing  the  durability 
analysis  extension  described  in  Section  III. 

The  durability  analysis  computer  software,  programmed  in 
"GWBASIC",  can  be  implemented  on  an  IBM  compatible  PC.  A 
"user's  guide"  is  available  for  implementing  the  computer 
software  in  Volume  V  [€].  This  software  was  used  extensively 
in  Volume  II [3]  to  evaluate  and  refine  the  advanced  durabil¬ 
ity  analysis  methods  developed  in  this  Volume  (I) . 
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SECTION  III 


DERIVATION  OF  CRACK  EXCEEDANCE  PROBABILITY,  p 
FOR  DURABILITY  ANALYSIS  EXTENSION 

The  main  obj active  of  the  durability  analysis  extension 
is  to  determine  the  probability  of  crack  exceedance  in  the 
large  crack  size  region  to  cover  the  functional  impairment 
due  to  fuel  leakage,  ligament  breakage,  etc.  The  general 
procedure  is  to  compute  the  distribution  of  crack  size  a(r) 
at  any  service  time  r  from  the  distribution  of  EIFS  through 
the  transformation  of  the  service  crack  growth  damage 
accumulation  equation.  Hence,  the  main  issue  is  to  establish 
the  accurate  EIFS  distribution  and  a  service  crack  growth 
behavior.  The  determinations  of  EIFS  distribution  have  been 
described  in  detail  previously.  In  this  section  we  will 
emphasize  the  accurate  description  of  the  service  crack 
growth  behavior  and  the  crack  exceedance  prediction  in  the 
large  crack  size  region. 

The  durability  analysis  methodology  is  simplified  tre¬ 
mendously  if  the  service  crack  growth  damage  accumulation  can 
be  modeled  or  approximated  analytically.  Depending  on  the 
treatment  of  the  crack  growth  damage  in  service,  the  follow¬ 
ing  three  approaches  are  investigated  (1)  deterministic  crack 
growth  approach  (DCGA) ,  (2)  stochastic  crack  growth  approach 
(SCGA) ,  and  (3)  two-segment  deterministic-stochastic  crack 
growth  approach.  Each  approach  will  be  described  in  the 
following  for  predicting  the  probability  of  crack  exceedance 
in  the  large  crack  size  region. 

3.1  SUMMARY  OF  EIFS  DISTRIBUTIONS 

Three  possible  distribution  functions  for  the  equivalent 
initial  flaw  size  (EIFS)  are:  (1)  Weibull  compatible,  (2) 
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lognormal  and  (3)  two-parameter  Weibull,  raapactivaly,  as 
givan  in  tha  following. 

Waibull  Compatibly  EIFSD  Function 

Tha  Waibull  compatible  distribution  function,  proposed 
by  Yang  and  Manning  [13,18],  is  givan  in  Eq.  20 f 

in(x,  ,/x) 

pa(0)<*>  *  — j - ]*}  \  0  1  *  <.  *u  (20) 

*  1  I  *  >  *u 

where  Fa^Qj (x)  -  cumulative  distribution  of  EIFS,  xu  is  the 
upper  bound  of  tha  EIFS  which  is  related  to  tha  reference 
crack  size,  aQ,  and  the  corresponding  lower  bound,  €  ,  of  the 
time-to-crack-initiation  (TTCI)  as  given  in  Eq.  21. 

-Qc 

xu  a  a°*  (21) 


Lognormal  EIFSD  Function 

Fa(0)(x)  *  )  »  *  >  0 

,  n  (22) 

.0  ;  x  <  0 

In  Eq.  22,  Fa(o)  (x)  “  cumulative  distribution  of  EIFS,<&(  )  is 
the  standardized  normal  distribution  function,  and  fi  and 
r  are  the  mean  value  and  standard  deviation  of  jfria(O), 
respectively,  where  a(0)  -  EIFS. 


Two-Parameter  Weibull  distribution 


Fa(0) 


(x) 


-( x/ 80)a° 

C 


x  >  0* 


a  0 


X  <  0 


(23) 


In  Eq.  23,  Fa(0) (x)  =  cumulative  distribution  of  EIFS,  a0  and 
/3  are  the  shape  parameter  and  the  scale  parameter  of  EIFS, 
a(0),  respectively. 
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3.2  DETERMINISTIC  CRACK  GROWTH  APPROACH 


In  this  approach  tha  equivalent  initial  flaw  sizaa 
(EIFSa)  ara  determined  by  back-extrapolating  fractographic 
data  in  tha  crack  aiza  ranga  of  time-to-crack-initiation 
(TTCI) ,  auch  aa  AL-AU ,  using  a  deterministic  crack  growth 
rata  aquation.  Such  EIFS  valuaa  ara  rafarrad  to  aa  “deter¬ 
ministic-based  EIFSa."  Tha  statistical  variability  of  tha 
fractographic  data  involvea  tha  inharant  statistical  dispar- 
sior.  of  tha  ganaric  EIFS  valuaa  aa  wall  aa  that  of  tha  crack 
growth  rata.  Hence,  tha  deterministic-based  EIFS  valuaa ,  in 
fact,  also  include  tha  statistical  variability  of  tha  crack 
growth  rata.  Consequently,  when  tha  datarministic-basad  EIFS 
distribution  is  grown  forward  under  tha  expected  service  load 
spectra,  a  deterministic  crack  growth  rate  model  must  be  used 
up  to  the  crack  size  range  of  TTCI  to  be  consistent  with  the 
fractographic  test  results.  Hence,  the  deterministic  crack 
growth  approach  for  the  durability  analysis  described  in  this 
section  will  be  based  on  the  deterministic-based  EIFS 
distribution.  The  durability  analysis  approach  is  schemati¬ 
cally  shown  in  Fig.  14. 

3.2.1  One-Segment  Deterministic  Crack  Growth  Approach 

In  this  approach,  the  entire  service  crack  growth  master 
curve  (SCGMC)  is  fitted  by  either  one  of  the  following  simple 
but  versatile  crack  growth  rata  models  recommended  by  Yang  and 
Manning  [13,18]; 

da/dt  *  Qt(a)bl  (24) 

or 

da/dt  »  Qta  (25) 
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DISTBiBUTIOII  OF  TTCI 


While  Eqa.  24  and  25  appear  to  ba  simple,  axtanaiva  tast 
raaulta  [16,  58],  including  3  inch  wida  7475-T7351  aluminum 

faatanar  hola  apaoimana  undar  bombar  or  fighter  load  apaotra 
[4],  indioata  that  thay  ara  quita  raaaonabla  [13,18,2,19,21, 
23-30,42] .  Moraovar,  Eq.  24  has  baan  ahown  to  ba  a  raaaonabla 
modal  for  crack  growth  damaga  accumulation  in  critical  loca- 
tiona,  auch  aa  bolt  holaa,  cooling  air  holaa,  ate.,  of  gaa 
turbina  angina  diaka  mada  of  nickal-baaad  auparalloya[40,  41]. 

intagrating  Eq.  25  from  t  -  0  to  any  aarvica  time  t-  r, 
ona  obtaina 

a(0)  «  a(r)  e‘Q*T  (26) 

Tha  distribution  function  of  tha  crack  size  a (r)  can  be 
derived  from  that  of  the  EIFS,  a(0),  given  by  Eqa.  20  -  23 
through  the  transformation  of  Eq.  26  aa  follows: 

F,(t)(»)  -  PC.(t)  <  .  p[,(0)  < 

The  probability  that  the  crack  size,  a(r),  at  any  ser¬ 
vice  time r in  the  ith  stress  region  will  exceed  any  value 
referred  to  as  the  probability  of  crack  exceedance,  is  ob¬ 
tained  as 


p(l,t)  *  PU(t)>xJ  a  1  -  ffa(T)(xJi)  ^28) 

Substituting  Eq.  27  into  Eq.  28  yields 

p(l,r)  a  1  -  pa(0)£y<*i!  t)]  (29) 

where  (x)  is  the  distribution  function  of  EIFS  given  by 

Eq.  20,  22  or  23  and 

y(xt;  t)  a  (3°) 
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Whan  the  more  general  crack  growth  rata  aquation  givan 
by  Eq.  24  can  ba  used  to  approximate  tha  SCOMC  with  euffici- 
ant  accuracy,  tha  intagration  raault  can  ba  axpraaaad  aa  fol- 
lowa: 

a(0) 

«(t)  a . . - . . . -  . .  (31) 

(1  -  [«(0)]°‘  o,0,t|1/e‘ 


in  which 


Than,  tha  crack  axcaadanca  probability,  p(i,r),  ia  givan 
by  Eg.  29,  in  which  y(x1;r)  becomes 

y(*iJT)  *  [Ut)*0*  ♦  (33) 

Tha  one-segment  datarminiatic  crack  growth  approach  dea 
cribad  abova  ia  achamatically  ahown  in  Fig.  14(b). 

3.2.2  Two-Segment  Datarminiatic  Crack  Growth  Approach 

In  soma  cases,  the  service  crack  growth  master  curve 
(SCGMC)  may  not  ba  fitted  by  a  single  aquation,  e.g.,  Eq.  24 
or  25,  with  sufficient  accuracy.  Then,  the  SCGMC  may  be  sep¬ 
arated  into  two  regions;  one  with  the  crack  size  smaller  than 
tha  reference  crack  size,  aQ,  at  crack  initiation,  and  the 
other  with  the  crack  size  larger  than  aQ. 

In  tha  crack  aiza  region  smaller  than  aQ,  the  following 
two  equations,  identical  to  Eqs.  24  and  25,  can  be  used 

da/dt  *  Qt(a)bi  for  a  <  a0  (34) 


or 


da/dt  »  Qxa  for  a  <  a0 


(35) 
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In  the  second  ragion  where  the  crack  aiza  la  largar  than  aQ, 
Gq.  36  can  ba  uaad  for  approximation 

da/dt  ■  QaU)*2  for  a  >  «o  (36) 

in  which  b2  may  ba  equal  to  unity.  Such  an  approach  ia  sche- 
matically  ahown  in  Fig.  14(c).  Tha  crack  growth  rata  para- 
matara  (Q^/b^)  in  tha  email  crack  aiza  ragion  ara,  in 
ganaral,  diffarant  from  (Q2'b2^  in  Eq*  36 • 

Lat  tha  time-to-crack-initiation  for  tha  rafaranca  crack 
aiza  aQ  ba  danotad  by  T,  i.e.,  a(T)  *  aQ.  Whan  tha  crack 
aiza  of  intereat,  ,  ia  amallar  than  aQ,  tha  crack  ex- 
caadance  probability,  p(i,r),  ia  identical  to  that  given  by 
Eqs.  29  and  30.  It  ia  mentioned  that  tha  service  crack 
growth  master  curve  (SCGMC)  in  this  ragion  should  ba  "tuned" 
or  "curve  fitted"  to  ba  consistent  with  the  EIFS  master  curve 
and  than  best-fitted  by  Eq.  34  or  35. 

When  the  crack  size  of  interest  is  greater  than  aQ, 
the  crack  size  a (T)  at  any  service  time  r>  T  can  be  ob¬ 
tained  as  follows.  Integrating  Eq.  35  from  t  -  0  to  t  -  T 
(at  which  a(T)  -  aQ)  leads  to  the  following  expression 

T  *  Q r*  Zn[a0/a(O)]  (37) 

In  the  region  where  a (r)  >  aQ  (orT>  T)  ,  Eq.  36  is  inte¬ 
grated  with  b2  »1  from  t  -  T  to  t  *  T;  with  the  result 

T  *  t  -  Q-1  ZnCa(t)/a03  l  for  t)  >  a0  (38) 

Equating  Eqs.  37  and  38,  one  obtains  the  relationship  between 
the  crack  sizes  a (r)  and  a(0)  as  follows 
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«(0>  ■ 


•*P<A  -  Qtr) 


(39) 


in  which 


A  •  Cl  -  (Qt/Qailmao 


(40) 


Thu*  the  distribution  function,  F, (x),  of  the  crack 
size  a (  r  )  can  ba  obtained  from  that  of  a(0)  through  the 
transformation  of  Eq.  39,  and  tha  probability  of  crack  ax- 
ceadanca,  p(i,*),  can  ba  obtained  in  a  similar  manner  as  Eq. 
29,  i.a. , 

p(i,T)  a  1  -  Fa<0)£y<*li*)] 

in  which 

y(xl;t)  »  (*4)®i^2  exp(A  -  Qtr) 


(41) 


(42) 


where  A  is  given  by  Eq.  40. 


The  solution  obtained  above  holds  for  b1  *  b2  ■  1.  For 
the  case  where  bL  -  1  and  b2  +  1,  the  solution  is  given  by 
Eq.  41  except  that 

y(xl?T)  a  exp  {a  -  Qtr  -  Q^c^)-4^)"®2}  \ 

A  »  Cn«Q  *  Q;(c2Q2)*l(ao)'C2  >  (43) 

c2  *  b2  -  1  / 

Furthermore,  for  the  case  where  b1  t  1  and  b2  +  1,  the  crack 
exceedance  probability  is  also  given  by  Eq.  41  in  which 


y(xtjr)  *  (Xi)C2/ci/rilli  ♦  (x A)®2  c1Q1(A  ♦ 

/  1^2 

A  »  (ctQl)-l(a0)'Ci  -  (c2Q2)-4(a0)'C2 
C  l  *  b  t  "  1  ,  c2  a  b2  -  1 


(44) 
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3.2.3  General  Deterministic  Crack  Growth  Approach 


Tha  aarvica  crack  growth  master  curve  (SCGMC)  can  ba 
obtained  numerically  uaing  an  analytical  crack  growth 
computer  program,  [e.g.,  38,39],  Note  that  tha  SCGMC  variaa 
from  one  atraaa  region  to  another.  Whan  tha  SCGMC  cannot  ba 
approximated  by  two  or  more  segments  of  analytical  models 
(aquations)  with  sufficient  accuracy,  tha  numerically  defined 
SCGMC  can  be  used  directly  for  the  durability  analysis.  The 
solution  for  the  crack  exceedance  probability,  p(i,T)»  can  be 
obtained  numerically  in  the  following  manner  [13,19]. 

The  SCGMC  can  be  represented  symbolically  by  the  follow¬ 
ing  form 


a(tt)  a  yCa(t2)  ;t2-tj  (45) 

in  which  the  crack  size,  a(t1) ,  at  t1  is  expressed  as  a  func¬ 
tion  of  the  crack  size,  a(t2),  at  t2  with  t2  >  ^  and  the 

time  difference  tj-t^  The  function  y  is  a  monotonically  in¬ 
creasing  function  of  a(t2)  and  tj-t^ 

Strictly  speaking,  this  relation  may  not  be  valid  for 
the  general  case  within  one  flight  of  a  particular  mission 
because  of  the  load  sequence  effect  and  other  contributors. 
Thus,  within  one  flight  of  a  particular  mission,  the  rela¬ 
tionship  between  a(t1)  and  a(t2)  depends  on  both  Tx  and  t2* 
Since,  however,  the  design  loading  spectra  in  one  lifetime 
consists  of  many  repeated  missions  and  flights,  it  is  reas¬ 
onable  to  assume  that  tha  relationship  between  a(tx)  and 
a(t2)  depends  only  on  the  time  difference  t2-t1  of  the  ser¬ 
vice  time.  Such  an  expedient  approximation  appears  to  be 
acceptable. 

For  the  special  case  in  which  -  0  and  t2  *  t,  Eq.  45 
becomes 
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t(0)  ■  yC«( T)| t] 


(46) 


The  symbolic  equation,  such  as  Eq.  45  or  46,  will  serve, 
mathematically,  as  a  transfer  function  for  transforming  the 
statistical  distribution  of  the  crack  population  from  one 
service  time  (such  as  t  ■  0)  to  another. 


The  distribution  function  of  the  crack  size  a(T)  at  any 
service  time  T  can  be  derived  from  that  of  a(0)  through  the 
transformation  of  Eq.  46  as  follows 

Fa(  T)  *  Pta(  t)  £  *]  »  P[a(0)  £  y(x;  t)]  =»  pa(Q)Cy(x;r)] 
or  (47) 

F.(t)U)  *  F.(0)t)'<xiT)1 

in  which  the  property  that  y(x;T)  is  a  monotonically  increas¬ 
ing  function  has  been  used.  The  probability  of  crack  exceed¬ 
ance,  p(i,T),  is  obtained  as 


p(i,T)  «  PCa(r)  >  xt]  »  1  •  Fa(T)(*i) 
*  1  *  F*(0)Cy(*iJT)] 


(48) 


where  the  distribution  of  EIFS,  Fa^0)  (x) ,  is  given  by  Eq.  20, 
22,  or  23.  The  argument  y(x1?T)  can  be  determined  graphic¬ 
ally  (or  numerically  with  interpolation  procedures)  from  the 
SCGMC  as  shown  schematically  in  Fig.  15.  The  numerical  pro¬ 
cedures  for  the  determination  of  y(x1;7")  are  described  as 
follows;  (1)  Given  the  crack  size  of  interest,  x^,  determine 
the  corresponding  time  t^  from  the  SCGMC  by  interpolation, 
and  (2)  Determine  the  crack  size  y(x1;T)  corresponding  to  the 
time  -  T  from  the  SCGMC  again  using  the  interpolation  pro¬ 
cedures  (Fig.  15) . 


The  computation  of  the  crack  exceedance  probability, 
p(i,T),  described  above  is  straight-forward  and  simple. 
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Figure  15.  General  Deterministic  Creek  Growth  Approach 
for  Service  Crack  Growth  Master  durve. 
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3.3  STOCHASTIC  CRACK  GROWTH  APPROACH 


Analysis  results  of  fractographic  data  indicate  that  the 
crack  growth  damage  accumulation  in  the  small  crack  size  re¬ 
gion  involves  considerable  statistical  variability  [e.g.,  27 
42-44].  Theoretically,  the  statistical  variability  of  the 
fractographic  data  in  the  crack  size  region  close  to  TTCI  or 
the  TTCI  fractographic  data  is  contributed  by  the  statistical 
dispersions  of  both  the  equivalent  initial  flaw  size  and  the 
crack  growth  damage  accumulation  up  to  aQ.  In  the  determin¬ 
istic  crack  growth  approach,  the  distribution  of  the  EIFS  is 
determined  by  back-extrapolating  the  fractographic  data  using 
a  deterministic  crack  growth  model.  The  statistical 
dispersion  of  the  EIFS  obtained  also  includes  the  statistical 
dispersion  of  the  crack  growth  rate  in  the  small  crack  size 
region.  Such  an  approach  is  quite  reasonable  as  long  as  the 
deterministic  crack  growth  rate  model  is  used  to  grow  the 
EIFS  values  up  to  the  reference  crack  size  aQ.  Such  an 
expedient  approach  simplifies  the  durability  analysis 
methodology  tremendously,  and  the  simplification  makes  the 
proposed  durability  analysis  very  attractive. 

Another  approach  [28,29]  is  to  treat  the  crack  growth 
rate  as  a  stochastic  process,  and  the  EIFS  values  are  obtain¬ 
ed  by  back-extrapolating  the  fractographic  data  stochastic¬ 
ally  as  shown  in  Fig.  16(a).  This  implies  that  the 
statistical  variability  of  the  crack  growth  rate  is  filtered 
out  in  the  back-extrapolation  process.  The  EIFS  values  thus 
obtained  are  referred  to  as  "stochastic-based  EIFSs."  Then, 
under  the  expected  service  load  spectra,  the  forward  crack 
propagation  from  EIFS  must  be  stochastic.  Such  an  approach 
will  be  addressed  in  this  section  and  the  stochastic-based 
EIFS  will  be  used  for  the  computation  of  the  crack  exceedance 
probability. 
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Figure  16.  Stochastic  Crack  Growth  Approach  for 
Durability  Analysis  Extension. 


Th«  stochastic  crack  growth  approach  depends  on  a 
simple,  yet  reasonable,  stochastic  crack  growth  rate  model. 
Recently,  a  simple  and  plausible  stochastic  crack  growth 
model  proposed  by  Vang,  et  al  [27,40-42,45-51]  can  be  used 
conveniently  for  durability  analysis  as  follows 

da/dt  a  XQ^a)1^  (49) 

in  which  X  is  a  lognormal  random  variable  with  a  median  value 
equal  to  unity,  i.e.,  X  takes  values  around  unity.  Thus,  the 
deterministic  crack  growth  rate  model  given  by  Eq.  24  repre¬ 
sents  the  median  crack  growth  behavior,  whereas  the  statisti¬ 
cal  variability  of  the  crack  growth  rate  is  taken  care  of  by 
the  lognormal  random  variable  X.  Such  a  stochastic  crack 
growth  model  is  referred  to  as  the  lognormal  random  variable 
model  [27,50,51].  The  stochastic  crack  growth  approach  for 
the  durability  analysis  is  schematically  shown  in  Fig.  16. 

While  many  different  stochastic  crack  growth  models  have 
been  suggested  in  the  literature  [e.g.,  43,44,52-56],  the 
lognormal  random  variable  model,  Eq.  49,  is  the  simplest 
stochastic  model  most  suitable  for  practical  applications, 
such  as  the  durability  analysis.  It  has  been  demonstrated 
[e.g.,  27,50,51]  that  Eq.  49  can  be  used  satisfactorily  for 

the  fatigue  crack  propagation  in  7475-T7351  fastener  holes 
under  bomber  and  fighter  loading  spectra  in  both  the  small 
crack  size  and  large  crack  size  regions. 

Since  X  is  a  lognormal  random  variable  with  a  median 
value  equal  to  unity,  the  random  variable 

Z  *  log  X  (50) 

is  normal  with  zero  mean  and  standard  deviation  <r .  The 

z 

probability  density  function  of  X  is  given  by 


? 


u  >  0 


f.(u)  .  -i°S  *-  .J-  1  (  i*SLH  )* 

ST*  x  9t  (2  at 
*  0 


u  <  0 


(51) 


3.3.1  One-Segment  Stochastic  Crack  Growth  Approach 


Suppoaa  tha  atochaatic  aarvica  crack  growth  master  curve 
can  bo  represented  by  Eq.  49.  Then,  integrating  Eq.  49  from 
t  »  0  to  t  ■  T,  for  ■  l,  one  obtains 


a(t)  -  «<0)eXQiT 


(52) 


or 


a(0)  a  a(t)e‘XQ*T  (53) 

The  crack  size,  a(r),  at  any  aarvica  time  Tie  a  function 
of  two  random  variables,  a(0)  and  X.  To  derive  the  distribu¬ 
tion  function  of  a(T),  a  particular  value  u  is  assigned  to  X, 
i.e.,  X  -  u,  and  then  the  theorem  of  total  probability  is 
used.  The  conditional  distribution  function  of  a (T)  given  X 
*  u  is  derived  as  follows 

*  RC«<  T)<Xt|X«u]  ( 


»  P(a(0)<y(xt;  tjxau)]  *  Fa(0)C>'<*L}T|x*u)] 
in  which  it  follows  from  Eq.  53  that 

ytx^TjXau)  s  xte'u<3lT  (55) 

The  conditional  crack  exceedance  probability,  denoted  by 
p(i,  r|x-u) ,  given  X  ■  u,  is  obtained  as 


p(i,tjXsu)  =  1  -  fra(  T)(*;|Xau> 

3  1  *  Fa(0)Cy(*l|X3u)] 


(56) 


in  which  Eq.  54  has  been  used. 
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The  unconditional,  crack  exceedance  probability  is  ob- 
tainad  from  ths  conditional  one  using  ths  thaoram  of  total 
probability  as 

m 

■  /  p(i,  r|x«u)fx(u)du 

*  1  •  /  Fa^Q)Cy(xAi  t|x*u)]fx(u)du  *  ^ 

in  which  y^/Tjx-u)  and  fx(u)  ara  givan  by  Eqs.  55  and  51, 
raspactivaly. 


Whan  b ^  +  1 ,  Eg.  49  can  ba  intagratad  analytically  and 
the  crack  axcaadanca  probability,  p(i,r),  is  idantical  to 
that  prasantad  in  Eq.  57,  axcapt  that  y(xx;  t|x  -  u)  is  givan 
by 


y(*i » tjx*u) 


_ Kj _ 

Cl  -  <xt)C*  ClQttu]1/Ci 


(58) 


in  which  C1  «  bx  -1. 


The  solution  darivad  abova,  Eq.  is  exprassad  in  tarrns 
of  the  distribution  function  of  EIFS,  Fa(0)^x^  that  is  given 
by  Eq.  20  in  the  case  of  the  Weibull  compatible  distribution, 
by  Eq.  22  in  the  case  of  lognormal  distribution,  and  by  Eq. 
23  in  the  case  of  the  two-parameter  Weibull  distribution. 
For  instance,  when  the  Weibull  compatible  distribution  is 
used,  the  crack  exceedance  probability  given  by  Eq.  57 
becomes 


p(l,t)  *  1 


ir\ 


xu 

y(*i?  t 


X»g) 


f^(u)du 


(59) 


In  the  numerical  computation  using  Eq.  59,  it  should  be  noted 
that  the  term  in  the  integrand  e*p(- (inCxu/y(*As  t|x»u)]/4|°)  *  1  If 
y ( x A { t | X*u )  >  xu. 
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The  crack  axcaadanca  probability  praaantad  in  Eq.  57 
should  ba  avaluatad  numerically,  ainca  tha  analytical  inte¬ 
gration,  in  general,  is  not  faaaibla.  Howavar ,  tha  numarical 
intagration  is  simple  and  straight-forward.  Tha  one-segment 
stochastic  crack  growth  approach  is  schamatically  shown  in 
Fig.  16(b). 


3.3.2  Two-Segment  Stochastic  Crack  Growth  Approach 

Whan  tha  service  crack  growth  master  curva  cannot  ba  ap¬ 
proximated  by  a  single  crack  growth  rata  aquation,  Eq.  49, 
with  sufficient  accuracy,  a  two-segment  crack  growth  approxi¬ 
mation  analogy  to  that  described  in  Section  3.2.2  should  ba 
used. 

In  the  crack  size  region  smaller  than  the  reference 
crack  size  aQ,  the  crack  growth  rate  equation  is  approximated 
by 

da/dt  *  XQt(a)bl  i  for  i  <  i0  or  t  <  T 

(60) 

In  the  crack  size  region  larger  than  aQ,  the  following  ap¬ 
proximation  is  used 

da/dt  a  XQ2(a)b2  .  for  a  >  a0  or  t  >  T  (61) 

Integrating  Eq.  60  from  t  ■  0  to  t  >  T  (at  which  a(T)  - 
aQ)  ,  with  *  1,  leads  to  the  following  expression 


a(T)  ■  a(0)e*^lT  a  ag 


from  which 


T  «  _L  m[  -ffi-  1 

XQi  a(0) 


(62) 


(63) 


Equation  61  can  be  integrated  from  T  (i.e.,  aQ)  to  any 
service  time  T >  T, 
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«(  t)  ■  «0« 


(64) 


KQ2(  t-T) 


Substituting  Eq.  63  into  Eq.  64  lsads  to  ths  relation  bstwaan 
a (T)  and  a(0)  as  follows 

«(0)  -  CaCrJl^o)1^  #-XQ2ty  (65) 


in  which 


Y  *  Oi/Qj  (66) 

whsn  ths  crack  sizs  of  intarast,  x. ,  is  smaller  than  tha 
rafaranca  crack  siza  aQ ,  i.e.,  x1  <  aQ ,  tha  probability  of 
crack  axcaadanca  is  idantical  to  Eqs.  57  and  55.  Whan  the 
crack  siza  of  interest,  x1#  is  larger  than  aQ,  however,  Eq. 
65  should  be  used  to  derive  tha  crack  exceedance  probability; 
with  the  result 

m 

p(ifT)-1-^/  Pa(0)Cy(*i;t)|x»u)fx(u)du  ;  for  x,  >  a0  (67) 
in  which 


y(x1jT|x*u)  ■  xJ(a0)1*Y  e’u®*TY  (68) 

Again,  for  the  Weibull  compatible  distribution  given  by 
Eq.  20,  the  term  Fa(0) Cy(x1;  Tjx-u) ]  in  the  integrand  of  Eq. 
67  is  equal  to  unity  if  y(xx;  r|x»u)  is  greater  than  xu,  i.e., 

Fa(0)  Cy(xl;  rlX“u)  3  “  1  for  y(x1'*'T|x“u)  >  xu- 

When  bx  and  b2  in  Eqs.  60  and  61  are  not  equal  to  1.0, 
i.e.,  b1  f  and  b2  +  1,  the  crack  exceedance  probability  de¬ 
rived  in  Eq.  67  is  still  valid  except  that  y(x1;r|x«u)  is 
different;  however,  it  can  be  derived  easily.  The  two-seg¬ 
ment  stochastic  crack  growth  approach  is  schematically  shown 
in  Fig.  16(c) . 
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3.3.3  Genoral  stochastic  Crack  Growth  Approach 

Whan  tha  aarvica  crack  growth  master  curva  cannot  ba 
approximated  by  two  or  mora  sagmants  of  stochastic  crack 
growth  rata  aquations  with  sufficient  accuracy,  numarically 
dafinad  general  stochastic  service  crack  growth  master  curves 
should  be  used.  The  derivations  and  numerical  procedures  for 
computing  tha  crack  exceedance  probability  are  similar  to 
those  described  in  Section  3.2.3.  The  only  difference  is 
that  the  crack  exceedance  probability  obtained  in  Section 
3.2.3  corresponds  to  a  particular  percentile  with  certain 
probability,  and  a  sufficiant  number  of  repetitive 
computations  should  ba  made  to  cover  a  wide  range  of 
percentiles. 


3.4  TWO-SEGMENT  DETERMINISTIC-STOCHASTIC 
CRACK  GROWTH  APPROACH 

Two  approaches  for  establishing  the  equivalent  initial 
flaw  size  have  been  considered;  namely,  the  deterministic- 
based  EIFS  and  the  stochastic-based  EIFS.  For  the  determin¬ 
istic-based  EIFS,  a  deterministic  crack  growth  rate  model 
should  be  applied  to  grow  the  entire  EIFS  population  up  to 
the  reference  crack  size  aQ,  whereas  a  stochastic  crack 
growth  rate  model  should  be  used  to  grow  the  stochastic-based 
EIFS  population.  It  has  been  shown  [28,29]  that  both  ap¬ 
proaches  are  reasonable  for  durability  analysis.  However, 
the  deterministic-based  EIFS  concept  appears  to  be  more  att¬ 
ractive  for  two  reasons:  (1)  The  mathematical  expressions 
and  computations  are  much  simpler,  and  (2)  within  the  short 
crack  size  range,  e.g.,  crack  size  smaller  than  0.005",  the 
crack  growth  behavior  of  short  cracks  is  not  well  understood 
and  the  applicability  of  linear  elastic  fracture  mechanics  is 
not  certain.  Thus,  the  EIFS  is  simply  a  hypothetical  concept 
of  convenience  for  conducting  various  types  of  analyses,  such 


as  ths  durability  analysis.  Hence,  it  is  not  absolutsly 
nsosssary  to  uss  ths  stochastic-  based  EIFS  concept  for  the 
purpose  of  practical  applications. 


Because  of  simplicity  in  the  durability  analysis  proce¬ 
dures,  the  deterministic-based  EIFS  concept  will  be  consider¬ 
ed  in  this  section,  and  a  deterministic  crack  growth  rate 
model 


da/dt  »  Qt(a)bl 


for  a  <  aQ  or  t  <  T 


(69) 


will  be  used  to  grow  the  deterministic-based  EIFS  population 
up  to  the  reference  crack  size  aQ  to  be  consistent  with  the 
fractographic  results. 


When  the  crack  size  of  interest  is  extended  beyond  the 
reference  crack  size  aQ,  such  as  that  associated  with  the 
functional  impairment,  the  effect  of  the  statistical  variab¬ 
ility  of  the  crack  growth  rate  influences  the  durability 
analysis  results.  In  particular,  such  an  effect  on  the  pre¬ 
diction  of  the  crack  exceedance  probability  must  be  accounted 
for.  As  a  result,  the  stochastic  crack  growth  rate  model 
described  previously  will  be  employed  to  grow  the  cracks  in 
the  crack  size  region  larger  than  the  reference  crack  size 
8q/  i . e . , 

da/dt  *  XQ2(a)b*  ;  for  a  >  a0  (70) 


in  which  X  is  a  lognormal  random  variable  with  a  median  value 
equal  to  unity.  The  approach  described  above  is  referred  to 
as  the  two-segment  deterministic-stochastic  approach  [30]. 

Integrating  Eq.  69  from  t  ■  0  to  t  «  T  for  b1  »  1,  one 
obtains 


T  a  Q*ltn[ao/a(0)3 


(71) 
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in  which  it  is  understood  that  a(T)  *  aQ . 

In  the  region  where  a {*r)  >  a Q  (or  T>  T) ,  Eq.  70  is  inte¬ 
grated  with  b2  ■  1  from  t  ■  T  to  t  »T( or  from  a(t)  *  aQ  to 
a(t)  -  a (T) ) ;  with  the  result 

T  «  t  -  (XQ2)-itnC*(T>/a03  >  a(  t)  >  a0  (72) 

Equating  Eqs.  71  and  72  leads  to  the  following  relation 
between  a(T)  and  a(0) 

a<0>  .  a0e’QiT[a(T)/a0]r/x  for  a<  t)  >  (73) 

in  which 


Y  *  Qi/02  (74) 

When  the  crack  size,  a (T),  at  any  service  time  is 
smaller  than  aQ,  the  relation  between  a(T)  and  a(0)  is  ob¬ 
tained  by  integrating  Eq.  69  from  t  ■  0  to  t  ■  ^ 

a<0)  »  a(t)e"^iT  ;  a(  t)  <  a0  (75) 

The  crack  exceedance  probability,  p(i,T),  can  be  comput¬ 
ed  depending  on  the  crack  size  of  interest  x1  in  the  follow¬ 
ing  manner. 

(1)  When  the  crack  size  of  interest  x1  is  smaller  than 
the  reference  crack  size  aQ,  the  distribution  function 

Fa(  r )  (xl)  of  the  cracJc  size  aCn  *or  xi  <  ao  is  d®^ivad  from 
the  distribution  function  of  a(0)  through  the  transformation 

of  Eq.  75, 

f.(t)u‘)  •  F.<0)C>'<’,‘|,)5  (76) 
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in  which 


ytXiJT)  a  Xte’Q*T  (77) 

and  the  crack  exceedance  probability,  p(i,T),  is  given  by 

p<i,t)  *  1  -  F4(0)Cy<*i,t)3  for  xA£a0  (78) 

(2)  When  the  crack  size  of  interest  x1  is  larger  than 
aQ,  the  conditional  distribution  function  of  a (T)  at  any  ser¬ 
vice  time  T  given  X-u  can  be  derived  from  that  of  a(0) 
through  the  transformation  of  Eq.  73.  Then,  the  uncondition¬ 
al  distribution  function,  F^^fx^),  of  a (T)  can  be  obtained 
by  using  the  theorem  of  total  probability;  with  the  result 

Fa(r)^x^  a  Qf  Fa(0)  IG<xi*  t|x»u)  ]fx(u)du  (79) 


in  which 


CCx^jTjXau)  a  ao«"^tT  (xt/a0)Y^U  (80) 

The  crack  exceedance  probability  for  x1  >  aQ  is  given  by 
p(i,t)  ■  1  -  /  Fa(0)[G<xt}t|X-u)]fx(u)du  (81) 

in  which  G(x1?  7'jx-u)  is  given  by  Eq.  80  and  fx(u)  is  the  pro¬ 
bability  density  function  of  the  lognormal  random  variable  X 
given  by  Eq.  51. 

When  the  Weibull  compatible  distribution,  Eq.  20,  is 
used  for  the  EXFS,  the  condition  that  Fa^j  [G(x1;  t|x*u)  ]  -  1 
for  G(x1;  t|x-u)  >  xu  should  be  accounted  for  in  the  computer 
program  for  computing  the  crack  exceedance  probability  P(i,T)» 
Eq.  81.  The  two-segment  deterministic-stochastic  crack 
growth  approach  described  above  is  schematically  shown  in 
Fig.  17. 
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3.5  DISTRIBUTION  OF  TIME  TO  REACH  ANY  CRACK  SIZE  x2 

The  determination  of  the  crack  exceedance  probability, 
p(i,T),  in  both  the  large  and  email  crack  aize  regions  has 
been  presented.  Another  quantity  of  interest  in  the  durabil¬ 
ity  analysis  is  the  distribution  of  service  time  to-  reach  any 
specific  crack  size  x^.  Let  T(x^)  be  the  time  for  a  crack  to 
reach  any  size  x^  and  FT^x  ^(T')be  the  corresponding  cumulative 
distribution  function,  i\e.,  FT^X  jT)  ■  P[T(x1)<T3-  The 
distribution  function  of  T(x1)  is  tiie  probability  that  the 
crack  will  reach  a  size  x^  before  the  service  time  T.  Such  a 
probability  is  equal  to  the  probability  that  the  crack  size 
a(T)  at  service  time  Twill  exceed  x^,  which  is  simply  the 
probability  of  crack  exceedance.  Hence, 

FT(Xl)(r)  *  P[T(x1)-T1  “  PtaCDSxJ  -  p(i,T)  (82) 

Consequently,  the  distribution  function  of  the  service  time 
to  reach  any  crack  size  xx  is  obtained  by  computing  the  crack 
exceedance  probability,  p(i,T),  at  different  values  of  T. 

3.6  SERVICE  CRACK  GROWTH  MASTER  CURVE  FOR 
DURABILITY  ANALYSIS  EXTENSION 

A  general  procedure  for  determining  a  suitable  service 
crack  growth  master  curve  (SCGMC)  for  growing  the  EIFSD  for¬ 
ward  is  described  and  discussed  in  this  section.  The  SCGMC 
is  determined  using  a  general  crack  growth  computer  program 
[e.g.,  38,39]  and  a  suitable  crack  growth  model  (e.g.,  Eq. 
35) .  It  is  not  essential  that  a  single  SCGMC,  composed  of 
two  segments,  actually  be  generated  to  implement  the  durabil¬ 
ity  analysis  extension.  The  SCGMC  can  be  represented  by 
crack  growth  rate  models  (e.g.,  Eq.  34  and  35),  and  the  ex¬ 
pressions  for  p(i;T)  and/or  FT(t)  contain  the  applicable 
model  parameters.  However,  it  is  important  to  understand  how 
the  SCGMC  is  developed  and  how  it  works  so  that  the  needed 
crack  growth  model  parameters  can  be  determined. 
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In  the  following  the  two-segment  DCGA-DCGA  will  be  em¬ 
phasized  and  applications  for  the  DCGA-SCGA  will  be  discuss¬ 
ed.  The  first  segment  covers  the  small  crack  size  range 
(i.e.,  0<a(t)<AU)  and  the  second  segment  covers  the  large 
crack  size  range,  AU-AU1  as  shown  in  Fig.  18. 

1.  A  suitable  analytical  crack  growth  program  [e.g., 
38,39]  is  "tuned”  or  "curve  fitted"  to  the  basis  for  the 
EIFSD.  This  involves  curve  fitting  the  analytical  crack 
growth  program  to  the  EIFS  master  curves  reflected  in  the 
EIFSD,  in  the  AL-AU  range.  Refer  to  Step  2  in  Section  2.4.1. 

2.  Use  the  tuned  analytical  crack  growth  program  to 
predict  the  crack  growth  (i.e.,  a(t)  versus  t)  from  AL  to  AUr 
for  the  desired  analysis  conditions  (i.e.,  stress  level,  load 
spectrum,  %  bolt  load  transfer,  etc.).  For  example,  such 
crack  growth  predictions  are  shown  in  Fig.  18  for  three  dif¬ 
ferent  stress  levels,  in  which  >0^  In  Fig.  18  the 

analytical  crack  growth  predictions  cover  a  crack  size  range 
where  LEFM  principles  apply  (i.e.,  a(t)>  AL) . 

3.  The  SCGMC  for  segment  1  (i.e.,  0<a(t)<AU)  is  deter¬ 
mined  using  a  suitable  crack  growth  model  (e.g.,  Eg.  35)  to 
fit  the  analytical  crack  growth  predictions  in  the  AL-AU 
range.  in  Eq.  35  can  be  determined  using  the  procedure 
described  in  Appendix  C.  Similar  methods  for  determining  (^ 
and  b1  in  Eq.  34  are  available  [15]. 

4.  In  a  similar  manner,  the  SCGMC  for  segment  2  (i.e., 
AU<a(t)<Au'  )  can  be  determined  using  a  suitable  crack  growth 
model  as  depicted  in  Fig.  18  to  fit  the  analytical  crack 
growth  predictions  in  the  AL-AU  range.  Use  the  same  methods 
described  previously  to  determine  the  model  parameter (s) . 

5.  The  two  crack  growth  master  curve  segments  for  a 
given  stress  level  can  be  physically  combined  into  a  single 
SCGMC  as  illustrated  in  Fig.  19.  At  point  0  segments  1  and 
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A|  9  -  Crack  size  -  time  predictions 


2  have  the  same  a(t),t  values  but  not  necesaarily  the  sane 
slopes.  The  two-segment  SCGMC  is  functionally  the  same  as 
that  described  in  Fig.  15.  For  example,  the  EIFS,  yu(;-)  • 
corresponding  to  crack  size  x^^  at  any  service  time  T  is  shown 
in  Fig.  19. 


Actually,  the  two-segment  SCGMC  for  the  DCGA-DCGA  can 
also  be  used  for  the  DCGA-SCGA  for  the  same  crack  growth  an¬ 
alysis  conditions  (e.g.,  stress  level,  load  spectra,  etc.). 

The  only  difference  is  that  the  parameter  is  required  to 

implement  the  DCGA-SCGA.  <j*_  is  the  standard  deviation  of 

z 

da/dt  with  respect  to  the  plot  of  In  da/dt  «  In  a(t)  +  Qt.  <r 

z 

can  be  determined  from  suitable  fractographic  results,  if 
available.  Recommended  values  of  <rz  and  guidelines  for  both 
straight  bore  and  countersunk  fastener  holes  are  presented 
elsewhere  [3,7]. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 
4.1  CONCLUSIONS 

An  advanced  durability  analysis  (ADA)  methodology  has 
been  developed  for  metallic  aircraft  structures.  This  meth¬ 
odology  is  very  general  and  it  can  be  used  to  evaluate  dur¬ 
ability  design  requirements  for  functional  impairment  due  to: 
(1)  excessive  cracking,  and  (2)  fuel  leakage/ ligament  break¬ 
age.  The  ADA  methodology  accounts  for  the  initial  fatigue 
quality  variation  of  structural  details,  the  crack  growth 
accumulation  for  a  population  of  structural  details  under 
specified  durability  analysis  conditions  (e.g.,  material, 
load  spectra,  stress  level,  %  bolt  load  transfer,  etc.)  and 
structural  properties. 

An  "engineering  approach",  based  on  engineering  prin¬ 
ciples,  rather  than  mechanistic-based  theories  for  micro- 
structural  behavior,  is  used  to  make  the  methodology  useful 
for  durability  design  applications.  Hence,  the  ADA  method¬ 
ology  is  considered  to  be  a  "durability  design  tool."  This 
perspective  is  very  important  for  durability  analysis  appli¬ 
cations. 

The  initial  fatigue  quality  of  a  structural  detail  (e.g. 
fastener  hole,  cutout,  fillet,  lug,  etc.)  is  represented  by 
an  equivalent  initial  flaw  size  distribution  (EIFSD)  func¬ 
tion.  The  EIFDS  is  based  on  the  back-extrapolation  of  frac- 
tographic  results  to  time  zero  using  a  deterministic  crack 
growth  approach  (DCGA)  .  Analytical  method;*  hav*  been  devel¬ 
oped  and  verified  for  determining  a  suitable  EIFSD  for  the 
durability  analysis,  of  clearance-fit  fastener  holes  in  metal¬ 
lic  aircraft  structures. 

The  ADA  methodology  reflects  a  probabilistic  approach 
for  predicting:  (1)  the  probability  of  crack  exceedance, 


71 


p(i,'7')#  At  any  service  time  and/or  (2)  tha  cumulative  distri¬ 
bution  of  TTCI ,  rT(t) ,  at  any  crack  siza.  Datarminiatic  and 
stochastic  crack  growth  mathods  have  baan  davalopad  for  mak¬ 
ing  p(i,T)  and  rT(t)  predictions  for  any  crack  siza  and  ser¬ 
vice  tima.  For  example,  tha  ADA  methodology  applies  to  tha 
small  crack  siza  range  associated  with  excessive  cracking 
(a.g.,  <  0.05")  and  to  large  through-  the-thicknass  cracks 
(e.g.,  0.50"  -  0.75")  associated  with  fuel  leakage/ligament 
breakage . 

The  ADA  methodology  presented  in  this  report  (Vol.  I)  is 
evaluated  and  verified  in  Volume  II  [3).  The  following  con¬ 
clusions  and  understandings  are  based  on  the  work  of  this  re¬ 
port  as  well  as  Volume  II. 

1.  Actual  initial  flaws  in  the  bore  of  manufactured 
fastener  holes  in  metallic  aircraft  structures  usually  con¬ 
sist  of  random  scratches,  burrs,  microscopic  imperfections, 
etc.  Different  forms,  locations  and  combinations  of  flaws 
exist  in  each  fastener  hole.  Such  flaws,  except  for  gross 
mamifacturing  defects,  cannot  be  reliably  detected  and  quan¬ 
tified  by  NDE  for  production  aircraft  structures.  In  real¬ 
ity,  the  actual  initial  flaws  in  fastener  holes  produced  by 
manufacturing  and  assembly  are  not  physical  "cracks"  in  the 
usual  sense  associated  with  the  linear  elastic  fracture  mech¬ 
anics  approach.  Whatever  the  source  of  fatigue  cracking  may 
be,  a  practical  method  for  representing  the  reality  of  the 
as-manufactured  condition  is  needed  for  durability  analysis. 
This  is  taken  care  of  by  the  equivalent  initial  flaw  size 
concept  described  in  the  following. 

2.  An  equivalent  initial  flaw  (EIFS)  is  an  artificial 
crack  size  which  results  in  an  actual  crack  size  at  an  actual 
point  in  time  when  the  initial  flaw  is  grown  forward.  It  is 
determined  by  back-extrapolating  fractographic  results.  It 
has  the  following  characteristics:  (1)  an  EIFS  is  an  arti- 
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ficial  crack  assumed  to  rapraaant  tha  initial  fatigue  quality 
of  a  atructural  datail  in  tha  as -manufactured  condition  what- 
avar  tha  aourca  of  fatiqua  cracking  may  ba,  (2)  it  haa  no  di- 
ract  relationship  to  actual  initial  flawa  in  faatanar  holes 
such  as  scratches,  burrs,  microdefects,  ate.,  and  it  cannot 
ba  verified  by  NDI,  (3)  it  has  a  universal  crack  ahapa  in 
which  tha  crack  size  is  measured  in  tha  direction  of  crack 
propagation,  (4)  EIFSs  are  in  a  fracture  mechanics  format  but 
they  are  not  based  on  linear  elastic  fracture  mechanics 
(LEFM)  principles  and  are  not  subject  to  such  laws  or  limita¬ 
tions  such  as  the  "short  crack  effect,"  (5)  it  depends  on  the 
fractographic  data  used,  the  fractographic  crack  size  range 
used  for  the  back-extrapolation  and  the  crack  growth  rate 
model  used,  (6)  it  must  be  grown  forward  in  a  manner  consis¬ 
tent  with  the  basis  for  the  EIFS,  (7)  EIFSs  are  not  unique  - 
a  different  set  is  obtained  for  each  crack  growth  law  used 
for  the  back-extrapolation,  and  (8)  EIFSs  are  not  necessarily 
comparable  with  minimum  crack  sizes  obtained  from  fracto¬ 
graphic  results  from  tear-down  inspections. 

3.  An  EIFS  is  not  strictly  "generic"  because  it  depends 
on  the  following:  (1)  crack  growth  rate  model  used  to  back- 
extrapolate  fractographic  results,  (2)  conditions  reflected 
in  the  fractographic  results  (e.g.,  material,  type  fastener/ 
hole/fit,  load  spectra,  stress  level,  %  bolt  load  transfer, 
etc.),  (3)  crack  growth  approach  used  (deterministic  or 
stochastic),  (4)  fractographic  crack  size  range  used  (i.e., 
AL-AU) ,  and  (5)  goodness-of-f it  criterion.  However,  the  real 
issue  is  not  whether  the  EIFSs  or  EIFSD  is  generic  or  not. 
The  important  question  is:  "Can  an  EIFSD,  based  on  the  frac¬ 
tographic  results  for  one  or  more  data  sets,  be  used  to  make 
reasonable  durability  analysis  predictions  for  a  different 
set  of  conditions  (e.g.,  similar  material,  same  type  of  load 
spectra  (e.g.,  fighter,  bomber  or  transport),  similar  type 
fastener/hole/fit  but  different  stress  levels  and/or  %  bolt 
load  transfers)?"  The  answer  to  this  question  is  "yes"!  A 
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fractographic  data  pooling  procadura  has  baan  davalopad  and 
verifiad  for  determining  auitabla  EIFSD  for  durability  analy¬ 
sis. 


4.  An  EIFS  and/or  EIFSD  should  ba  usad  in  tha  propar 
context  for  durability  analysis  applications.  EIFSs  roust  ba 
grown  forward  in  tha  saroa  nannar  as  such  EIFSs  wars  defined 
in  the  first  place.  For  example,  an  artificial  initial  flaw 
size  or  EIFS  should  not  be  indiscriminately  grown  forward  to 
any  desired  crack  size  using  a  LEFM-based  crack  growth  pro¬ 
gram  without  first  considering  the  "basis"  for  the  EIFS  or 
EIFSD.  If  this  principle  is  not  strictly  followed,  durabil¬ 
ity  analysis  predictions  will  not  be  consistent  with  the 
EIFSD  and  the  predictions  may  not  be  reasonable. 

5.  Initial  flaws  can  be  selected  from  the  EIFSD  for  a 
given  probability  of  crack  exceedance  (e.g.,  1/10000  or 
0.0001) .  Such  initial  flaws  are  artificial  -  not  actual  phy¬ 
sical  flaws  that  can  be  compared  with  the  minimum  crack  size 
obtained  from  the  fractographic  evaluation  of  tear-down  in¬ 
spection  results. 

6.  The  ADA  methodology  includes  a  fracture  mechanics 
framework.  However,  an  empirical-based  crack  growth  model  is 
used  in  the  small  crack  size  region  rather  than  linear  elas¬ 
tic  fracture  mechanics  (LEFM)  principles.  LEFM  analytical 
crack  growth  methods  are  used  only  for  crack  sizes  and  stress 
intensity  ranges  where  such  principles  apply.  Therefore, 
existing  LEFM  methods  do  not  have  to  be  modified  to  account 
for  the  "short  crack  effect"  or  threshold  stress  intensity 
for  crack  propagation.  Since  an  empirically-based  crack 
growth  model  is  used,  mechanistic  theories  for  microstructur- 
al  behavior  are  not  required.  Since  the  ADA  methodology  is 
intended  as  a  "durability  design  tool,"  practical  considera¬ 
tions  for  achieving  reasonable  durability  analysis  results 
are  stessed  rather  than  mechanistic  theories  for  microstruc- 
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tural  behavior.  This  philosophy  in  no  way  diminishes  ths  im¬ 
portance  of  continuing  research  in  microstructural  behavior. 
Such  research  is  important  to  develop  a  fundamental  under¬ 
standing  of  microstructural  behavior  and  appropriate  mech¬ 
anistic-based  theories. 

7.  A  step-by-step  data  pooling  procedure,  based  on  the 
combined  least  square  sums  approach  (CLSSA) ,  has  been  devel¬ 
oped  for  estimating  the  EIFSD  parameters  using  one  or  more 
fractographic  data  sets.  A  scaling  method  is  provided  for 
normalizing  the  fractographic  results  for  multi-hole  test 
specimen  to  an  equivalent  single  hole  basis.  Therefore, 
fractographic  results  for  different  data  sets  can  be  "mixed 
and  matched"  to  obtain  the  EIFSD  parameters.  Data  pooling 
increases  the  size  of  the  data  base.  Since  EIFSD  parameters 
are  estimated  in  a  "global  sense,"  the  resulting  EIFSD  can  be 
justified  for  more  general  durability  analysis  applications. 
The  EIFSD  parameters  can  be  optimized  for  a  given  fracto¬ 
graphic  data  set (s)  by  minimizing  the  total  standard  error. 

8.  EIFSD  parameters  for  the  functions  considered  (i.e., 
Weibull  compatible,  lognormal  compatible,  lognormal  and  two- 
parameter  Weibull)  can  be  estimated  using  either  a  "TTCI  fit" 
or  "EIFS  fit."  It  doesn't  matter  which  fit  is  used  to  deter¬ 
mine  the  EIFSD  parameters  because  the  same  results  are  ob¬ 
tained  anyway.  However,  from  a  philosophical  standpoint,  the 
"TTCI  fit"  may  be  preferred  because  TTCIs  are  physical  quan¬ 
tities  that  can  be  obtained  directly  from  the  fractographic 
results.  On  the  other  hand,  the  EIFSs  are  artificial  flaws 
derived  from  the  TTCIs  using  a  suitable  crack  size-time  rela¬ 
tionship.  For  every  TTCI  there  is  a  corresponding  EIFS. 

9.  Each  fractographic  data  set  should  ba  screened  and 
plotted  before  the  data  set  is  used  to  estimate  the  EIFSD  pa¬ 
rameters  for  a  given  EIFSD  function.  Screening  is  recommend¬ 
ed  to  check  the  consistency  of  the  fractographic  results  so 
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that  questionable  results  can  be  censored  out.  There  are  no 
"hard  and  fast"  rules  for  censoring  but  une^  ained  extreme 
values  should  be  scrutinized. 

10.  This  program  has  further  advanced  the  understanding 
(e.g.,  meaning  and  limitations)  of  an  EIFS  and  EIFSD  for  dur¬ 
ability  analysis  applications.  The  durability  analysis  meth¬ 
odology  for  determining  a  suitable  EIFSD ,  including  optimiza¬ 
tion  of  the  EIFSD  parameters,  has  been  advanced  Therefore, 
confidence  in  the  EIFSD  for  durability  analysi*  has  been  in¬ 
creased.  Durability  analysis  methods  have  also  been  develop¬ 
ed  and  verified  for  covering  both  the  small  and  large  crack 
size  regions. 

11.  Methods  have  been  developed  for  determining  the 
EIFSD  and  initial  flaw  size  for  a  selected  probability  of 
crack  exceedance.  Such  methods  could  be  used  to  standardize 
the  way  aerospace  contractors  define  initial  flaw  size  from 
fractographic  results  for  teer-deem  inspections  and  other 
fatigue  tests.  A  standard  method  for  detinieeg  initial  flaw 
size  from  fractographic  results  is  needed  *e  that  comparable 
initial  flaw  sizes  can  be  obtained  by  different  aerospace 
contractors. 

12.  The  DCGA  is  recommended  for  durability  analyses  in 
the  small  crack  size  region  (e.g. ,  <  0.10*  and  the  EIFSD 
should  be  defined  using  the  DCGA.  A  fractographic  crack  size 
range  of  AL-AU  -  0.01"  -  0.05"  is  recommended  for  defining 
the  EIFSD.  The  reference  crack  siae  for  TTCIs,  x1 ,  should 
lie  in  the  AL-AU  range  (i.e.,  AL  <  x^  <  AU) .  A  different 
EIFSD  could  be  specified  for  applications  to  different  crack 
size  ranges.  However,  a  single  EIFSD  for  all  crack  sizes  is 
preferred.  If  a  single  EIFSD  for  any  crack  size  is  used,  the 
service  crack  growth  master  curve  (SCGMC)  for  p(i,T)  and 
FT(t)  predictions  should  be  approximated  by  at  least  two 
curve  segments.  This  aspect  will  be  discussed  later.  In  the 
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small  crack  sizs  region  a  single  segment  DCGA  SCGMC  can  be 
used. 


13.  A  two-segment  approach  for  durability  analysis  is 
recommended  to  cover  the  small  and  large  crack  size  region. 
The  cumulative  distribution  of  TTCI,  FT(t),  can  be  obtained 
for  a  convenient  reference  crack  size  in  the  AL-AU  range,  by 
growing  the  EIFSO  forward  using  a  single  segment  DCGA  SCGMC. 
For  example,  x1  -  0.05"  would  be  a  reasonable  reference  crack 
size  to  use  for  TTCIs  for  two  reasons:  (1)  TTCIs  for  x^  » 
0.05"  can  be  directly  verified  by  fractographic  results  and 
(2)  an  0.05"  crack  size  in  a  fastener  hole  is  a  reasonable 
maximum  crack  size  that  can  be  cleaned  up  by  reaming  the  hole 
to  the  next  fastener  size.  The  two  segment  approach  is  at¬ 
tractive  because  the  total  fatigue  life  can  be  divided  into 
two  parts,  i.e.,  crack  initiation  and  crack  propagation. 
Crack  initiation  could  be  used  for  small  cracks  (e.g., 
<0.05")  and  crack  propagation  could  be  used  for  cracks  > 
0.05".  The  DCGA  and  the  stochastic  crack  growth  approach 
(SCGA)  can  be  used  to  evaluate  functional  impairment  due  to 
fuel  leakage/ligament  breakage.  Reasonable  durability  anal¬ 
ysis  predictions  were  obtained  in  the  large  crack  size  region 
(e.g.,  0.50"  -  0.75"  through-  the-thickness  cracks)  using  the 
two  segment  DCGA.  However,  based  on  results  presented  in 
Volume  II  [3],  the  two-segment  DCGA-SCGA,  in  which  the  first 
segment  (or  crack  initiation)  is  based  on  the  DCGA  and  the 
second  segment  (crack  propagation)  is  based  on  the  SCGA,  is 
recommended  for  durability  analysis  in  the  large  crack  size 
region. 


14.  A  simple  crack  growth  rate  model  (da(t)/dt  ■  Qa(t)) 
was  used  to  develop  crack  size-time  relationships  between 
EIFSs  and  TTCIs  and  to  transform  TTCID  functions  into  a  EIFSD 
function  and  vice  versa.  This  model  has  been  found  to  be 
very  reasonable  for  durability  analysis  applications.  Other 
crack  growth  models  could  also  be  used  to  define  EIFSs  and/or 
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EIFSD  functions.  Ths  crack  growth  model  ussd  should  provids 
ths  propsr  crack  size-time  trends  for  a  wide  range  of  fracto- 
graphic  results.  Also,  the  model  parameters  for  different 
fractographic  data  sets  should  be  determined  in  such  a  way 
that  the  parameters  for  different  data  sets  are  directly  com¬ 
parable  and  that  they  are  on  the  same  baseline.  The  compat¬ 
ible  type  EIFSD  function  (i.e.,  EIFSD  function  is  compatible 
with  the  TTCID  function)  must  be  consistent  with  the  crack 
growth  model  used.  Furthermore,  the  durability  analysis  pre¬ 
dictions  in  the  small  crack  size  region  must  be  consistent 
with  the  basis  for  the  EIFSD.  These  guidelines  emphasize  the 
generality  and  flexibility  of  the  ADA  methodology  developed. 

15.  Compatible  and  non-compatible  type  EIFSD  functions 
were  considered  under  this  program.  The  "compatible"  EIFSD 
function  is  derived  from  the  TTCID  function  using  a  suitable 
deterministic  crack  growth  law.  The  EIFSD  function  can  also 
be  assumed  with  no  consideration  given  to  the  TTCID  function 
whether  it  is  physically  meaningful  or  not.  This  type  of 
EIFSD  function  is  referred  to  as  "non-compatible".  The  typ¬ 
ical  non-compatible  type  EIFSD  functions  (e.g.,  two-  para¬ 
meter  Weibull  and  lognormal)  used  will  generally  have  an 
"open"  upper  tail  which  implies  that  infinite  initial  flaw 
sizes  are  possible  for  some  probability  of  crack  exceedance. 
On  the  other  hand,  a  compatible  type  EIFSD  function  has  a 
"closed  upper  end"  which  imposes  an  upper  bound  limit  on  the 
initial  crack  size.  Compatible  type  EIFSD  functions  are  re¬ 
commended  for  durability  analysis  (e.g.,  Weibull  compatible 
and  the  lognormal  compatible  EIFSD  functions) . 

16.  An  EIFS  upper  bound  limit,  xu,  has  to  be  specified 
for  a  compatible  type  EIFSD  function.  This  limit  should  be 
established  in  conjunction  with  the  other  parameters  in  the 
EIFSD  function  (e.g.,  oc  and  ^  for  the  Weibull  compatible 
EIFSD  function)  and  the  economical  repair  limit  for  fastener 
holes.  The  EIFSD  parameters  should  be  optimized  to  minimize 
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the  error  in  tha  FT(t)  and  p(i,7)  predict iona  in  the  araaa  of 
most  interest.  For  example,  for  FT(t)  predictions  the  area 
of  most  interst  is  the  tail  of  the  TTCID  where  TTCIs  are  the 
smallest.  For  p(i,T)  predictions  the  upper  tail  of  the  crack 
size  distribution  where  crack  sizes  are  the  largest  is  of 
most  interest.  It  has  been  determined  in  Volume  II  [3]  that 
different  EIFSD  parameter  values  are  obtained  for  a  given  xu< 
The  accuracy  of  FT(t)  and  p(i,T)  predictions  seem  to  be  in¬ 
sensitive  to  significant  changes  in  the  xu  value.  Conse¬ 
quently,  reasonable  predictions  for  FT(t)  and  p(i ,T)  were  ob¬ 
tained  using  different  xu  values.  For  example,  different 
values  of  oC  and  4>  can  be  obtained  for  the  Weibull  compatible 
EIFSD  function  for  xu  -  0.01"  and  0.03".  In  this  case,  one 
xu  value  is  three  times  larger  than  the  other  xu  value  -  yet, 
both  xu,  OC  and  £  combinations  can  still  provide  comparable 
durability  analysis  results.  An  EIFS  upper  bound  limit  range 
of  xu  »  0.03"  -  0.05"  is  considered  reasonable  for  durability 
analysis.  This  xu  limit  is  based  on  NDI ,  considerations,  the 
economical  repair  limit  for  fastener  holes,  and  the  initial 
flaw  size  for  damage  tolerance  requirements.  Ideally,  the 
probability  of  exceeding  the  economical  repair  limit  should 
be  zero  at  time  zero  (i.e.,  as -manufactured  state).  However, 
if  xu  is  greater  than  0.05",  the  probability  of  exceeding  the 
economical  repair  limit  will  not  be  zero. 

17.  "User  friendly"  computer  software  has  been  develop¬ 
ed  for  an  IBM  compatible  PC  for  implementing  the  durability 
analysis  methods  developed  under  this  program.  This  includes 
software  for:  (1)  storing  fractographic  data  on  floppy  disk, 
(2)  screening  and  plotting  fractographic  results,  (3)  opti¬ 
mizing  the  EIFSD  parameters  for  one  or  more  fractographic 
data  sets,  (4)  checking  and  plotting  goodness-of-fit  of  FT(t) 
and  p(i,T)  predictions  based  on  given  EIFSD,  (5)  making  dur¬ 
ability  analysis  predictions  based  on  DCGA  and/or  SCGA.  A 
software  user's  guide  is  available  in  Volume  V  [6]. 
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4 . 2  RECOMMENDATIONS 


The  following  recommendations  are  based  on  the  work  per¬ 
formed  under  this  program: 

1.  The  "durability  analysis  tools"  need  to  be  incorpor¬ 
ated  into  the  overall  durability  design  process  for  metallic 
aircraft  structures  so  that  aerospace  contractors  can  effect¬ 
ively  use  the  tools  developed. 

2.  The  durability  analysis  methodology  should  be  ex¬ 
tended  to  interference  fit  and  cold  worked  fastener  holes  and 
other  structural  details  such  as  cutouts,  lugs,  fillets,  etc. 
Appropriate  test  specimens  should  be  designed  and  demonstrat¬ 
ed  for  acquiring  the  fractographic  data  needed  to  define  the 
EIFSD  for  other  structural  details.  Guidelines  and  procedur¬ 
es  are  needed  for  applications. 

3.  Further  experience  is  needed  in  developing  service 
crack  growth  master  curves,  compatible  with  the  basis  for  the 
EIFSD  used,  using  an  analytical  crack  growth  program  [e.g., 
38,39].  Specifically,  the  analytical  crack  growth  program 
used  should  be  tuned  or  "curve  fitted"  to  the  EIFS  master 
curves  for  selected  fractographic  data  sets.  The  "tuned"  an¬ 
alytical  crack  growth  program  should  be  used  to  make  crack 
exceedance  predictions  for  different  durability  design  condi¬ 
tions.  Then,  the  predictions  should  be  verified  using  actual 
fractographic  results.  More  experience  is  needed  in  setting 
up  the  crack  growth  analysis  for  durability  applications,  in¬ 
cluding  the  applicable  assumptions  that  go  into  a  damage  tol¬ 
erance  analysis.  Increased  confidence  in  the  durability  an¬ 
alysis  predictions  for  a  wide  range  of  design  conditions  is 
needed. 


4.  The  ADA  methodology  should  be  used  to  investigate  an 
existing  aircraft  structure.  Durability  design  trade-offs 
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such  as  fores  structural  maintenance  requirements,  risk  of 
functional  impairment  due  to  excessive  cracking  and  fuel 
leaks  or  ligament  breakage,  life  cycle  costs,  aircraft  readi¬ 
ness,  structural  weight  and  life  should  be  quantified  and  the 
trade-offs  evaluated  and  discussed. 

5.  The  initial  flaw  concepts  and  methods  developed  un¬ 
der  this  program  for  durability  analysis  should  also  be  in¬ 
vestigated  and  evaluated  for  damage  tolerance,  risk  analysis, 
and  force  structural  maintenance  applications. 

6.  Probability  of  crack  detection  aspects  should  be  in¬ 
corporated  into  the  method  for  defining  an  EIFSD  for  struc¬ 
tural  details. 

7.  The  effects  of  crack  growth  interactions  on  durabil¬ 
ity  analysis  predictions  for  p(i,T)  and  FT(t)  should  be  in¬ 
vestigated  -  particularly  for  large  cracks  in  adjacent  fast¬ 
ener  holes  or  other  structural  details.  Durability  analysis 
research  so  far  has  not  addressed  this  issue.  As  long  as  the 
dominant  fatigue  crack  in  adjacent  fastener  holes  is  rela¬ 
tively  small  (e.g.,  <  0.05”),  the  growth  of  one  fatigue  crack 
should  not  significantly  effect  the  growth  of  a  neighboring 
crack  and  vice  versa. 

8.  The  effect  of  bolt  load  transfer  on  the  EIFSD  para¬ 
meters  and  the  service  crack  growth  master  curve  need  to  be 
investigated  further.  Fatigue  test  results  should  be  acquir¬ 
ed  using  test  specimens  where  the  amount  of  bolt  load  trans¬ 
fer  is  controlled.  The  double-reverse  dog  bone  specimens 
used  under  thi?  program  and  the  previous  program  had  a  vari¬ 
able  amount  of  bolt  load  transfer  due  to  the  nature  of  the 
specimen  design.  Therefore,  the  effects  of  bolt  load  trans¬ 
fer  on  the  EIFSD  parameters,  initial  flaw  size,  and  SCGMC 
could  not  be  properly  evaluated.  The  effect  of  the  %  bolt 
load  transfer  can  be  accounted  for  in  the  LEFH  crack  growth 
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program  used.  However,  further  work  is  needed  to  tune  the 
SCGMC  to  the  same  basis  as  the  EIFSO  with  the  %  bolt  load 
transfer  properly  accounted  for. 

9.  Fatigue  test  results,  based  on  test  specimens  with 
a  preflawed  fastener  hole,  should  be  investigated  for  use  in 
tuning  the  analytical  crack  growth  program.  More  experience 
is  needed  in  developing  the  SCGMC  for  specified  design  condi¬ 
tions  in  conjunction  with  the  EIFSD. 

10.  Tear  down  inspection  results  from  various  aircraft 
should  be  used  to  determine  the  EIFSD  for  different  types  of 
materials  and  fastener  hole  types.  EIFSDs  should  be  deter¬ 
mined  from  applicable  fractographic  results.  Initial  flaw 
sizes  should  be  determined  for  selected  crack  exceedances. 
Then,  the  initial  flaw  sizes  for  different  aircraft  and  mate¬ 
rials  should  be  compared.  The  initial  flaw  sizes  should  be 
developed  using  the  standard  procedures  described  in  this  re¬ 
port  so  that  the  initial  flaw  sizes  will  be  compatible. 

11.  The  ADA  methodology  should  be  extended  to  advanced 
composites  applications. 
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1 .  Combined  Least  Square  Sums  Approach  (CLSSA)  -  the 
least  square  sums  for  individual  fractographic  data  sets  are 
combined  to  estimate  the  EIFSD  parameters  in  a  "global 
sense."  This  approach  is  used  in  conjunction  with  tho  data 
pooling  philosophy. 

2 .  Compatible  Equivalent  Initial  Flaw  Size  Distribution 
Function  -  this  is  a  distribution  function  for  equivalent  in¬ 
itial  flaw  sizes  (EIFS)  which  is  derived  using  a  physically 
meaningful  cumulative  distribution  of  time-to-crack  initia¬ 
tion  (TTCI)  function  and  a  suitable  deterministic  crack 
growth  law. 

3 .  Crack  Size  -  is  the  length  of  a  crack  in  a  structur¬ 
al  detail  in  the  direction  of  crack  propagation. 

4.  Data  Pooling  -  is  a  concept  for  estimating  the  EIFSD 
parameters  for  a  given  EIFSD  function  for  one  or  more  fracto¬ 
graphic  data  sets  in  a  "global  sense."  A  data  pooling  proce¬ 
dure  is  used  to  justify  the  EIFSD  for  more  general  durability 
analysis  applications. 

5.  Deterministic  Crack  Growth  Approach  (DCGA)  -  Crack 
growth  parameters  are  treated  as  deterministic  values  result¬ 
ing  in  a  single  value  prediction  for  crack  length. 

6.  Durability  -  is  a  quantitative  measure  of  the  air¬ 
frame's  resistance  to  fatigue  cracking  under  specified  ser¬ 
vice  conditions.  Structural  durability  is  concerned  with  the 
prevention  of  functional  impairments  due  to:  (1)  excessive 
cracking  and  (2)  fuel  leakage/ligament  breakage.  Excessive 
cracking  is  concerned  with  relatively  small  subcritical  crack 
sizes  (e.g.,  <  0.05")  which  affect  functional  impairment, 
structural  maintenance  requirement  and  life-cycle-costs. 
Such  cracks  may  not  pose  an  immediate  safety  problem.  How¬ 
ever,  if  the  structural  details  containing  such  cracks  are 
not  repaired,  economical  repairs  cannot  be  made  when  these 
cracks  exceed  a  limiting  crack  size.  Functional  impairment 
due  to  fuel  leakage/ligament  breakage  is  typically  concerned 
with  large  through-the-thickness  cracks  (e.g.,  0.50"  - 
0.75").  Although  such  cracks  are  usually  subcritical,  they 
affect  the  residual  strength  and  may  required  increased  main¬ 
tenance  action. 

7.  Durability  Analysis  -  is  concerned  with  quantifying 
the  extent  of  structural  damage  due  to  fatigue  cracking  for 
structural  details  (e.g.,  fastener  hole,  fillet,  cutout,  lug, 
etc.)  as  a  function  of  service  time.  Results  are  used  to 
ensure  design  compliance  with  Air  Force's  durability  design 
requirements . 
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3.  Economic  Llfs  -  is  that  point  in  tima  whan  an  air¬ 
craft  structure ' a  damage  atata  dua  to  fatigue,  accidantal 
damage  and/or  anvironmantal  datarioration  raachaa  a  point 
whara  oparational  raadinaaa  goala  cannot  ba  praaarvad  by 
aconomically  accaptabla  maintananca  action. 

9.  Economic  Li fa  Critaria  -  ara  guidalinaa  and  formats 
for  defining  quantitative  aconoaic  lifa  raquiramanta  for  air¬ 
craft  atructura  to  satisfy  U.  S.  Air  Fores  Durability  dasign 
requirements.  Tha  aconomic  lifa  critarion  providas  tha  basis 
for  analytically  and  axparimantally  ansuring  dasign  com¬ 
pliance  of  aircraft  structure  with  durability  dasign  raquira¬ 
mants.  Two  racommandad  formats  for  aconomic  lifa  critaria 
are: 


o  probability  of  crack  exceedance 

o  cost  ratio:  repair  cost/replacement  cost 

10.  Economic  Repair  Limit  -  is  tha  maximum  damaga  size 
that  can  be  economically  repaired  (e.g.,  repair  0.03"  - 
0.05"  radial  crack  in  fastener  holes  by  reaming  hole  to  next 
size)  . 


11.  Equivalent  Initial  Flaw  Size  (EIFS)  -  is  an  artifi¬ 
cial  crack  size  which  results  in  an”  actual  crack  size  at  an 
actual  point  in  time  when  tha  initial  flaw  is  grown  forward. 
It  is  determined  by  back-extrapolating  fractographic  results. 
It  has  the  following  characteristics:  (1)  an  EIFS  is  an  ar¬ 
tificial  crack  assumed  to  represent  tha  initial  fatigue  qual¬ 
ity  of  a  structural  detail  in  the  as-manufactured  condition 
whatever  the  source  of  fatigue  cracking  may  be,  (2)  no  direct 
relationship  to  actual  initial  flaws  in  fastener  holes  such 
as  scratches,  burrs,  microdefects,  etc.,  and  it  cannot  be 
verified  by  NDI ,  (3)  a  universal  crack  shape  in  which  the 
crack  size  is  measured  in  the  direction  of  crack  propagation, 
(4)  it's  in  a  fracture  mechanics  format  but  EIFSs  are  not 
based  on  linear  elastic  fracture  mechanics  (LEFM)  principles 
and  are  not  subject  to  laws  or  limitations,  such  as  tha 
"short,  crack  effect"  [e.g.,  31-37],  (5)  it  depends  on  the 
fractographic  data,  the  fractographic  crack  size  range  for 
the  back-extrapolation,  and  tha  crack  growth  rate  model  used, 
(6)  it  must  be  grown  forward  in  a  manner  consistent  with  the 
basis  for  the  EIFS,  (7)  EIFSs  are  not  unique  -  a  different 
sat  is  obtained  for  each  crack  growth  law  used  for  the  back- 
extrapolation,  and  (8)  EIFSs  are  not  necessarily  comparable 
with  minimum  crack  sizes  obtained  from  fractographic  results 
from  tear-down  inspections. 


12.  Equivalent  Initial  Flaw  Size  Distribution — _ 
is  used  to  represent  the  initial  fatigue quality  variation  of 
a  structural  detail.  An  EIFS  is  a  random  variable,  and  the 
EIFSD  statistically  describes  the  EIFS  population. 

13.  EIFS  Master  Curve  -  is  a  curve  (e.g.,  equation, 
tabulation  of — alt)  vi7“t  or  curve  without  prescribed  func¬ 
tional  form)  used  to  determine  the  EIFS  value  at  t-0  corres¬ 
ponding  to  a  given  TTCI  value  at  a  specified  crack  size. 
Such  a  curve  is  needed  to  determine  the  IFQ  distribution  from 
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the  TTCI  distribution.  Ths  EIFS  master  curvs  depends  on  sev¬ 
eral  factors,  such  as  ths  fractographic  data  bass,  ths  frac- 
tographic  crack  sizs  rangs  ussd,  ths  functional  fora  of  ths 
crack  growth  aquation  ussd  in  ths  curvs  fit,  stc.  (Raf. 
EIFS)  . 


14.  Extant  of  Damage  -  is  a  guantitativs  measure  of 
structural  durability  at  a  given  ssrvics  tias.  For  sxaapls, 
ths  nuabar  of  structural  dstails  (a.g.,  fastsnsr  holss,  cut¬ 
outs,  fillsts,  stc.)  or  psrcsntags  of  dstails  sxcssding 
spscifisd  crack  sizs  liaits.  Crack  lsngth  is  ths  fundaasntal 
asasurs  for  structural  daaags.  Ths  prsdictsd  sxtsnt  of  dam- 
ags  is  coaparsd  with  ths  spscifisd  sconoaic  lifs  critsrion 
for  snsuring  dssign  coapliancs  with  U.  S.  Air  Fores  durabil¬ 
ity  rsquirsasnts . 

15.  Genarlc  EIFS  Distribution  -  An  EIFS  distribution  is 
"generic"  if  it  depends  only  on  ths  material  and  aanufactur- 
ing/ fabrication  processes.  An  EIFSD  is  not  strictly  "gen¬ 
eric"  because  it  is  based  on  fractographic  results  which  re¬ 
flect  given  conditions  (e.g.,  load  spectra,  stress  level,  % 
bolt  load  transfer,  etc.).  For  durability  analysis,  an  EIFSD 
is  established  using  the  fractographic  results  for  one  or 
more  data  sets,  and  the  resulting  EIFSD  is  justified  for  a 
different  set  of  conditions. 

16.  Initial  Fatigue  Quality  (IFQ)  -  characterizes  the 
initial  manufactured  state  of  a  structural  detail  or  details 
with  respect  to  initial  flaws  in  a  part,  component,  or  air¬ 
frame  prior  to  service.  Actual  initial  flaws  in  a  fastener 
hole  are  typically  random  scratches,  burrs,  microscopic  im¬ 
perfections,  etc.  Such  flaws  are  not  cracks  per  se  like 
those  associated  with  linear  elastic  fracture  mechanics.  The 
IFQ  is  represented  by  an  equivalent  initial  flaw  size  distri¬ 
bution  (EIFSD) . 

17 .  Probability  of  Crack  Exceedance  (p(i,T))  -  refers 
to  the  probability  of  exceeding  a  specified  crack  size  , 

at  a  given  service  time,  T .  It  can  be  determined  from  the 
statistical  distribution  of  crack  sizes  and  can  be  used  to 
quantify  the  extent  of  damage  due  to  fatigue  cracking  in 
fastener  holes,  cutouts,  fillets,  lugs,  etc. 

18*  Reference  Crack  Size  (a  J_  -  This  is  the  specified 
crack  size  In  a  detail  used  to  reference  TTCI's.  The  IFQ 
distribution  'is  based  on  a  selected  reference  crack  size. 

19.  Service  Crack  Growth  Master  Curve  (SCGMC)  -  This 
curve  is  used  to  determine  the  EIFS,  y.WT),  corresponding  to 
an  exceedance  crack  size  x,  at  timeV'.  The  probability  of 
crack  exceedance,  p(i,T),  can  be  determined  from  the  EIFS 
cumulative  distribution  for  a  given  y, . ( T  ) .  The  SCGMC  is 
defined  for  the  applicable  design  variables  (e.g.,  stress 
level,  spectrum,  etc.),  and  it  can  be  determined  using  either 
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test  data  or  an  analytical  crack  growth  program.  All  sCGMC's 
must  bs  consistent  with  ths  corresponding  EIFS  master  curve 
and  the  fractographic  data  base.  The  SCGMC  must  be  consis¬ 
tent  with  the  basis  for  the  IFQ  distribution. 

20.  Statistical  Scaling  -  is  used  to  normalize  the 
fatigue  cracking  resistance  or  all  structural  details  in  a 
test  specimen  and  fractographic  data  set(s)  to  an  equivalent 
single  detail  population  basis.  Statistical  scaling, 
minimizes  the  fractographic  reading  requirements  because 
fractographic  results  are  required  for  only  the  largest 
fatigue  crack  per  test  specimen. 

21.  Stochastic  Crack  Growth  Approach  (SCGA)  -  an  app¬ 
roach  which  directly  accounts  for  the  crack  growth  rate  dis¬ 
persion  in  the  durability  analysis. 

22.  Structural  Detail  -  is  any  element  in  a  metallic 
structure  susceptible  to  fatigue  cracking  (e.g.,  fastener 
hole,  fillet,  cutout,  lug,  etc.). 

23.  Time-To-Crack-Initiation  (TTCI)  -  is  the  time  or 
service  hours  required  to  initiate  a  specified  (observable) 
fatigue  crack  size,  a0  ,  in  a  structural  detail  (with  no  init¬ 
ial  flaws  intentionally  introduced) . 

24.  TTCI  Lower  Bound  Limit  (g )  -  is  a  cutoff  value 
for  TTCI ' s  reflected  in  the  IFQ  model.  It  varies  for  a  given 
a0  and  it  depends  on  the  EIFS  upper  bound  limit,  x  ,  and  the 
EIFS  master  curve.  TTCIs  for  a  given  crack  size,  xlT,  should  be 

6  .  This  Weibull  distribution  parameter  provides  a  basis 
for  quantifying  the  EIFS  distribution  for  different  TTCI 
crack  sizes  on  a  common  baseline. 

25.  Upper  Bound  EIFS  Limit  (x  )  -  defines  the  largest 
EIFS  in  the initial  fatigue  quality  distribution.  Con¬ 
straints  on  x  for  fatigue  holes:  largest  EIFS  in  data  set  < 
xu  (0.03"-0.05") . 


ACRONYNS 


ADA  -  Advanced  Durability  Analysis 

CLSSA  -  Combined  Least  Square  Suns  Approach 

DADTA  -  Durability  and  Damage  Tolerance  Assessment 

DCGA  -  Deterministic  Crack  Growth  Approach 

EIFS  -  Equivalent  Initial  Flaw  Size 

EIFSD  -  Equivalent  Initial  Flaw  Size  Distribution 

FHQ  ■  Fastener  Hole  Quality 

HEIFS  *  Homogeneous  EIFS 

IFQ  -  Initial  Fatigue  Quality 

LEFM  -  Linear  Elastic  Fracture  Mechanics 

LT  *  Load  Transfer  Through  the  Fastener 

MM  -  Method  of  Moments 

NDE  -  Non  Destructive  Evaluation 

NDI  -  Non  Destructive  Inspection 

NLT  ■  No  Load  Transfer  Through  the  Fastener 

SCGA  ■  Stochastic  Crack  Growth  Approach 

SSE  «  Sum  Squared  Error 

TSE  -  Total  Standard  Error 

TTCI  -  Time-to-Crack  Initiation 


LIST  07  SYMBOLS 


a 

*o 

»(0) 

a(0) 

*(t) 


-  Crack  Siza 

»  Reference  crack  aiza  for  given  TTCZa 
•  EZFS  -  Crack  aiza  at  t-0 

-  Maan  EZFS 

•*  Crack  ai<sa  at  any  service  time  t 


a(t) ,  a(t2) ,  a(t2) 

a  (T) 


Crack  aiza  at  time  t ,  t^  and  t2,  res¬ 
pectively 

Crack  aiza  at  aarvica  tina  T 


a  (T) 


Crack  siza  at  any  aarvica  tine  T 


AL,  AU  ■  Lowar  and  upper  bound  fractographic 

crack  size,  respectively,  used  to  de- 
fine  the  EZFSD  parameters.  Also  used 
in  conjunction  with  the  SCGMC  to  de¬ 
fine  crack  size  limits  for  the  small 
crack  siza  region. 

AU1  »  Upper  bound  crack  size  limit  for  the 

large  crack  size  region 


b,  Q  -  Crack  growth  parameters  in  the  equation 

da (t) a ^ jb.  Used  in  conjunction 
dt 

with  the  1FQ  model. 


b^Qj  ■  Crack  growth  rata  parameters  in  the 

aquation  da/dt  *  a )  associated  with 

the  one-segment  DCGA  or  1st  segment  of 
the  two  segment  approach. 
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B 

C 


Crack  growth  rata  parameters  in  the 
aquation  da/dt  -  Q3(a)fet  for  segment 
two  of  tha  two-segment  DCS A. 

-  of  Inf  or  In  (  a(0)  ) 

b  -  If  Uaad  in  conjunction  with  tha  XFQ 
modal  whan  tha  crack  growth  law, 

5&CfelMQ[a(t)  r°  ia  uaad  and  b  >  l.o. 
dt 


ci 


COV(x) 

COV(?) 


b^  -  1;  Uaad  in  conjunction  with  tha 
SCGMC  whan  da  ft) mq  £a(t)]bi  uaad. 
dt  i 


Tha  aubacript  "i"  rafara  to  tha  ith 
straaa  ragion. 


Coefficient  of  Variation  of  x 
Coefficient  of  Variation  of  F 


daft)  ■  Crack  growth  rate  as  a  function  of  time 

dt 

*a(0) ^  “  EIrs  probability  density  function  « 

dr.(0)  <*> 

dx 

VtJ  -  dFT(t) 

dt 


fx(u) 


Probability  density  function  of  X. 
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*»(0) 


(X) 


Vo> 


(X) 


ra^(0)  (xij} 


Fa(t) 


(x) 


V> 


(x) 


FT(t) 

*J 


-  TITS  cumulative  distribution  function 
for  an  "equivalent  single  hola  popula¬ 
tion.  " 

■  Cumulative  distribution  of  EXFS  bassd 
on  ths  largest  fatigue  crack  per  test 
specimen. 

“  Subscripted  notation  used  for  Fft  ^ (x) 
used  in  conjunction  with  data  pooling, 
where:  j  denotes  the  jth  crack  in  the 
ith  data  set. 

■  Cumulative  distribution  of  crack  size 
at  any  service  time  t  based  on  an 
"equivalent  single  hole  population" 
EXFSD. 

-  Cumulative  distribution  of  crack  size 
at  any  service  time  t  based  on  the 
largest  fatigue  crack  per  test  specimen 
basis  for  the  EIFSD. 

■  TTCI  cumulative  distribution  function 

■  Cumulative  distribution  of  minimum 
TTCIs  based  on  the  largest  fatigue 
crack  per  test  specimen. 


rT„(ti3> 


Subscripted  notation  used  for  FT^(t) 
used  in  conjunction  with  data  pooling, 
where:  j  «  jth  TTCI  value  in  the  ith 
data  set. 
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Predicted  cumulative  distribution  of 
minimum  TTCZ  based  on  a  aalaetad  xu  and 
initial  valuaa  of  0^  and  ^  for  tha 
Waibull  compatible  EIF5D  function. 

Used  in  con junction  vith  non- linear 
procedure  for  estimating  the  EIFSD  pa¬ 
rameters  . 

No.  of  fastener  holes  per  test  specimen 
in  which  only  the  largest  fatigue  crack 
per  specimen  is  included  in  the  fracto- 
graphic  data  set.  Used  in  conjunction 
with  the  "statistical  scaling"  techni¬ 
que. 

Same  as  £  with  the  subscript  " j"  re¬ 
ferring  to  the  ith  fractographic  data 
set. 

Total  and  average  number  of  details, 
respectively,  in  the  entire  component 
having  a  crack  size  >x1  at  any  service 
time  T 

Load  transfer  through  the  fastener 

Total  number  of  fractographic  data  sets 
used  to  estimate  the  EIFSD  parameters 

Number  of  TTCI  or  EIFS  values  for  the 
ith  data  set  used  in  conjunction  with 
the  combined  least  square  sums  app¬ 
roach/  or  total  number  of  details  in  the 
ith  stress  region. 


N(i,r)#  N(i/n 


P(irT) 


Q 


Q 


i 


S  (?) 
S(2) 


-  Total  and  avaraga  number  of  details, 
respectively,  having  a  crack  air a  ax- 
caading  x1  at  any  aarvica  time  T 

-  Probability  that  a  datail  in  tha  ith 
atraaa  ragion  will  hava  a  crack  size 
>x^  at  tha  aarvica  tima  T 

-  Crack  growth  rata  parameter  in  da(t)/dt 
■  Qa(t)  rafarrad  to  aa  "pooled  Q"  valua 
whan  parameter  represents  tha  typical 
crack  growth  for  a  given  fractographic 
data  eat. 

•  Crack  growth  rata  parameter  in  da(t)/dt 
-  Q.a(t).  In  conjunction  with  fracto¬ 
graphic  data  pooling,  danotaa  tha 
crack  growth  rata  parameter  for  tha  ith 
data  aet  in  M  data  seta  to  be  used  to 
estimate  the  EIFSD  parameters.  In  con¬ 
junction  with  SCGA-based  EIFSs,  Q(-  is 
the  crack  growth  rate  parameter  for  the 
ith  fatigue  crack  in  a  data  set. 

«  Standard  deviation  of  ft 

«  Standard  deviation  of  x 

»  Flight  hours  at  t,  t^ ,  t2,  respective¬ 
ly. 


t^  -  Time  to  initiate  crack  size  x^ 

T,  TTCI  ■  Time-to-crack- initiation 
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A  particular  value  of  X 


u 


x 


X 


1 


X 


u 


X 


yli(r) 


Y 


ij 


•  crack  size 

*  Crack  size  used  for  p(i,7*)  pradictions 

«  Uppar  bound  limit  for  El FS 

■  Lognormal  random  variabla  with  a  madian 
of  1.0. 

-  In  ln(xy/x^j ) ;  EIFS  Pit 

-  ln[ln(xu/x1)  +  Q.t^j,  TTCI  m 
.  In  ln(xu/x^)  -  In  Qj  ;  EIFS  Pit 

-  In  [ln(*u/Xx)  +  Q|tij]  -  •" Q(.  7  TTCI  Fit 


An  EIPS  in  tha  EIFSD  corresponding  to  a 
crack  size  x^  at  time  T  in  the  ith 
stress  region.  Value  determined  using 
the  SCGMC. 


in 


-  (1/^  In  Fo  (0)(x  jj ) 
1 


7  EIFS  Fit 


in 


-  (1/i.)  in 


1  -  ft,  V 

i 


;  ttci  Fit 


-j 


1- 

B^O)  y  jj 

7  EIFS  Fit 


$ 


-1 


1  - 


..  ut 

(/..)!  ‘  ;  TTCI  Fit 

Te‘V J 


-  4,1.  ezS  (-  )  [ln{XJX^ )  + 
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-«t-l 


#1 


zij 


n  ) 


<r 


-  Log  X 

-  3/(H.  +  l)-i  +  •npl-jSj^') 

■  Weibull  distribution  parameters  for 
shape,  scala,  and  lowar  bound  TTCX, 
raspactivaly. 

■  Gamma  function 

-  Empirical  constants  in  the  equation: 

Q,  ■  where  C  *  stress 

■  Stress  or  standard  deviation 


<5 


Standard  deviation  of  da(t)/dt  with 
respect  to  plot  of  jB*?da(t)/dt  -  /r»a(t)  +  Qt 


T 


A  particular  service  time 


oC 


k 


«,4> 


Initial  oC  value  for  the  kth  data  set 
used  in  conjunction  with  data  pooling 
and  Method  of  Moments 

Two-parameter  Weibull  distribution 
shape  and  scale  parameters  for  EXFS, 
a(0),  respectively. 

Weibull  compatible  shape  and  scale 
EIFSD  parameters,  respectively 


etftfl4>o  ■  Initial  values  of  ot  and  ^  used  in  con¬ 

junction  with  the  non-linear  approach 
for  estimating  EIFSD  parameters  for  the 
Weibull  compatible  EIFSD  function. 
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Incramants  of  Of  and  &  ,  raspactivaly 

0  *0 


*,<r 


Maan  valua  and  standard  daviation  of 
In  a(0)  (or  In  x) ,  raspactivaly 


4>(  ) 


ij 


Standardizad  normal  distribution 
function 

k 

4  e.'tj 


V°2 


fa 

[1  -  <Qi/Q2)]toa0  j  b, 
An«o  *  Qi(c2Q2)*l(*o)‘C2>; 


bj  *  1 


3 ;  bj  *  i 


*  (eiM-'Uo)’61  -  (ejd2)-‘(,0)-eJ.  b 

s  1  m  *■>  b  *  1 
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APPENDIX  A 

DURABILITY  ANALYSIS  CONCEPTS ,  ISSUES  ,  AND  PHILOSOPHY 

A. 1  INTRODUCTION 

Durability  analysis  concepts,  key  issues  and  philosophy 
are  reviewed  and  discussed  in  this  section.  The  time-to- 
crack-initiation  (TTCI)  and  equivalent  initial  flaw  size 
(EIFS)  concepts  are  reviewed  in  terms  of  initial  fatigue 
quality  representation.  Key  issues  are  clarified  and  the 
probability  of  crack  exceedance  concept  is  discussed  in  terms 
of  deterministic  crack  growth. 

Initial  fatigue  quality  (XFQ) ,  the  cornerstone  of  any 
durability  analysis,  is  represented  by  an  equivalent  initial 
flaw  size  distribution  (EIFSD) .  An  "engineering  approach"  is 
used  to  represent  the  IFQ  of  the  structural  details  to  be  in¬ 
cluded  in  the  durability  analysis.  Although  engineering 
principles,  instead  of  mechanistic-based  theories,  are  used 
to  model  the  microcracking  behavior,  LEFM  principles  can 
still  be  used  to  perform  the  durability  analysis.  The  nature 
of  the  approach,  implications  and  perspective  for  durability 
analysis  need  to  be  reviewed  and  discussed.  Various  aspects 
are  discussed  based  on  the  advancements  and  understandings  of 
this  program  and  past  research.  A  conceptual  approach  is 
used  and  mathematical  details  are  provided  in  other  sections 
of  this  report  (Vol.  I). 


A. 2  TTCI  CONCEPT 

The  time  to  initiate  a  natural  fatigue  crack  in  a 
structural  detail  to  any  crack  size,  x1,  is  referred  to  as 
the  "time-to-crack-initiation"  (TTCI) .  Initial  fatigue 
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quality  of  structural  datails  can  ba  raprsaantad  using  aithar 
a  distribution  of  TTCI  or  aquivalant  initial  flaw  siza  (EIFS) 
[2/12/13,18].  TTCIs  can  ba  obtainad  from  fractographic  re- 
sults,  aithar  directly,  by  intarpolation  or  by  extrapolation. 
For  example,  in  Fig.  A.l  tha  TTCI  for  aach  of  tha  solid 
circlas  is  aqual  to  tha  corrasponding  tints  raquirad  to  init- 
iata  tha  spacific  crack  siza. 

A. 3  EIFS  CONCEPTS 

Actual  initial  flaws  in  tha  bora  of  manufactured 
fastener  holes  in  metallic  aircraft  structures  usually  con¬ 
sist  of  random  scratches,  burrs,  microscopic  imperfections, 
etc.  Under  service  conditions  fatigue  cracks  may  eventually 
initiate  from  such  flaws  and  then  propagate  to  crack  sizes 
that  can  be  detected  by  NDI .  However,  the  actual  flaws  in 
fastener  holes  produced  by  manufacturing  and  assembly  are  not 
typical  "cracks"  in  the  usual  sense  associated  with  the  LEFM 
approach.  Whatever  type  of  initial  flaw  may  be  in  a  fastener 
hole,  such  flaws  must  be  accounted  for  to  assure  that  the 
structure  will  be  durable  and  can  be  economically  maintained. 
Some  means  is  needed  to  quantify  the  initial  fatigue  quality 
of  actual  initial  flaws  in  such  a  way  that  their  effect  on 
structural  durability  can  be  accounted  for. 

An  equivalent  initial  flaw  is  an  artificial  crack  with  a 
characteristic  size  measured  in  the  direction  of  crack  propa¬ 
gation.  Such  a  flaw  is  not  an  actual  physical  flaw.  EIFSs 
are  determined  by  back-extrapolating  fractographic  results  to 
time  zero  (Fig.  A.l).  The  initial  fatigue  quality  is  repre¬ 
sented  by  an  equivalent  initial  flaw  size  distribution 
(EIFSD) .  This  is  an  "engineering  approach"  to  the  durability 
analysis  problem.  The  objective  of  such  an  approach  is  to 
provide  a  practical  "design  tool"  that  can  be  used  to  evalu¬ 
ate  Air  Force  durability  design  requirements  and  tradeoffs 
for  metallic  aircraft  structures. 


Notes  •-denotes  f ractographic  data  (i.e.,  a(t)#t) 


CRACK 
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EIFS=a (0) 


TIME 


Figure  A.l.  An  Equivalent  Initial  Flaw  is  Determined  b 


EIFSs  are  generally  dafined  using  either  a  deterministic 
crack  growth  approach  (DCGA)  or  a  stochastic  crack  growth  ap¬ 
proach  (SCGA) .  Tha  two  approaches  are  conceptually  described 
in  Fig.  A. 2  and  further  details  are  given  in  Appendix  E. 

In  Fig.  A. 2 ,  suppose  that  ElFSs  are  to  be  determined  for 
bo+'h  the  DCGA  and  tha  SCGA  using  the  same  fractographic  data. 
Using  the  DCGA,  the  EIFS  for  each  fatigue  crack  is  determined 
using  the  same  crack  si2a-time  relationship  or  "EIFS  master 
curve."  Each  EIFS  has  a  corresponding  time-to-crack-initia- 
tion  at  a  given  reference  crack  size  and  vice  versa.  As  il¬ 
lustrated  in  Fig.  A. 2 (a),  the  EIFSs  and  TTCIs  have  the  same 
rank  in  the  respective  populations.  With  the  DCGA,  the  vari¬ 
ance  in  the  crack  growth  rate,  is  directly  reflected  in  the 
EIFS  value  computed. 

Using  the  SCGA,  the  EIFS  for  each  fatigue  crack  is 
determined  using  a  crack  size-time  relationship  that  is  sep¬ 
arately  fitted  to  each  fatigue  crack.  In  other  words,  each 
fatigue  crack  is  back-extrapolated  to  time  t  «  0  using  a 
crack  size-time  relationship  that  is  tailored  for  each  fa¬ 
tigue  crack.  Whereas,  the  DCGA  uses  a  single  EIFS  master 
curve  to  determine  the  EIFSs,  the  SCGA  uses  a  different  EIFS 
curve  to  obtain  the  corresponding  EIFS  for  each  fatigue 
crack.  The  SCGA  is  conceptually  described  in  Fig.  A. 2(b). 
Unlike  the  DCGA,  the  crack  growth  trajectories  for  the  SCGA 
may  cross  over  each  other  as  depicted  in  Fig.  A. 2(b).  Since 
each  fatigue  crack  is  back-extrapolated  independently,  the 
crack  growth  rate  dispersion  is  filtered  out  of  the  EIFS 
computed. 

The  "basis"  for  an  EIFS  and/or  EIFSD  fracture  includes 
consideration  of  such  factors  as:  (1)  fractographic  crack 
size  range,  AL-AU,  used,  (2)  material,  (3)  fatigue  testing 
conditions  (e.g.,  spectrum,  stress  level,  %,  bolt  load 
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transfer,  etc.)/  (4)  manufacturing  process,  and  (5)  crack 
growth  law  and  procedure  for  back-extrapolation,  including 
goodness-of-fit  criterion.  An  EIFS  and/or  EIFSD  should  be 
used  in  the  proper  context  for  durability  analysis.  For 
example,  an  artificial  initial  flaw  size  or  EIFS  should  not 
be  indiscriminately  grown  forward  using  a  LEFM-basad  crack 
growth  program  without  considering  the  "basis*1  for  the  EIFS 
or  EIFSD.  If  this  principle  is  not  strictly  followed,  dur¬ 
ability  analysis  predictions  will  not  be  consistent  with  the 
IFQ,  and  the  results  will  be  questionable. 

An  EIFS  is  a  random  variable.  Each  structural  detail  to 
be  used  in  the  durability  analysis  is  assumed  to  have  a  sin¬ 
gle  dominant  flaw  which  can  be  represented  by  an  equivalent 
initial  flaw  size  (EIFS) .  The  population  of  EIFSs  is  repre¬ 
sented  by  an  equivalent  initial  flaw  size  distribution 
(EIFSD)  (Fig.  A. 3).  Initial  flaw  sizes  can  be  selected  from 
the  EIFSD  which  correspond  to  a  given  probability  of  crack 
exceedance.  However,  such  initial  flaws  are  artificial  -  not 
actual  physical  flaws.  For  example,  an  equivalent  initial 
flaw  from  an  EIFD,  for  a  given  probability  of  crack  exceed¬ 
ance  (e.g.,  1/10000  or  0.0001),  has  no  direct  relationship  to 
the  minimum  crack  size  obtained  from  the  fractographic  evalu¬ 
ation  of  tear-down  inspection  results  (Fig.  A. 4). 

The  EIFSD  reflects  a  fracture  mechanics  format.  How¬ 
ever,  if  the  EIFSD  is  based  on  an  empirical  crack  growth  law, 
the  growth  of  the  EIFSD  forward  must  be  consistent  with  the 
crack  growth  law  used.  A  LEFM  crack  growth  program  can  be 
used  to  grow  the  EIFSD  forward  but  only  in  terms  of  the  ap¬ 
plicable  crack  growth  law  reflected  in  the  EIFSD.  An  artifi¬ 
cial  initial  flaw  should  not  be  grown  forward  using  a  LEFM 
program  without  considering  the  basis  for  the  EIFS.  (Fig. 
A.  5)  . 


^Measured  in  Direction  of  Crack  Propagation 


STRUCTURAL  DETAILS 
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CRACK 
SIZE  * 
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Figure  A. 3. 


An  Equivalent  Initial  Flaw  Size  Distribution  Represents 
the  Initial  Fatigue  Quality  of  Structural  Details. 


TEAR-DOWN  INSPECTION 


a 


Fractographic  Results 


CRACK 


r«vi 


“  minimum  crack  size  from  fractographic  evaluation 


Figure  A. 4.  An  Artificial  Initial  Flaw  Size  is  Not  Necessaril 
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Figure  A. 5.  An  EIFSD  or  Initial  Flaw  Size  Must  be  Grown  Forward 
Consistent  with  the  Basis  for  the  EIFSD . 
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Sinc4  an  aquivalant  initial  flaw  ia  not  an  actual  physi¬ 
cal  flaw  which  can  ba  varifiad  by  NDE,  than  how  can  an  &X7SD 
ba  juatifiad  for  durability  analysis?  An  EIFSD  can  ba  justi- 
fiad  by  showing  that  raasonabla  pradictions  for  tha  cumula- 
tlva  distribution  of  crack  or  cumulative  distribution  of  TTCX 
can  ba  obtained  using  tha  given  EXFSD.  Such  pradictions  can 
ba  correlated  with  actual  fractographic  results  for  given 
conditions  (a.g.,  load  spectrum,  stress  level,  %  bolt  load 
transfer,  etc.)  (Fig.  A. 6). 

Tha  concept  of  a  single  dominant  flaw  in  a  fastener 
hole(s)  is  used.  (Fig.  A. 7).  This  concept  greatly  simpli¬ 
fies  the  durability  analysis  and  is  an  important  aspect  of 
the  overall  method. 

Fatigue  cracks  originate  in  fastener  holes.  Multiple 
micro-cracks  with  different  geometries  and  locations  may 
occur  in  fastener  holes.  During  microcrack  initiation  and 
crack  growth,  there  can  be  multiple  branching  with  complex 
interactions  between  the  various  cracks  affecting  the  crack 
growth  rates  and  damage  accumulation.  Xn  general,  the  micro¬ 
cracks  eventually  merge  into  a  single  dominant  crack.  The 
crack  growth  rates  for  the  dominant  crack  can  ba  determined 
from  the  fractographic  results  (Fig.  A. 7).  Each  fastener 
hole  or  holes  sharing  a  common  fastener  have  a  single  domi¬ 
nant  crack  in  the  fastener  hole  (Fig.  A. 7).  The  effects  of 
the  various  cracks  in  the  fastener  hole  on  the  crack  growth 
rate  is  accounted  for  in  the  experimental  data.  Therefore, 
the  effects  of  microcracking  and  crack  interactions  on  the 
crack  growth  rate  and  damage  accumulation  are  accounted  for. 
The  concept  of  a  single  dominant  flaw  in  a  fastener  hole 
greatly  simplifies  the  analysis.  Each  fastener  hole  is  as¬ 
sumed  to  have  a  single  dominant  flaw,  and  the  crack  growth  of 
the  dominant  flaw  in  each  fastener  hole  is  assumed  to  be  in¬ 
dependent  of  the  crack  growth  of  the  dominant  flaws  in  sur¬ 
rounding  fastener  holes. 
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Fatigue  C racks  in  Fastener  Holes  May  Have  Different 
Initiation  Sites  and  Exhibit  Crack  Branching. 


A. 4  CLARIFICATION  OF  EIFS  ISSUES 

The  following  issues  have  been  raised  about  the  EIFS 
approach  for  durability  analysis:  (1)  the  EIFS  approach  does 
not  account  for  limitations  of  LEFM  and  the  "short  crack 
effect,"  (2)  a  threshold  stress  intensity  may  exist  below 
which  cracks  do  not  propagate  (3)  an  EIFS  and/or  EIFSD  is  not 
a  "generic"  property  of  the  material  and  manufacturing  pro¬ 
cesses  (4)  EIFS  does  not  necessarily  account  for  crack  geome¬ 
try  and  shape.  These  issues  are  discussed  and  clarified  in 
this  section. 

Before  addressing  the  above  issues,  we  will  clarify  the 
use  of  an  EIFSD  for  representing  the  initial  fatigue  quality 
of  a  structural  detail.  There  are  two  basic  durability  de¬ 
sign  requirements:  (2)  excessive  cracking  and  (2)  fuel  leak¬ 
age/ligament  breakage.  The  excessive  cracking  requirements 
are  concerned  with  relatively  small  fatigue  cracks  (e.g., 
<0.05").  For  example,  to  predict  the  probability  of  exceed¬ 
ing  an  0.05"  crack  size  in  a  fastener  hole  at  any  time,  t,  a 
reasonable  estimate  of  the  initial  flaw  size  distribution  and 
a  means  for  growing  the  population  of  initial  flaws  forward 
is  needed. 

Ideally,  we  would  like  to  define  the  geometry,  shape  and 
location  of  actual  initial  flaws  in  the  bore  of  a  fastener 
hole  due  to  material  and  manufacturing  processes  using  con¬ 
ventional  NDI  methods.  Such  initial  flaws  are  too  small  to 
detect  and  accurately  quantify  in  a  physical  sense  by  direct 
NDI  means.  If  actual  initial  flaws  could  be  accurately  quan¬ 
tified  (i.e.,  location,  geometry,  shape)  in  fastener  holes, 
two  major  problems  would  still  have  to  be  resolved:  (1)  how 
to  accurately  characterize  actual  initial  flaws  mathematical¬ 
ly  and  (2)  how  to  predict  crack  initiation  by  such  flaws. 
Actual  initial  flaws  in  the  bore  of  a  fastener  hole  can  have 
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different  shapes,  geometry  and  location.  Moreover,  multiple 
micro-flaws  may  exist,  and  such  flaws  may  have  different 
crack  initiation  sites  under  service  conditions.  Further¬ 
more,  during  the  initial  stage  of  the  crack  initiation, 
growth  and  interaction  of  such  cracks  is  very  complex.  For 
example,  microcrack  branching  can  occur  in  which  the  growth 
of  one  crack  affects  the  growth  to  another  crack  and  vice 
versa.  With  multiple  crack  initiation  sites  and  crack 
branching,  cracks  tend  to  merge  into  a  common  crack  front 
(Fig.  A. 7) . 

A  suitable  mechanistic  model  for  characterizing  the 
crack  initiation  process  in  the  small  crack  size  region  is 
simply  not  available  to  implement  the  probabilistic-based 
durability  analysis  approach.  Therefore,  an  "engineering 
approach"  is  used  to  define  the  initial  fatigue  quality.  The 
"engineering  approach"  is  considered  to  be  a  "design-tool" 
for  evaluating  airframe  durability  design  requirement  and 
tradeoffs.  The  goal  of  the  "engineering  approach"  is  to  pro¬ 
vide  reasonable  durability  analysis  results.  Research  work 
should  be  continued  to  develop  a  better  understanding  of  the 
crack  initiation  process  and  an  effective  model.  In  the 
meantime,  an  "engineering  approach"  is  recommended  for  dura¬ 
bility  analysis. 

A. 4.1  LEFM  Limitations  and  "Short  Crack  Effect" 

It  is  well  known  that  linear  elastic  fracture  mechanics 
(LEFM)  laws  break  down  for  small  crack  sizes.  Although  the 
minimum  crack  size  for  valid  LEFM  applications  is  debatable, 
a  minimum  size  of  0.005"  [59]  is  considered  reasonable.  The 
crack  growth  rate  is  typically  faster  for  "short  cracks"  than 
for  "long  cracks"  for  the  same  stress  intensity  range  [e.g., 
31-37] . 
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The  EIFS  approach  and  philaophy  for  durability  analysis 
in  viaw  of  tha  abova  conaidarations  will  now  ba  axplainad  in 
a  concaptual  manner  using  Fig.  A. 8.  Thraa  diffarant  ap¬ 
proaches  are  illustrated  in  Fig.  A. 8  for  determining  the  EIFS 
value  using  actual  fractographic  results  in  a  selected  AL-AU 
crack  size  range.  In  each  of  the  three  cases  the  fracto¬ 
graphic  results  are  extrapolated  backwards  to  time  zero  to 
obtain  an  EIFS  value.  An  empirical-based  crack  growth  law  is 
assumed  for  case  1.  In  case  2,  it  is  assumed  that  the  LEFM 
program  is  simply  a  "curve  fitting  tool"  applicable  to  any 
crack  size.  In  this  case  the  "short  crack  effect"  is  ignored 
entirely.  In  case  3,  the  back-extrapolation  procedure 
accounts  for  both  the  "short  crack  effect"  and  limits  the  use 
of  LEFM  to  crack  sizes  where  most  consider  LEFM  methods  to  be 
valid.  It  is  assumed  in  case  3  that  the  mechanistic-based 
model  is  used  to  accurately  characterize  the  crack  initiation 
process  in  the  small  crack  size  region. 

In  Fig.  A. 8  it  is  assumed  that  the  empirical  crack 
growth  law  is  fitted  on  a  least  squares  fit  basis  to  the 
fractographic  results  in  the  AL-AU  range.  For  cases  2  and  3, 
it  is  assumed  that  a  trial  and  error  curve  fitting  method  is 
used  and  the  goodness-of-fit  is  based  on  a  subjective  deci¬ 
sion  by  the  analyst.  The  resulting  EIFS  value  for  cases  1, 
2,  and  3  is  denoted  in  Fig.  A. 8  by  x1#  x2,  and  x3,  respect¬ 
ively.  Observations  and  interpretation  of  the  results  de¬ 
picted  in  Fig.  A. 8  are:  (1)  a  different  EIFS  value  is  ob¬ 
tained  for  each  case  considered ,  (2)  using  the  fractographic 
results  for  a  given  AL-AU  range,  the  same  EIFS  value  for  case 
1  should  be  obtained  by  any  analyst  if  the  least  squares  fit 
criterion  and  the  same  fitting  procedures  are  used,  and  (3) 
since  the  goodness-of-fit  is  subjective  for  cases  2  and  3,  no 
two  analysts  will  necessarily  get  the  same  initial  flaw 
sizes. 
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Considered  as  "curve  fitting  tool"  to  justify  in  "short  crack"  region 


Figure  A. 8.  Equivalent  Initial  Flaw  Size  (EIFS)  Depends  on  Approach  Used 
to  Back-extrapolate  Fractographic  Results  In  A  Selected  AL^AU 


No  standard  analytical  crack  growth  program  is  used  by 
most  asrospaca  companias  today.  Most  companies  usually 
dsvslop  thsir  own  working  versions  for  fracturs  mechanics 
applications.  Soms  companies  may  try  to  incorporate  the 
'•short  crack  effect"  into  their  program  while  others  do  not. 
Different  analytical  crack  growth  programs  are  in  use  today. 

The  back-extrapolation  procedure  initially  developed  for 
determining  EIFSs  for  the  F-4  DADTA  program  [57]  is  based  on 
an  analytical  crack  growth  program  and  a  subjective  decision 
by  the  analyst  is  required  about  goodness-of-f it .  In  view  of 
the  different  versions  of  analytical  crack  growth  that  are 
currently  being  used  and  the  subjective  decision  factor,  it's 
not  surprising  that  EIFSs  for  the  same  material  and  manufact¬ 
uring  process  differ  so  much  from  company-to-company .  There 
has  been  a  real  need  for  some  time  now  to  develop  a  universal 
procedure  for  defining  EIFSs  and  EIFSDs  so  that  the  same  in¬ 
itial  fatigue  quality  results  will  be  obtained  by  any  analyst 
or  aerospace  company  for  the  same  fractographic  data  set. 

A. 4. 2  Consistent  EIFSs  and  Crack  Growth  Predictions 

In  any  durability  analysis,  it  is  very  important  that 
EIFSs  be  defined  in  a  consistent  manner.  Furthermore,  the 
EIFSs  must  be  grown  forward  in  the  same  manner  that  the  EIFSs 
were  determined.  If  this  principle  is  not  strictly  followed 
inconsistent  EIFSs  and  durability  analysis  results  will  be 
obtained.  To  conceptually  illustrate  this  crucial  point  con¬ 
sider  that  the  EIFS  value  x2  shown  in  Fig.  A. 8  will  be  grown 
forward  using  each  of  the  three  approaches  indicated.  In 
Fig.  A. 9,  crack  growth  predictions  are  shown  for  each  of  the 
three  approaches  based  on  an  initial  flaw  size  of  x2<  The 
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Figure  A. 9.  Crack  Growth  Predictions  Must  Be  Consistent  with  the  Basis 
of  the  Equivalent  Initial  Flaw  Size . 


actual  fractographic  data  in  the  AL-AU  crack  size  range  used 
in  the  back-extrapolation  procedure,  is  denoted  by  a  solid 
circle  (•).  The  crack  growth  prediction  shown  in  Fig.  A. 9 
reflects  the  same  EIFS  master  curve  shown  in  Fig.  A. 8. 

The  following  observations  are  made  about  the  results 
depicted  in  Fig.  A. 9:  (1)  x2  was  defined  based  on  approach  2 
and  since  x2  was  grown  forward  the  same  way  the  fractographic 
results  were  back-extrapolated  the  predictions  are  consistent 
with  the  basis  of  the  EIFS  (note:  in  this  case  the  EIFS  mas¬ 
ter  curve  and  the  SCGMC  are  identical) ,  (2)  if  the  EIFS  x2 
is  grown  forward  using  the  SCGMCs  for  either  approach  ®  or  (§), 
the  crack  growth  predictions  are  not  consistent.  For  ex¬ 
ample,  the  crack  growth  predictions  for  approaches  ®  and  @ 
are  unconservative  and  conservative,  respectively,  compared 
to  the  observed  fractographic  results. 

A. 4. 3  Generic  EIFS  and  EIFSD  Considerations/Philosophy 

Ideally,  EIFSs  and/or  EIFSDs  should  be  "generic  quanti¬ 
ties"  that  depend  only  on  the  material,  actual  initial  flaws 
(e.g.f  scratches,  voids,  imperfections),  and  manufacturing 
processes.  Moreover,  such  initial  flaws  should  have  the  same 
geometry,  shape  and  details  as  actual  physical  flaws  that  may 
be  present  in  as-manufactured  fastener  holes.  Unfortunately, 
most  as-manufactured  initial  flaws  cannot  be  economically  de¬ 
tected  or  quantified  using  state-of-the-art  NDE  techniques. 
If  actual  initial  flaws  could  be  quantified,  other  problems 
would  still  remain:  (1)  how  to  mathematically  describe  and 
adequately  characterize  actual  initial  flaws  and  (2)  how  to 
model  the  crack  initiation  process  in  a  mechanistic  sense. 

The  crack  initiation  behavior  in  the  microstructure  is 
very  complex.  Such  factors  as  plasticity,  persistent  slip 
bands,  grain  boundaries,  crack  branching,  interactions  be¬ 
tween  microcracks,  etc.  are  important  considerations  for  de- 
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veloping  a  physically  meaningful  model  for  crack  initiation. 
Research  is  continuing  in  this  important  technical  area.  A 
suitable  mechanistic-based  crack  initiation  model  is  not  cur¬ 
rently  available  for  practical  durability  design  and  analysis 
applications.  In  lieu  of  such  a  model,  an  "engineering  ap¬ 
proach"  is  used  for  durability  analysis. 

Are  EIFSs  and/or  EIFSDs  generic?  Since  EIFSs  are  deter¬ 
mined  by  back-extrapolating  fractographic  results,  they  are 
not  strictly  "generic",  because  the  fractographic  results  are 
obtained  from  fatigue  cracks  generated  under  specified  test 
conditions  (e.g.,  load  spectrum,  stress  level,  %  bolt  load 
transfer,  fastener  hole  type/fit,  etc.).  Therefore,  the 
EIFSs  depend  on  the  test  conditions  reflected.  This  implies 
that  EIFSs  and/or  EIFSDs  have  to  be  developed  from  fracto¬ 
graphic  data  reflecting  the  applicable  conditions  not  only 
for  the  EIFSs  but  also  for  the  crack  growth  predictions. 
Since  there  are  many  possible  combinations  of  durability 
analysis  variables  (i,e.,  load  spectrum,  stress  level,  %  bolt 
load  transfer,  fastener  hole  type/fit) ,  this  could  result  in 
a  prohibitively  expensive  acquisition  program. 

The  real  issue  of  the  EIFS  approach  is  not  whether  the 
EIFSs  or  EIFSDs  are  "generic"  or  not  but  rather  can  reason¬ 
able  durabilty  analysis  predictions  be  obtained  for  a  given 
set  of  conditions  using  EIFSs  based  on  another  set  of  condi¬ 
tions.  A  data  pooling  procedure  has  been  developed  for  quan¬ 
tifying  EIFSD  parameters  for  one  or  more  fractographic  data 
sets.  Different  conditions  can  be  reflected  in  the  defini¬ 
tion  of  the  EIFSD.  This  provides  a  way  for  using  the  result¬ 
ing  EIFSD  for  variables  outside  the  range  of  the  baseline 
conditions.  The  derived  EIFSD,  based  on  the  data  pooling 
procedure,  is  justified  for  more  general  durability  analysis 
by  showing  that  predictions  for  the  distribution  of  TTCI  at  a 
givan  crack  size  or  the  distribution  of  crack  size  at  a  given 
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tin*  correlate  reasonably  wall  with  the  observed  fractograph- 
ic  data  baaa  uaad  to  define  the  EIFSD.  The  analyst  must  un¬ 
derstand  the  philosophy  and  justification  for  using  the  EIFS 
approach. 


A. 4. 4  Shape  and  Geometry  of  EIFSs 

EIFSs  are  determined  by  back-extrapolating  fractographic 
results  for  fastener  holes.  The  shape  and  geometry  as  well 
as  crack  initiation  sites  may  vary  depending  on  the  type  of 

fastener  hole/fit  as  well  as  the  amount  of  load  transfer 

through  the  fastener.  For  example,  it  has  been  observed  that 
the  dominant  fatigue  crack  in  an  unloaded  fastener  hole  typi¬ 
cally  originates  in  the  bore  of  a  hole  [e.g.,  16,  58]. 

whereas,  the  dominant  fatigue  crack  in  a  bolt  hole,  where  the 
bolt  is  loaded  in  single  shear,  usually  occurs  at  the  bolt 
hole  interface  with  the  mating  part  [e.g.,  16,  58].  For 

crack  growth  analysis,  the  fatigue  crack  in  the  bore  of  the 
hole  might  be  represented  by  a  through-the-thickness  initial 
flaw.  In  a  similar  manner,  the  fatigue  crack  at  the  bolt 
hole  interface  might  be  represented  by  an  initial  corner 
flaw.  The  shape  of  the  flaw  and  the  number  of  flaws  in  a 
fastener  hole  affects  the  stress  intensity  factor.  For  ex¬ 
ample,  the  stress  intensity  factor  for  a  single  corner  flaw 
is  different  than  that  for  two  corner  flaws  on  opposite  sites 
of  a  hole  or  for  a  single  through-the-thickness  flaw.  Deci¬ 
sions  about  flaw  shape,  size  and  locations  are  common  con¬ 
siderations  in  damage  tolerance  analysis.  Such  decisions  may 
vary  depending  on  the  analyst.  This  simply  points  out  that 
judgments  frequently  have  to  be  made  about  the  input  to  the 
LEFM  crack  analysis.  Similar  judgments  and  assumptions  have 
to  be  made  to  implement  the  probabilistic-based  durability 
analysis  methodology. 

The  "crack  size"  reflected  in  the  fractographic  data 
used  in  the  EIFS  determi nation  is  typically  measured  in  the 
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direction  of  crack  propagation.  In  thia  context,  the  shape 
of  the  crack  is  not  directly  considered  in  the  determination 
of  EIFSs  or  EIFSOs.  An  equivalent  initial  flaw  is  assumed  to 
be  a  "characteristic  crack  size"  which  applies  to  different 
flaw  shapes  in  a  crack  growth  analysis.  Using  the  data  pool¬ 
ing  procedure,  EIFSs  for  different  fractographic  results  and 
possibly  different  actual  flaw  shapes  are  treated  as  a  single 
population  of  EIFSs  for  purposes  of  determining  the  EIFSD. 

A. 4. 5  Initial  Flaw  sizes  for  Fuel  Leakage/ 

Ligament  Breakage  Requirements 

One  of  the  objectives  of  this  program  was  to  develop  a 
method  for  defining  an  initial  flaw  size  for  a  given  crack 
exceedance  probability.  In  the  beginning  we  hoped  that  a 
single  initial  flaw  size  could  be  used  to  make  crack  growth 
predictions  for  any  crack  size  associated  with  fuel  leakage/ 
ligament  breakage  requirements.  The  purpose  of  this  section 
is  to  discuss  the  determination  of  initial  flaw  size  and  un¬ 
derstandings  reached  under  this  program,  regarding  meaning 
and  use  of  an  EIFS. 

Three  key  issues  are:  (1)  how  to  determine  an  EIFS  for 
fuel  leakage/ ligament  breakage  applications  for  a  given  crack 
exceedance  probability  and  confidence  level,  (2)  what  limit¬ 
ations,  if  any,  are  placed  on  the  growth  of  the  EIFS  to  vari¬ 
ous  crack  sizes,  and  (3)  how  to  interpret  EIFSs  in  terms  of 
actual  in-service  cracks.  An  EIFS  for  a  given  probability  of 
crack  exceedance  can  be  defined  directly  from  the  EIFSD.  The 
size  of  the  equivalent  initial  flaw  depends  on:  (1)  the 
fractographic  crack  size  range,  AL-AU,  used  to  define  the 
EIFSD,  (2)  the  approach  used  to  back-extrapolate  the  fracto¬ 
graphic  results  (e.g.,  empirical  crack  growth  law  or  fracture 
mechanics-based  analytical  crack  growth  program) ,  (3)  crack 
growth  approach  (deterministic  on  stochastic),  (4)  goodness- 
of-fit  criterion,  (5)  statistical  distribution  function  used 
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to  define  the  EIFSD,  (6)  EIFS  upper  bound  limit  (xu) ,  and  (7) 
desired  confidence  level.  Tha  following  fundamental  princi- 
pla  must  ba  strictly  followad  whan  crack  growth  pradictiona 
ara  made  using  an  EIFS:  crack  growth  pradictiona  must  ba 
conaistant  with  tha  baaia  fcr  tha  EIFS  uaad.  For  axample,  if 
tha  EIFSO  ia  defined  for  a  selected  fractographic  crack  aize 
range  (e.g.,  AL-AU  «  0.01"  -  0.05")  and  given  empirical  crack 
growth  law  (e.g.,  da(t)/dt  »  Qa(t),  tha  EIFS  uaad  for  crack 
growth  predictions  is  also  limited  to  applications  using  the 
same  crack  size  range  and  empirical  crack  growth  law. 

Either  of  the  two  following  approaches  can  ba  used  to 
define  an  EIFS  from  the  EIFSO  and  to  grow  the  EIFS  forward  to 
crack  sizes  associated  with  fuel  leakage/ ligament  breakage 
(e.g.,  0.50"  -  0.75"):  (1)  use  a  single  EIFSD  and  two  or 
more  crack  growth  segments  to  make  crack  size  predictions  in 
the  desired  crack  size  range  and  (2)  define  the  EIFSD  based 
on  the  fractographic  results  for  the  crack  size  range  where 
crack  growth  predictions  are  described;  then  define  tha  ser¬ 
vice  crack  growth  master  curve  considering  the  same  AL-AU 
range  used  to  define  the  EIFSD.  The  two  approaches  above  are 
conceptually  described  and  compared  in  Fig.  A. 10. 

Any  initial  flaw  size  based  on  an  EIFD  must  be  consider¬ 
ed  in  the  same  context  as  the  EIFSD  was  determined.  This 
means  that  the  EIFS  should  be  interpreted  considering  every¬ 
thing  involved  in  its  definition  (e.g.,  fractographic  data 
used,  AL-AU  crack  size  range,  crack  growth  law  used  for  back- 
extrapolations,  goodness-of-fit  criterion,  etc.).  EIFSs  are 
not  physical  cracks  per  se  and  should  not  be  treated  as  such. 
Equivalent  initial  flaws  are  artificial  cracks  whose  charac¬ 
teristics  size  has  no  direct  relationship  to  actual  minimum 
crack  sizes  associated  with  fractographic  results  from  tear- 
down  inspections  or  other  fatigue  test  results.  This  per¬ 
spective  is  extremely  important  to  assure  that  initial  flaws 
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Fiaure  A. 10.  Comparison  of  Two  Deterministic  Crack  Growth  Approaches^for  Making 
Crack  Growth  Predictions  Based  on  EIFSP  and  a  Single  EIFS. 


based  on  an  EIFSD  are  proparly  usad  in  any  durability  analy¬ 
sis. 

A. 5  PROBABILITY  OF  CRACK  EXCEEDANCE/ 

CUMULATIVE  DISTRIBUTION  OF  TTCI 

Tha  two  most  important  predictions  in  any  durability  an¬ 
alysis  are:  (1)  probability  of  crack  axcaadanca  at  any  time 
t,  p(i,  t  ) ,  and  (2)  cumulative  distribution  of  TTCI  for  any 
rafaranca  crack  size  x^,  FT(t) .  Assume  tha  EIFSD  has  bean 
defined  and  that  a  deterministic  crack  growth  approach  is 
usad  to  grow  tha  EIFSD  forward  using  a  service  crack  growth 
master  curve  (SCGMC)  which  reflects  selected  design  vari¬ 
ables.  For  example,  in  Fig.  A. 11  predictions  for  p(i,T)  and 
FT(t)  are  illustrated  for  three  different  stress  levels, 
where  o  >  ^  >  a  All  SCGMCs  must  be  consistent  with  the 
basis  for  the  EIFSD.  In  Fig.  A. 11 (a),  the  probability  of  ex¬ 
ceeding  crack  size  x^  x2  and  x3  is  equal  to  the  cross-hatch¬ 
ed  area  under  the  crack  size  distributions  shown.  Therefore, 
for  equal  p(i,  t  )  >  x2  >  x3> 

Predictions  for  the  cumulative  distribution  of  TTCI, 
FT ( t ) ,  are  illustrated  in  Fig.  A. 11(b).  In  this  case,  the 
cross-hatched  area  under  the  TTCID  represents  (FT(t) .  Also 
note  in  Fig.  A. 11(b)  that  for  the  same  FT(t)  value,  T^  <  T2 
<  T3 .  This  trend  would  be  expected  since  >  o2  >  °3’ 

SCGMCs  can  be  determined  using  either  suitable  fracto- 
graphic  results  or  a  LEFM  analytical  crack  growth  program 
[e.g.,  38,  39].  Details  and  guidelines  for  determining 
SCGMCs  which  are  compatible  with  the  basis  for  the  EIFSD  are 
given  in  Sections  II  and  III. 
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APPENDIX  B 


STATISTICAL  SCALING  TECHNIQUE 

The  EIFS  distribution,  Fa(0)(x),  for  faotanar  holas  is 
dafinad  for  a  "single  hole  population.”  Tast  specimens  for 
acquiring  fatigua  crack  growth  data  may  hava  one  or  more 
fastanar  holas  par  specimen.  Soma  spacimans  may  not  be  fati¬ 
gue  tested  to  failure.  Also,  every  fastanar  hola  in  each  re¬ 
plicate  tast  specimen  may  not  contain  a  measurable  fatigue 
crack  or  else  tha  crack  is  too  small  or  complex  (o.g.,  multi¬ 
ple  crack  origins  and  branching)  for  fractographic  analysis. 
Hence,  a  statistical  scaling  technique  has  bean  developed  for 
determining  the  equivalent  initial  flaw  size  distribution  for 
a  single  hole  population  using  the  largest  fatigue  crack  per 
specimen  [2,  13].  This  minimize  the  fractographic  reading 
requirements,  permits  a  maximum  utilization  of  the  available 
fractographic  data  and  allows  for  "mixing  and  matching”  of 
fractographic  data  for  specimens  with  a  different  number  of 
holes. 

A  technique  is  described  in  this  section  for  statisti¬ 
cally  scaling  the  TTCID  or  EIFSD.  This  general  technique, 
conceptually  described  in  Fig.  B.l  applies  to  any  distribu¬ 
tion  function  used  to  represent  the  TTCID  or  EIFSD.  Various 
aspects  are  discussed. 

Although  every  fastener  hole  in  replicate  test  specimen 
may  be  equally  stressed,  there's  no  guarantee  that  every  hole 
will  have  a  measurable  fatigue  crack  when  the  specimen  fails 
or  when  the  fatigue  test  is  stopped  before  failure.  In  any 
case,  the  cracking  resistance  of  each  fastener  hole  in  each 
IFQ  test  specimen  should  be  accounted  for  when  determining 
the  EIFSD. 
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Statistical  Scaling  Concept  for  "Single  Hole  Population 


There's  another  important  aapact  about  scaling.  With 
scaling,  ths  EIFSD  can  bs  determined  using  only  ths  largsst 
fatigus  crack  par  spacimsn.  Tsst  apscimsns  for  diffsrsnt  re- 
plicate  data  sets  may  contain  one  or  more  fastener  holes  per 
specimen.  For  example,  the  test  specimens  in  one  replicate 
data  set  may  have  only  one  fastener  hole  per  specimen  and 
those  in  another  data  set  may  contain  four  fastener  holes  per 
specimen.  Scaling  accounts  for  the  number  of  fastener  holes 
in  each  specimen  and  is  the  common  denominator  for  obtaining 
homogeneous  results  for  different  fractographic  data  sets. 

A  "single  hole  population"  accounts  for  the  actual  num¬ 
ber  of  fastener  holes/specimen  in  the  fractographic  data  base 
used  to  define  either  the  TTCID  or  the  EIFSD.  There  are  two 
ways  to  do  this:  (1)  use  the  fractographic  data  for  the 
largest  fatigue  crack  in  each  fastener  hole/ specimen,  or  (3) 
use  only  the  largest  fatigue  crack  in  1  of  JL  holes  per  speci¬ 
men  and  a  scaling  technique  to  account  for  1  of  JL  holes  pur 
specimen  (Fig.  B.l). 


Let  the  cumulative  distribution  of  EIFS  for  a  single 
hole  population  be  denoted  by  Fa(0) (x)  and  that  of  the  maxi¬ 
mum  EIFS,  based  on  the  largest  fatigue  crack  per  specimen 
with  fastener  holes,  be  denoted  by  Fft  ^ (x) .  Assuming 
that  fatigue  cracking  in  each  equally  stressed  fastener  hole 
of  a  specimen  is  statistically  independent  of  that  of  the 
other  holes,  Fa  (0) (x)  is  related  to  Ffl  ^ (x)  as  shown  in 
Eq.  B-l, 


(B-l) 


where  £  =  number  of  equally  stressed  fastener  holes  per 
specimen. 
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Ltt  the  cumulative  distribution  of  crack  size  at  any 

a(t) 


time  T  for  a  "single  hole  population"  be  denoted  by  F^^fx) 


and  that  for  a  "1  of  l  hole  population"  be  denoted  by 
F^(t)  ^  Fa(t)  and  Fa  (t)  ^  can  b®  determined  from 

Fa(0)  ^  and  Fa  (0)^x^'  respectively,  using  a  service  crack 
growth  master  ^curve,  (SCGMC)  denoted  by  y^T).  Similarly, 
let  FT(t)  and  F^  (t)  denote  the  cumulative  distribution  of 
TTCI  for  a  "single  hole  population"  and  for  a  "1  of  JL  hole 
population",  respectively.  The  notations  and  meanings  are 
illustrated  in  Fig.  B.2. 


The  cumulative  distribution  of  crack  size  for  any  time, 
7;  can  be  determined  from  the  EIFSD  using  the  applicable 
SCGMC.  When  the  initial  flaw  size,  y^JT),  is  grown  forward 
using  the  SCGMC,  a  crack  size  x1  is  obtained  at  time,  T  , 
(Fig.  B.3). 


The  cumulative  distribution  of  crack  size  at  any  time, 
7t  is  obtained  from  the  respective  EIFSDs  as  follows.  In  Egs. 
B-2  and  B-3,  y^fT)  is  the  initial  flaw  size  corresponding  to 
a  crack  size  x1  at  time,  T ,  as  determined  from  the  SCGMC.  An 
expression  for  y^CT)  is  given  in  Eq.  B-4. 


*  f«o. 


(8-2) 


(8-3) 


Yu(t)  -  x{exp{>-  QtT)  ' 

The  cumulative  distribution  of  crack  size  is  assumed  to 
be  statistically  compatible  with  the  cumulative  distribution 
of  TTCI.  Equations  B-5  through  B-8  reflect  this  assumption. 


0  For  "1  of  X  Holt  Population" 

®  For  "single  hole  population"  (or  1  of  £  hola  population) 


(a)  Notations  for  Cumulative  Distribution  of  Crack  Size 


(b)  Notations  for  Cumulative  Distribution  of  TTCX 


Figure  B.2,  Notations  and  Concepts  for  Cumulative  Distribution 
of  Crack  Size  and  TTCI.  “  ~ 
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If  the  largest  fatigue  crack  in  each  fastener  hole  Is 
statistically  independent  of  that  of  the  other  holes,  FT(t) 
and  FT^(t)  ar*  related  as  follows  in  Egs.  B-9  and  B-10. 
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in 
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(B-10) 


Equations  B-l,  B-9  and  B-10  are  the  key  relationships 
for  scaling.  These  equations  are  useful  when  the  data  pool¬ 
ing  procedure  is  used  to  estimate  the  EIFSD  parameters  for  an 
equivalent  single  hole  population. 


In  some  cases,  fractographic  results  may  not  be  avail¬ 
able  for  each  fastener  hole  per  specimen  in  a  replicate  data 
set.  Whenever  scaling  is  used  to  define  the  EIFSD,  it  must 
also  be  accounted  for  when  the  EIFSD  is  grown  forward  to 
check  the  goodness-of-fit  of  the  cumulative  distribution  of 
crack  size  and/or  TTCI. 
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APPENDIX  C 


DERIVATION  OF  EQUATIONS  FOR  DETERMINING  Q^  POOLED  Q 
AND  STANDARD  DEVIATION 

C.l  EQUATIONS  FOR  Q^  AND  POOLED  Q 

Equations  for  determining  Q1  and  pooled  Q  are  derived  in 
this  section.  Q^  and  Q  are  empirical-based  crack  growth 
parameters  in  a  crack  size-time  relationship.  Q^  describes 
the  crack  growth  for  a  single  crack  in  a  data  set,  and  Q  de¬ 
scribes  the  typical  crack  growth  for  a  complete  fractographic 
data  set.  Since  Q  is  used  to  characterize  the  typical  crack 
growth  for  all  cracks  in  a  data  set,  it  is  referred  to  as  a 
"pooled  Q"  value. 

A  suitable  crack  size-time  relationship  is  needed  to 
back-extrapolate  fractographic  results  to  time  zero  and/or  to 
grow  on  EIFSD  forward  to  make  durability  analysis  predic¬ 
tions.  A  useful  crack  size-time  relationship,  given  in  Eq. 
C-l,  has  been  derived  from  the  crack  growth  law  given  in  Eq. 
C-l  [13],  In  eqs.  C-l  and  C-2,  a(t)  ■  crack  size  at  any  time 
t,  a(0)  -  crack  size  at  t  ■  0  or  EIFS,  Q  *  crack  growth  para¬ 
meter  (empirical  constant)  and  da(t)/dt  -  crack  growth  rate. 
The  relationship  given  in  Eqs.  C-l  and  C-2  have  been  shown  to 
very  reasonable  for  durability  analyses  [2,13,15,24-30). 

a(t)  »  a(0)exp  (Qt)  (C-l) 

daUVdt  *  Qa(t)  (C-2) 

For  a  single  crack  application  the  notations  in  Eqs.  C-l 
and  C-2  are  modified  as  shown  in  Eqs.  C-3  and  C-4. 

ait)  *  a(0) exp  (Qft)  ( C-3 ) 


C-l 


daitVdt  m  Q  ait) 


(C-4) 


Qj  in  Eq.  C-3  can  ba  determined  using  a  iaast  square  fit 
procedure  and  fractographic  data  or  analytical  predictions 
based  on  a  LEFM  crack  growth  program  [e.g.,  38]. 

Equation  C-3  can  ba  transformed  into  a  linear  least 
square  fit  form  by  taking  the  natural  log  of  both  sides  of 
Eq.  C-3  as  shown  in  Eq.  C-5.  A  least  squares  fit  expression 


In  ait)  =  ^naiO)  +•  Q  t 


(C-5) 
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X, 


for  Qi  is  given  in  Eq.  C-6,  where  N  -  number  of  X,  Y  pairs 
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considered  in  the  least  squares  fit.  X^  and  Y^  are  defined 
in  Eq.  C-5. 


The  expression  for  Q  or  "pooled  Q"  is  derived  as 

2 

follows.  An  expression  for  sum  squared  error  (SEE),  E  , 
based  on  Eq.  C-2,  is  given  in  Eq.  C-7.  In  Eq.  C-7, 


.v 

SSE  =  E2  *  V 

i-i 

B  «  lnQ,  and  N  ■  number  of  cracks  in  the  data  set.  After 
minimizing  E2  with  respect  to  B  (i.e.,  csE2/dB  «  0)  and 
simplifying  Eq.  C-8  is  obtained.  The  resulting  expression 
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for  Q  is  given  in  Eg.  09.  Fractographic  data  oan  ba  uaad  to 
dafino  tha  da(t)/dt  and  e(t)  valuaa  in  Eq.  C-9. 


Q  -  exp 


(C-9) 


An  expression  for  Q  can  ba  darivad  from  Eq.  C-9  in  terms 
of  Q  as  follows,  Assume:  (1)  tha  fractographic  data  to  ba 
used  to  determine  covers  tha  salactad  crack  size  range 
AL-AU  and  (2)  each  crack  in  tha  data  sat  is  homogeneous.  Tha 
fatigue  cracks  are  made  "homogeneous H  by  artificially  equally 
spacing  the  a(t)s  for  each  crack  in  tha  data  sat.  Thus,  each 
fatigue  crack  in  tha  data  set  has  an  equal  number  of 
eqr illy-spaced  a(t) 's  in  tha  AL-AU  range.  We  will  show  that 
it  doesn't  matter  how  many  equally  spaced  cracks  are  used. 


Tha  tarms  of  the  numerator  of  Eq.  C-9  can  ba  generalized 
to  sum  results  for  individual  fatigue  cracks  in  the  data  sat 
as  follows.  In  Eqs.  C-10  through  C-12,  ai(t  )  «  crack  size 
for  crack  i  at  time  t  and  -  crack  growth  parameter  for  a 
single  crack  as  defined  by  Eq.  C~6.  The  desired  expression 
for  Q,  given  in  Eq.  C-13,  is  obtained  by  substituting  Eq. 
C-12  into  Eq.  C-9. 
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Equation  C-13  ia  very  convenient  for  computing  Q  lor  a 
given  fractographic  data  set.  Zn  other  words,  Q  is  equal  to 
the  exponential  of  the  average  lnQi  for  the  data  set. 

Equation  C-13  was  derived  assuming  that  each  fatigue 
crack  in  the  data  set  had  an  equal  number  of  equally-spaced 
a(t)s  per  crack.  This  is  equivalent  to  assuming  that  n  \,q 
fixed  for  the  data  set.  Equation  C-13  is  also  valid  when  in 
varies  for  each  fatigue  crack  in  the  data  set.  For  examp?,e\ 
in  Eqs.  C-10  and,  C-11,  substitute  n^  for  N  in  this  case; 
where  -  the  number  of  a(t)s  for  the  ith  crack.  With  this 
substitution,  the  last  two  terms  in  Eq.  C-11  still  cancel  out 
leaving  the  term  shown  in  Eq.  C-12. 

C-2  EQUATION  FOR  STANDARD  DEVIATION 

When  the  stochastic  crack  growth  approach  (SCGA)  is  used 

(see  Section  III)  the  standard  deviation,  is  required. 

z 

Expressions  for  o*z  are  derived  and  discussed  in  this  sec¬ 
tion. 


Equation  C-2  is  transformed  by  taking  the  natural  log  of 
each  side  of  the  equation  and  the  result  is  given  in  Eq. 
C-14. 


In  daltVdt  =  (nQ  +•  (n  alt)  (C-14 ) 

Suppose  in  da(t)/dt  versus  In  a(t)  is  plotted  based  on 
the  fractographic  results  for  a  given  data  set.  Equation 
C-14  can  be  used  to  determine  the  relationship  between 
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d»(t)/dt  and  a(t)  and  is  a  "bast  fit"  straight  lina  as  illus¬ 
trated  in  Fig.  C.l.  By  definition,  is  equal  to  the  sum 
squared  error  (SSE)  between  the  straight  line  given  in  Eq. 
C-14  divided  by  the  number  of  da(t)/dt,  a(t)  pairs  in  the 
data  set.  The  SSE  in  this  case  is  given  in  Eq.  C-15.  In  Eq. 
C-15 ,  N  -  total  number  of  da(t)/dt,  a(t)  pairs  in  the  data 
set.  An  expression  for  oi /  obtained  using  Eq.  C-15,  is  given 
in  Eq.  C-16 .  In  Eq.  C-16  the  summation  is  from  i  ■  1  to  N, 
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where  N  is  the  total  number  of  a(t)  points  in  the  data  set  or 
in  the  selected  crack  size  range,  AL-AU. 


Equation  C-16  can  be  simplified  to  obtain  an  expression 

for  crt  in  terms  of  and  Q.  To  simplify,  Eq.  C-4  is 

substituted  into  Eq.  C-16  and  a  summation  notation  is  used  to 

obtain  the  expressions  for  or,  as  shown  in  Eq.  C-17. 
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In  Eq.  C-17,  M  -  number  of  cracks  in  the  data  set. 
for  each  crack  can  be  obtained  using  Eq.  C-6,  which  accounts 
for  the  number  of  a(t),  t  pairs  in  the  selected  AL-AU  range 
for  a  given  crack. 
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Figure  C.l  Conceptual  Illustration  for  Standard 
Deviation  ~~~~ 
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APPENDIX  D 


DURABILITY  ANALYSIS  EQUATIONS  BASED  ON  SELECTED 
EIFSD  FUNCTIONS  AND  DCOA 

D.l  INTRODUCTION 

Expression  for  ths  cumulative  distribution  of  crack  size 
at  anytime  t,  Fa^)(x),  and  the  cumulative  distribution  of 
TTCI ,  FT(t)  are  presented  in  this  section  for  these  different 
EIFSD  functions:  (1)  two-parameter  Weibull,  (2)  lognormal 

and  (3)  lognormal  compatible  [12]. 

Expressions  for  Fa^tj (x)  and  FT(t)  are  obtained  using  a 
transformation  of  the  applicable  EIFSD  function,  Fa(o)^x^' 
and  the  deterministic  crack  size-time  ultimately  given  in  Eg. 
D~l.  Equation  D-l 

EIFS  s=  x  -  x^expi-Qt)  (D-l) 

is  based  on  the  crack  growth  rate  model  of  Eq.  C-2  and  the 
generalized  crack  size-time  relationship  of  Eq.  C-l  in  which 
a(t)  *  and  a(0)  -  EIFS  -  x.  In  Eq.  D-l,  t  ■  time  or  TTCI 
corresponding  to  crack  size  and  Q  “empirical  crack  growth 
constant.  The  cumulative  distribution  of  EIFS  is  separately 
represented  by  a  two-parameter  Weibull,  a  lognormal  and  a  log¬ 
normal  compatible  distribution  function.  Details  are  pro¬ 
vided  in  the  following. 

D.2  TWO-PARAMETER  WEIBULL  EIFSD  FUNCTION 

The  two-parameter  Weibull  distribution  can  be  used  to 
represent  the  cumulative  distribution  of  EIFS,  denoted  by 
F  (x),  as  given  in  Eq.  D-2 .  In  Eq.  D-2,  x  -  EIFS  random 
variable,  a  and  (3  are  empirical  constants. 
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The  cumulative  distribution  of  crack  size  at  any  tins  t 
can  bo  obtainad  through  a  transformation  of  Eg.  D-2  and  ths 
crack  size-time  ralationship  givan  in  Eq.  D-l  as  follows*  In 
Eq.  D-3 ,  tha  tarms  ara  dafinad  as  givan  in  Eqs.  D-l  and  D-2. 
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An  axprassion  for  tha  cumulativa  distribution  of  TTCI 
can  be  obtainad  from  Eq.  D-3  as  follows. 

FrU)  *  1  -  fatl(*>  (D-4 ) 
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The  two  parameters  Weibull  function  for  the  cumulative 
distribution  of  EIFS,  given  in  Eq.  D-2,  contains  two  "EIFSD 
parameters"  to  be  determined  (i.e.,  oc  and  fi  ).  These  same 
parameters  are  also  included  in  Eqs.  D-3  and  D-4.  Therefore, 
the  EIFSD  parameters  can  be  determined  using  either  Eq.  D-2, 
D-3  or  D-4.  Since  the  two-parameter  Weibull  EIFSD  function 
does  not  have  an  EIFS  upper  bound  limit  (xu) ,  infinite 
initial  flaw  sizes  are  possible. 

D. 3  LOGNORMAL  EIFSD  FUNCTION 


The  cumulative  distribution  of  EIFS  based  on  the  lognor¬ 
mal  distribution  function  is  given  in  Eq.  D-5,  where: - 
standard  normal  distribution  function,  x  -  EIFS  random  vari¬ 
able;  M  and  r  are  empirical  constants. 
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An  expression  for  the  cumulative  distribution  of  crack 
sice  at  any  time  t  can  ba  obtainad  through  a  transformation 
of  Eq.  D-S  using  Eq.  D-l.  Tha  tarms  in  Eq.  D-6  ara  dafinad 
in  Eqs.  D-l  and  D-S. 
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Tha  cumulativa  distribution  of  TTCI  can  ba  obtainad 
using  Eq.  D-7  as  follows. 
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The  lognormal  EIFSD  function  in  Eq.  D-5  contains  two 
EIFSD  parameters  to  be  determined  (i.e.,  m.  and  <r  ).  Since 
these  same  parameters  are  also  included  in  Eqs.  D-6  and  D-7, 
either  Eq.  D-5,  D-6,  or  D-7  can  be  used  to  estimate  the  EIFSD 
parameters.  Since  there  is  no  EIFS  upper  bound  limit,  xu, 
for  the  lognormal  EIFSD  function,  infinite  EIFS  values  are 
possible. 


D. 4  LOGNORMAL  COMPATIBLE  EIFSD  FUNCTION 

The  lognormal  compatible  EIFSD  is  given  in  Eq.  D-8,  in 
which  xu  -  EIFS  upper  bound  limit,  m  and  <r  are  empirical 

constants,  <p  *  standard  normal  distribution  function  and 
x  =  EIFS  random  variable.  Equation  D-8  was  obtained  by 
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transforming  ths  lognormal  cumulative'  distribution  of  TTC1 
using  ths  crack  size-time  rslationship  given  in  Eg.  D-l  [12]. 

An  sxprsssion  for  ths  cumulative  distribution  of  crack 
size  at  any  time  t,  denoted  by  Fa(t)C>0'  can  be  obtained 
through  a  transformation  of  Eq.  D-d  using  Eq.  D-l.  The 
resulting  expression  is  given  in  Eq.  D-9,  where:  xx  « 
reference  crack  size  for  TTCI,  t  -  TTCI;  xu,  a*  and  cr  are 
defined  in  Eq.  D-8. 
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An  expression  for  the  cumulative  distribution  of  TTCI, 
denoted  by  F_(t),  can  be  obtained  as  follows.  All  terms  in 
Eq.  D-10  have  already  been  defined. 
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The  lognormal  compatible  function  for  the  cumulative 
distribution  of  EIFS,  Eq.  D-8,  contains  three  "EIFSD  para¬ 
meters"  to  be  determined  (i.e.,  xu,  m  and  <r  ).  The  same 
EIFSD  parameters  are  also  included  in  Eqs.  D-9  and  D-10  for 
the  cumulative  distribution  of  crack  size  and  of  TTCI,  re¬ 
spectively.  Therefore,  the  EIFSD  parameters  can  determined 
using  either  Eq.  D-8,  D-9  or  D-10. 


APPENDIX  E 


METHODS  FOR  DETERMINING  POOLED  Q  AND  EIFSa  FOR 

DCGA  AND  SCGA 


E.l  INTRODUCTION 

The  purpose  of  this  section  is  to  review:  (1)  different 
methods  and  variation  for  determining  the  "pooled  Q"  value 
for  the  crack  growth  rate  model  shown  in  Eq.  E-l,  (2)  appro¬ 
aches  for  determining  DCGA-based  and  (3)  approaches  for 
determining  SCGA-based  EIFSs.  In  Eq.  E-l,  da(t)/dt  -  crack 

datiVdta  Qa(t)  (E-l) 

growth  rate,  Q  -  empirical-based  crack  growth  rate  constant 
and  a(t)  -  crack  size  at  time  t.  Methods  are  reviewed  and 
discussed  for  durability  analysis  applications. 

E.2  APPROACHES  FOR  DETERMINING  POOLED  Q 

Two  different  approaches  were  investigated  for  estimat¬ 
ing  the  "pooled  Q"  value  for  a  given  fractographic  data  set: 
(1)  direct  use  of  the  crack  growth  rate  model  given  in  Eq. 
E-l  and  (2)  use  of  the  generalized  crack  size-time  rela¬ 
tionship  given  in  Eq.  C-l.  The  two  approaches  are  conceptu¬ 
ally  compared  in  Fig.  E.l.  The  Q  parameter  in  Eqs.  E-l  and 
C-l  is  an  empirical  constant  usually  based  on  fractographic 
results.  Q  can  be  estimated  using  a  least  squares  fit  pro¬ 
cedure  using  either  Eq.  E-l  or  C-l.  Both  equations  are  con¬ 
sidered  separately  in  the  following. 


E-l 


E.2.1  Q  Based  on  da(t)/dt 


Q  can  bs  datsminsd  from  Eq.  E-l  using  the  transforma¬ 
tion  shown  in  Eq.  E-2,  the  least  squares  fit  procedure  and 
fractographic  data.  The  crack  growth  rate,  da(t)/dt,  for  any 
a(t) 


(n  daitVdt  »  (n  Q  +  ( n  ait) 


(E-2) 


can  be  estimated  from  the  fractographic  data  (i.e.,  a(t)  ver¬ 
sus  t  results)  using  either  the  modified  secant  method  [10] 
or  the  five-point  incremental  polynomial  method  [12].  Q  in 
Eq.  E-2  can  be  determined  using  Eq.  E-3  which  is  based  on  the 
least  squares  fit  procedure.  The  following  subscripted 
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subscripted  notation  is  used  in  Eq.  E-3:  M  -  number  of  fa¬ 
tigue  cracks  in  the  data  set;  •  number  of  a(t) 's  for  the 
ith  fatigue  crack;  subscripts  i  and  j  denote  the  ith  crack 
and  jth  value;  a^(t^)  «  jth  value  of  crack  size  for  the  ith 
fatigue  crack;  and  [da(t.)/dt]  ^  -  jth  crack  growth  rate  for 
the  ith  fatigue  crack.  Equation  E-3  can  be  interpreted  as 
follows.  The  "pooled  Q"  value  for  a  given  data  set  depends 
on  the  crack  growth  rates  for  each  a(t)  for  each  fatigue 
crack  in  the  fractographic  data  set. 

To  be  consistent,  Q  should  be  determined  using  the  frac¬ 
tographic  data  in  a  selected  fractographic  crack  size  range, 
AL-AU.  The  da(t)/dt  data  in  the  AL-AU  range  can  be  used  to 
compute  Q  various  ways:  (1)  use  only  the  actual  da(t)/dt, 
a(t)  data  that  lies  in  the  AL-AU  range  for  each  fatigue  crack 


in  the  data  sat,  and  (2)  artifically  aqualiza  tha  number  of 
da(t)/dt,  a(t) 's  in  tha  AL-AU  range.  With  tha  first  approach 
there  may  not  be  an  equal  number  of  a(t)'s  for  each  fatigue 
crack  in  tha  data  set  in  the  AL-AU  range.  This  means  that 
each  fatigue  crack  may  not  be  equally  weighted  in  the  deter¬ 
mination  of  Q.  One  way  to  equally  weight  each  fatigue  crack 
is  to  equalize  the  number  of  a(t)'s  in  the  selected  AL-AU 
range.  All  fatigue  cracks  in  the  data  set  may  not  cover  the 
same  crack  size  range.  Also,  the  smallest  and  largest  a(t) 
in  the  AL-AU  range  will  probably  be  different  for  each  fa¬ 
tigue  crack.  The  influence  of  the  crack  size  range  and  the 
fatigue  crack  weighting  on  the  resulting  Q  has  been  investi¬ 
gated  (see  Volume  II  [  3]). 

E.2.2  Q  Eased  on  Crack  Size-Time  Relationship 

Q  can  also  be  determined  for  a  given  fractographic  data 
set  using  the  crack  size-time  relationship  given  in  Eq.  C-l. 
The  main  advantage  of  this  approach  is  that  "Q"  can  be  deter¬ 
mined  directly  from  the  fractographic  results  (i.e.,  a(t) 
versus  t)  without  computing  da(t)/dt's.  Equation  C-l  can  be 
transformed  into  a  linear  least  squares  fit  form  as  shown  in 
Eq.  40.  Then,  Q  can  be  determined  using  Eq.  C-14,  based  on 
the  least  squares  fit  procedure.  In  Eq.  E-4,  the  subscripts 
i  and  j 
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denote  the  ith  fatigue  crack  and  the  jth  a(t),t  values  for 
the  ith  fatigue  crack  M  «  total  number  of  fatigue  cracks  in 
the  fractographic  data  set,  N^  -  total  number  of  a(t)'s  for 
the  ith  fatigue  crack;  and  are  defined  in  Eq.  E-5. 
In  Eq.  E-5,  t^  »  jth  time  for  the  ith  fatigue  crack  and 
«i(tj)  -  jth  crack  size  for  the  ith  fatigue  crack. 

y'JLu>  <e*5> 

V  1  J 

Using  Eq.  E-4,  the  "pooled  Q"  value  for  a  given  fracto¬ 
graphic  data  set  can  be  determined  using  the  a(t),t  values 
for  each  fatigue  crack  in  the  data  set.  This  approach  is 
recommended  for  durability  analyses.  Whatever  approach  is 
used  to  compute  Q,  it  is  important  to  recognize  that  "Q" 
depends  on  the  fractographic  crack  size  range  used. 

E . 3  APPROACHES  FOR  DETERMINING  DCGA-BASED  EIFSs 

The  following  three  different  approaches  for  determining 
EIFSs  based  on  the  DCGA  were  investigated:  (1)  averaging 
method,  (2)  least  squares  fit  method,  and  (3)  use  "pooled  Q" 
and  selected  x^t.  A  conceptual  comparison  of  the  three  ap¬ 
proaches  is  shown  in  Fig.  E.2. 

The  "pooled  Q"  value  for  a  given  fractographic  data  set 
is  required  to  implement  each  of  the  three  approaches.  Meth¬ 
ods  for  computing  "pooled  Q"  are  described  in  Section  E.2. 
EIFSs  for  each  of  the  three  approaches  should  be  determined 
using  the  same  fractographic  crack  size  range,  AL-AU,  used  to 
determine  the  "pooled  Q"  value.  Each  of  the  three  approaches 
for  computing  DCGA-based  EIFSs  are  described  and  discussed  in 
the  following. 
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£.3.1  Averaging  Mathod 


Tha  "avaraging  mathod”  ia  basad  on  tha  a(0)  ralationahip 
givan  in  Eg.  E-6. 
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This  aquation  is  obtainad  from  Eq.  C~1  by  substituting  a(t)  ■ 
x1  and  solving  for  a(0).  Tha  crack  size-time  ralationship 
given  in  Eqc  E-6  is  used  to  compute  tha  average  EIFS,  a(0), 
for  each  fatigue  crack  in  tha  data  sat  as  follows.  An  EXFS 
is  computed  for  each  a(t),t  pair  in  the  selected  AL-AU  range 
for  a  given  fatigue  crack.  Then,  an  average  EXFS  is  comput¬ 
ed.  This  approach  is  mathematically  described  in  Eq.  E-7. 


AuerageEIFS  =  a  (0)  = 
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in  Eq.  E-7,  h  -  total  numb.r  of  „(t),t  pairs  in  tha  a.lactad 

AL-AU  crack  size  rang.,  aft*)  -  1th  crack  size,  t,  -  ith 

time,  and  Q  .  "pooled  Q«  value  for  given  data  set  (see  sec- 
tion  E.2)  . 

E.3.2  Least  Squares  Fit  Method 


The  EIFS  for  a  given  fatigue  crack  is  determined  as 
follows.  The  crack  size-time  relationship  given  in  Eq.  E-6 
is  transformed  into  a  linear  least  squares  fit  form  as  shown 
in  Eq.  C-14,  An  expression  for  a(0)  or  EIFS,  given  in  Eq. 

E-8,  is  obtained  using  Eq.  c-14  and  a  least  squares  fit  pro¬ 
cedure.  in  Eq.  E-8,  N  «  total  number  of  a(t),t  pairs  in  the 


a(0)  =  EIFS  =  exp 


i  *  i 


i  ■  l 


X 


(E-8) 
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selected  AL-AL  prack  sin#  rang*  and  Q  »  "pooled  Q  valus"  for 
a  givsn  fractographic  data  sst.  X ,  and  Y,  ars  dsfined  in  Eq. 

i  A 

E-9,  where: 

Xi  -  ti  (E-9) 

Y*>  In  a(ti) 

-  ith  tins  corresponding  to  the  crack  size  a(t^)  in  the 
selected  AL-AU  crack  size  range. 

E.3.3  Selected  x^,t 

\ 

This  approach  is  a  simple  variation  of  the  averaging 
method  described  in  Section  E.3.1.  For  example,  the  E1FS  for 
a  given  fatigue  crack  can  be  determined  directly  from  Eq.  E-6 
for  selected  x^  and  t  values  (e.g. ,  see  Fig.  E.2(c)).  This 
approach  is  convenient  to  use  when  the  EIFSD  parameters  are 
to  be  estimated  using  a  "TTCI  fit."  In  this  case,  the  TTCIs 
for  a  given  data  set  are  defined  for  a  selected  reference 
crack  size,  Then,  the  EIFSs  for  each  fatigue  crack  in 

the  data  set  can  be  determined  using  Eq.  E-6. 

E.3.4  Discussion 

The  three  methods  for  determining  EIFSs  based  on  the 
DCGA  are  evaluated  in  Volume  II  [  3 ] .  No  significant  differ¬ 
ences  in  comparable  EIFS  values  were  found  for  any  of  the 
three  methods  described  in  this  section.  Whatever  method  is 
used  to  determine  the  EIFSs,  the  resulting  EIFSs  values  will 
depend  on  the  crack  growth  rate  model  used  as  well  as  the 
fractographic  crack  size  range  used.  Therefore,  a  different 
set  of  EIFSs  will  be  obtained  for  each  crack  growth  rate 
model  and  AL-AU  crack  size  range  used. 
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Different  crack  growth  rata  modal*  could  ba  usad  to 
determine  tha  EIFSs.  Tha  almpla  modal  dascrlbad  in  Eq.  E-l 
has  baan  found  to  ba  vary  raasonabla  for  describing  tha  crack 
size-time  ralationship  for  diffarant  materials,  load  spactra, 
stress  levels,  fastener  hole  type  (i.a.,  straight  bore  or 
countersunk),  %  bolt  load  transfer  anu  stress  level  [2,3,10,15, 
24 '30). 

E.4  APPROACHES  FOR  DETERMINING  SCGA-BASED  EIFSs 

Procedures  for  determining  SCGA-based  EIFSs  are  de¬ 
scribed  and  discussed  in  this  section.  Deterministic  and 
stochastic-based  EIFSs  can  be  determined  from  tha  same  frac- 
tographic  data  set  as  illustrated  in  Fig.  E*3.  For  tha  DCGA, 
the  EIFSs  are  determined  using  tha  same  EIFS  master  curve. 
On  the  other  hand,  SCGA-based  EIFSs  are  based  on  a  separate 
EIFS  curve  for  each  fatigue  crack  in  the  fractographic  data 
set.  DCGA  and  SCGA  based  EIFSs  are  discussed  in  Section  A. 3 
and  the  two  approaches  are  conceptually  compared  in  Fig.  A. 2. 

SCGA-based  EIFSs  are  used  in  the  durability  analysis  ex¬ 
tension  described  in  Section  III.  With  the  SCGA,  each  fa¬ 
tigue  crack  in  the  fractographic  data  set  is  back-extrapolat¬ 
ed  to  time  zero  using  the  fractographic  results  for  each 
crack  (separately)  and  a  reasonable  crack  growth  rate  model, 
such  as  that  described  in  Section  III. 

The  crack  size-time  relationship  shown  in  Eg.  C-3  can  be 
used  to  compute  the  SCGA-based  EIFS  as  follows.  Equation  C-3 
can  be  interpreted  as  follows  for  the  SCGA:  a(t)  ■  crack 
size  at  any  time  t,  a(0)  -  SCGA-based  EIFS,  and  ■  empiri¬ 
cal  crack  growth  constant  for  the  ith  fatigue  crack.  Equa¬ 
tion  C-3  is  transformed  into  a  linear  least  squares  fit  form 
as  shown  in  Eq.  C-5  so  that  a(0)  and  can  be  determined. 
The  empirical  constant  can  be  estimated  using  Eq.  C-6. 
Then,  a(0)  can  be  determined  using  either  Eq.  C-3  or  Eq.  E-6 
(with  Q  -  Q^)  as  conceptually  described  in  Fig.  A. 2.  Like 


the  DCGA,  the  SCGA-based  EIFSs  are  determined  using  a 
salactad  fractographic  crack  siaa  ranga  (e.g.,  AL-AU  -  o.oi" 
-  0.05").  To  ba  consistent,  tha  valua  for  aach  fatigua 
crack  must  raflact  tha  same  AL-AU  ranga. 

Tha  "avaraging  method"  described  in  Section  E.3.1  can 
also  be  used  to  compute  SCGA-based  EIFSs.  For  example,  in 
Eq.  E-7  simply  use  the  crack  growth  rate  parameter  for  a 
given  crack,  Q^,  instead  of  the  "pooled  Q"  value  for  the 
complete  fractographic  data  set. 

It  is  shown  in  Volume  II  [  3 ]  that  it  doesn't  make  a  lot 
of  difference  whether  Eq.  C-3,  E-6  or  E-7  is  used  to  deter¬ 
mine  the  EIFSs.  The  EIFS  values  for  a  given  fatigue  crack 
may  vary  depending  which  equation  is  used  to  determine  the 
EIFS.  However,  the  net  effect  on  the  EIFSD  parameters  and 
respective  initial  flaw  sizes  for  a  given  probability  is  not 
significant. 

Once  the  SCGA-based  EIFSs  have  been  determined,  the 
EIFSD  parameters  can  be  estimated  using  EIFS  data  and  the 
CLSSA  described  in  Appendix  F. 
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APPROACHES  FOR  ESTIMATING  EIFSD  PARAMETERS 

F.l  INTRODUCTION 

Five  approaches  ars  described  and  discussed  in  this  sec¬ 
tion  for  estimating  the  EIFSD  parameters  for  selected 
Fm,n, (x)  functions  including:  (1)  combined  least  square  sums 
approach  (CLSSA) ,  (2)  method  of  moments  (MM) ,  (3)  homogeneous 
EIFS  approach  (HEIFS) ,  (4)  upper  tail  fit  and  (5)  non-linear 
approach.  The  five  different  approaches  are  discussed  in 
terms  of  the  Weibull  compatible  EIFSD  function.  The  CLSSA 
is  also  described  for  the  lognormal  compatible  EIFSD 
function. 


F.2  COMBINED  LEAST  SQUARE  SUMS  APPROACH  (CLSSA) 

The  CLSSA  is  described  and  discussed  in  this  section  for 
two  different  EIFSD  functions:  (1)  Weibull  compatible  and 
(2)  lognormal  compatible.  This  approach  can  be  used  to  esti¬ 
mate  the  EIFSD  parameters,  in  a  global  sense,  for  one  or  more 
fractographic  data  sets.  Since  the  effects  of  scaling  are 
accounted  for,  fractographic  data  sets  can  be  "mixed  and 
matched"  to  estimate  the  EIFSD  parameters.  A  step-by-step 
procedure  for  implementing  the  approach  is  presented. 

The  EIFSD  parameters  for  a  given  EIFSD  function  can  be 
estimated  using  either  TTCI  or  EIFS  input  data  for  one  or 
more  fractographic  data  sets.  The  TTCIs  or  EIFSs  (whichever 
is  used)  for  each  fractographic  data  set  are  ranked  and  least 
square  sums  are  obtained  for  each  data  set  using  the  appli¬ 
cable  EIFSD  function  and  data  used  to  fit  the  EIFSD  paramet¬ 
ers.  Each  least  square  sum  accounts  for  the  scaling  factor 
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for  each  data  sat.  Tha  combined  laaat  aquara  suns  for  all 
data  sata  ara  obtainad  by  summing  tha  applloabla  valua  for 
aach  data  aat.  Than,  tha  EIFSD  parameters  ara  estimated 
uaing  tha  combinad  laaat  aquara  auma  and  a  global  fit. 

F.2.1  Ganaral  Procadura 

Tha  procadura  for  implamanting  tha  CLSSA  is  givan  balow: 

1.  Salact  fractographic  data  seta  to  ba  used  in  the 
pooling  procedure  to  estimate  the  EIFSD  parameters  in  a 
global  sense. 

2.  Acquire  the  TTCIs  or  EIFSs  for  each  data  set  to  be 
used  to  estimate  the  EIFSD  parameters. 

3.  Define  scaling  factor,  ,  for  each  data  set  to  be 

used. 

4.  Select  the  EIFSD  function  to  be  used  and  then  trans¬ 
form  it  into  a  linear  least  square  fit  form.  If  EIFSs  are 
used  for  the  fit,  use  the  expression  for  the  cumulative  dis¬ 
tribution  of  EIFS  for  the  fit.  On  the  other  hand,  if  TTCIs 
are  used  for  the  fit  use  the  expression  for  the  cumulative 
distribution  of  TTCI  based  on  the  transformation  of  the  EIFSD 
and  the  crack  growth  model  used  (e.g.,  Eq.  C-2). 

5.  Assume  a  value  for  the  EIFS  upper  bound  within  the 
constraint:  largest  EIFS  in  any  data  set  <  xu  <  0.05”.  com¬ 
pute  the  combined  least  square  sums  for  all  data  sets  com¬ 
bined.  Use  these  results  to  estimate  the  EIFSD  parameters. 
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6 .  Tha  SIFSD  parameters  can  he  optimizad  within  the  xu 
constraint.  A  rsasonahla  way  to  do  this  is  to  minimize  the 
total  standard  arror  for  all  data  sats  used  to  determine  tha 
EIFSD  parameters.  Tha  optimum  xu  can  ba  datarminad  by  trial 
and  arror  until  tha  total  standard  arror  has  baan  minimi zad. 


F.2.2  Waibull  Compatible  Function 


Tha  purpose  of  this  section  is  to  describe  how  tha  EIFSD 
parameters  can  be  estimated  for  the  Weibull  compatible  EIFSD 
function  using  the  CLSSA.  Equations  are  presented  and  dis¬ 
cussed  for  estimating  the  EIFSD  parameters  based  on  "TTCI 
fit"  and  on  "EIFS  fit."  The  effects  of  scaling  are  accounted 
for. 

F.2.2.1  EIFS  Fit.  Tha  Weibull  compatible  EIFSD  function 
Fa(0)(x)»  is  given  in  Eq.  F-l,  where  xu  -  upper  bound  EIFS 
limit, 
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x  -  EIFS  variable,  and  «  and  <t>  are  empirical  constants. 
Equation  F-l  is  modified  to  account  for  scaling  as  shown  in 


Eq.  F-2 .  In  Eq.  F-2, 


a- (0) 


(x)  -  cumulative  distribution  of 


maximum  EIFS  based  on  the  largest  fatigue  crack  per  specimen 
in  1  of  £  fastener  holes.  Other  terms  in  Eq.  F-2  have  already 
been  defined.  Since  Fft  ^ (x)  ■  (x)  when  there  is  no 

scaling,  Eq.  F-2 
cases  where  t  >  l. 


Fa*(0)  (x) 

is  'a  general  expression  that  applies  to 
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The  EIFSD  parameters  in  Eq.  F-2  (i.e.,  xu,  <*and$)  can 
be  estimated  using  the  CLSSA  as  follows.  Transform  Eq.  F-2 
into  the  linear  least  square  fit  form  given  in  Eq.  F-3.  For 
data  pooling  purposes  the  subscripted  notation  shown  in  Eq. 
F-4  is  used. 
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In  Eq.  F-4 ,  -  number  of  fastener  holes  for  the  repli¬ 

cate  specimen  in  the  ith  data  set,  *  jth  EIFS  values  for 
the  ith  data  set,  xu  -  EIFS  upper  bound  limit,  *  number  of 

EIFS  values  for  the  ith  data  set,  F_  ,rtv(x44)  -  cumulative 

a jt.  (0)  ij7 

distribution  of  maximum  EIFS  based  on  the  largest  fatigue 
crack  per  specimen  for  the  ith  data  set. 


The  EIFSD  parameters  and  ^  in  Eq.  F-3  can  be  estimated 
for  a  given  xu  value  using  a  least  squares  fit  procedure. 
Expressions  for  estimating  «  and  <p  using  the  CLSSA  are  shown 
in  Eq.  F-5  and  F-6,  respectively. 


V 

M  V  i 


V  V 

,Vf  •  t  Vf  t 


iN.IIVrlia  1  «r 


Ota 


»  -1  i«l  ;«l 


t » j  i«i _ t •<  ; ■  t 


.v 


.V 


V  M  "t  M  •  t 

2>, 1 1*;-  i  l*„ 

i«l  i •  1 


(F-5) 


<t > -  exp 


«  \  *  *. 

«I  I  I  I  Y„ 

t-1  ;■ l  i*l  ;«l 

M 

i-i 


(F-6) 


F-4 


unuunt  »  a.  kJLiOLtf  wnim  v*.  v*  ux  im.  vx  uvuxvx  umvx  w  w  ux  L'v wuv  lv  w  intuvs*  lv  l  v yv  _v  lv  l  v  l  v  lv  lv  \jv  •-•m 


Zn  Eqs.  P-5  and  F-6,  H  »  n unbar  of  EZF8  data  sets  uaad 
to  estimate  the  EZFSD  parameters,  K.t  -  numbar  of  EZFSs  in  tha 
ith  data  aat.  Equations  F-5  and  F-6  can  also  ba  used  whan  a 
TTCI  fit  is  usad.  If  an  EIFS  fit  is  usad,  than  and 
ara  dafinad  as  shown  in  Eq.  F-4. 

EIFS  Fit  Procadura 


1.  Salact  EIFS  data  sats  to  ba  usad  to  estinata  tha 
EIFSD  paramatars  oc  and 

2.  Rank  EIFSs  in  ascanding  order  for  aach  data  aat  sep¬ 
arately.  Use  Fa  (0)  (x)  -  j/(Mi  +  1). 

k 

3.  Assume  value  of  xu  within  constraint:  largest  EIFS 
in  any  data  set  <  xu  <  0.05'*.  Assume  tha  maximum  xu  corre¬ 
sponds  to  tha  initial  flaw  size  usad  for  damage  tolerance 
requirements. 

4.  compute  the  least  square  sums  shown  in  Eqs.  F-5  and 
F-6  for  each  EIFS  data  set  separately,  including  the  appro¬ 
priate  scaling  factor,  j^,  for  each  data  set. 

5.  Combine  the  least  square  sums  for  each  EIFS  data  set 
into  a  total  sum  for  all  data  sets  to  be  used  in  the  pooling 
procedure. 

6.  Use  the  iterative  procedure  described  in  Fig.  9  to 
optimize  the  EIFSD  parameters.  The  total  standard  error, 
TSE,  for  all  EIFS  data  sets  can  be  determined  using  Eq.  F-7. 
The  subscripted  notation  and  terms  in  Eq.  F-7  are  defined  the 
same  way  as  given  in  Eqs.  F-4  through  F-6. 
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F.2.2.2  TTCI  Fit.  The  cumulative  distribution  of  ttci, 
FT(t),  is  given  in  Eq.  F-8. 
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In  Eq.  F-8,  xu  «  EIFS  upper  bound  limit,  x,  *  reference 
crack  size  for  TTCIs,  t  *  TTCI,  Q  =  empirical-based  crack 
growth  rate  parameter,  and  a  and  <p  are  empirical  constants. 


The  EIFSD  parameters  for  the  Weibull  compatible  EIF-SD 
function  can  be  determined  using  a  "TTCI  fit"  as  follows. 
Firsts  the  expression  for  the  cumulative  distribution  of 
TTCI;  given  in  Eq.  F-8  is  modified  to  account  for  scaling  as 
shown  in  Eq,  F-9.  In  Eq.  F-9,  FT  (t)  ■  cumulative  distribu¬ 
tion  of  minimum  TTCI  baaed  on  the  largest  fatigue  crack  in  1 
of  i  holes  per  specimen.  Other  terms  in  Eq.  F-9 


V1-'-  V(x' 
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have  already  been  defined.  Since  FT^(t)  "  FT(t)  when  SL  «  1 
(no  scaling),  Eq.  F-9  applies  for  >  1. 

Equation  F-9  can  be  transformed  into  a  linear  least 
square  fit  form  as  shown  in  Eq,  F-XO. 
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To  implement  the  data  pooling  procadura,  tha  terms  in 
Eq.  F-10  ara  axpraasad  in  a  subscripted  form  as  shown  in  Eq. 
F-ll.  in  Eq.  F-ll,  »  numbar  of  fastanar  holes  for  the 
replicata  specimen  in  tha  ith  data  sat,  m  number  of  TTCI 
valuas  for  tha  ith  data  sat,  xu  -  EIFS  uppar  bound  limit,  xx 
-  aQ  -  rafaranca  crack  siza  for  TTCIs,  -  poolad  Q  value 
for  tha  ith  data  sat,  tij  -  jth  TTCI  value  for  tha  ith  data 
rTj.  ^i j ^  Vftlua  for  tha  cumulative  distribution  of 

TTCI  minimum?  (i.a.,  based  on  tha  largest  fatigue  crack  per 
specimen  in  e*ich  data  set)  for  the  ith  data  sat. 
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Tha  tflFSD  parameters  *  and  <P  in  Eq.  F-9 
mined  for  a  selected  xu  value  using  Eqs.  F-5, 


can  be  detar- 
F-6,  and  F-ll. 


TTCI  Fit  Procedure 


1.  Select  TTCI  data  sets  to  be  used  to  estimate  the 
EIFSD  parameters  <*and  ^  . 


2 .  Rank  the  TTCIs  in  ascending  order  for  each  data  set 

separately.  Use  Fr  it)  =  j/(N.  +  1)  • 

3.  Assume  value  of  xu  with  constraint:  largest  EIFS  in 
any  data  set  <  xu  <  0.05". 

4.  Compute  the  least  square  sums  in  Eq.  F-5  and  F-6  for 
each  TTCI  data  sat  separately,  including  the  applicable  scal¬ 
ing  factor,  JZj,  ,  for  each  data  sat.  Use  the  terms  defined  in 
Eq.  F-ll. 


5.  Combine  tha  least  square  sums  for  each  TTCI  data  set 
into  a  total  sum  for  all  data  sets  to  be  used  in  tha  pooling 
procedure. 
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6.  Us*  ths  iteration  procedure  described  in  Fig.  9  to 
optimize  the  EIFSD  parameters.  The  total  standard  error, 
TSE,  for  all  TTCI  data  sets  can  be  determined  using  Eq.  F-12. 
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In  Eq.  F-i:2,  H  -  number  of  fractographic  data  sets  used  in 
the  pooling  procedure. 

F.2.3  Lognormal  Compatible  EIFSD  Function 

The  purpose  of  this  section  is  to  describe  how  the  EIFSD 
parameters  can  be  estimated  for  the  lognormal  compatible 
function  using  the  CLSSA.  Equations  are  presented  and  dis¬ 
cussed  for  estimating  the  EIFSD  parameters  based  on  an  "EIFS 
fit"  or  "TTCI  fit."  The  effects  of  scaling  are  accounted 
for. 

F.2.3.1  EIFS  Fit.  The  cumulative  distribution  of  EIFS  for 
the  lognormal  capatible  distribution  is  given  in  Eq.  F-13, 
where 
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$  «  standard  normal  distribution  function,  xu  -  FIFS  upper 
bound  limit,  x  -  EIFS  variable,  Q  -  empirical  crack  growth 
rat*  constant,  and  m.  and  o'  are  empirical  constants. 

Equation  F-13  is  modified  to  account  for  scaling  as 
shown  in  Eq.  F-14.  Equation  F-14  is  transformed  into  a 
linear  least  squares  fit  form  as  shown  in  Eq.  F-15. 
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For  data  pooling  purposes,  the  terms  in  Eq.  F-15  are  repre¬ 
sented  in  a  subscripted  form  as  shown  in  Eq.  F-16.  In  Eq. 
F-16 ,  the  subscripts  i  and  j  refer  to  the  ith  EIFS  data  set 
and  the  jth  EIFS  value  for  the  ith  data  set.  Qi  *  "pooled  Q" 
value  for  the  ith  data  set,  inverse  of  the  standard  nor¬ 
mal  distribution  function,  and  *  number  of  EIFSs  .r,  the 
ith  data  set. 
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The  EIFSD  parameters  m  and  < r  in  Eq.  F-14  and  F-15  can 
be  estimated  for  a  given  xu  using  the  CLSSA.  Expressions  for** 
and  <r  are  given  in  Eqs.  F-17  and  F-18,  respectively.  The 
subscripted  notations  in  Eqs.  F-17  and  F-18  are  the  same  as 
those  in  Eq3.  F-5  and  F-6  for  a  global  fit. 
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S IFS  Fit  Procedure 


The  •tap-by-step  procedure  for  estimating  the  EIFSD 
parameters  using  an  "EIFS  fit'*  is  the  same  as  that  given  in 
section  F.2.2.1  except  as  follows: 

1.  In  step  4  compute  the  least  square  sums  for  Eqs. 
F-17  and  F-18  using  Eq.  F-16. 

2.  In  step  6  determine  the  total  standard  error ,  TSE, 
using  Eq.  F-19. 
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F.2.3.2  TTCI  Fit .  The  cumulative  distribution  of  TTCI , 
FT(t),  is  shown  in  Eq.  F-20.  To  account  for  "scaling"  Eq. 
F-20  is  modified  as  shown  in  Eq.  F-22.  For  data  pooling 
purposes,  a  subscripted  notation  is  used  as  shown  in  Eq. 
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in  Eq.  F-23,  the  subscripts  i  and  j  refer  to  the  ith 
TTCI  data  sat  and  the  j  the  TTCI  value  for  the  ith  data  set, 
*  "pooled  Q"  value  for  the  ith  fractographic  data  set,  t.. 
■  jth  TTCI  value  for  the  ith  data  set,£^  -  scaling  factor  for 
the  ith  data  set,  O'*  -  inverse  of  the  standard  normal  distri¬ 
bution  function  and  “  number  of  TTCIs  in  the  ith  data 
set. 


The  EIFSD  parameters  u  and  o'  in  Eq.  F-21  and  F-22  can 
be  estimated  for  a  given  x„  using  the  CLSSA.  Using  the  ex¬ 
pressions  of  Eq.  F-23,  and  can  be  determined  using  Eqs. 
F-17  and  F-18 ,  respectively. 

TTCI  Fit  Procedure 


The  step-by-step  procedure  for  estimating  the  EIFSD 
parameters  using  a  "TTCI  fit"  is  the  same  as  that  given  in 
section  F.2.2.2  except  as  follows: 

1.  In  step  4  compute  the  least  square  sums  for  Eqs. 
F-17  and  F-18  using  Eq.  F-19. 

2.  In  step  6  determine  the  total  standard  error.  (TSE) 
using  Eq.  F-24. 
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F. 3  METHOD  OF  MOMENTS  APPROACH 


The  method  of  moments  can  be  used  to  estimate  the  con** 
stants  o <  and  4  in  the  two-parameter  Weibull  cumulative  dis¬ 
tribution  function  shown  in  Eq.  F-25  [e.g.,  61 3  *  In  E9* 
F-25,  x  is  a  random  variable;  *  and  4  are  empirical  con¬ 
stants.  The  method  of  moments  approach  can  also  be  used  to 
estimate  the  EIFSD  parameters  for  the  Weibull  compatible 
EIFSD  function  using  either  a  "EIFS  fit"  or  a  "TTCI  fit." 
The  method  of  moments  approach  accounts  for  scaling  of 
different  fractographic  data  sets  and  can  be  used  to  estimate 
the  EIFSD  parameters  for  one  or  more  data  sets  in  a  global 
sense. 


F(x)  -  1  -  exp 


(F-25) 


The  usual  way  the  constants  <*•  and  A  in  Eq.  F-25  are 
determined  by  the  method  of  moments  is  as  follows: 

1.  compute  the  mean  of  the  random  variable,  x,  using 
Eq.  F-26,  where  N  «  number  of  x^  values  in  the  sample. 
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2. 


Compute  the  standard  deviation  of  x  using  Eq. 


F-27 . 

(F-27) 


3.  Compute  the  coefficient  of  variation  of  x  using  Eq. 


F-28 . 
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4.  at  can  be  determined  from  the  expression  for  CCV(x) , 
given  in  Eq.  F-29,  where  ro«  standard  gamma  function. 
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A  table  of  COV(S)  versus  <x  values  can  be  obtained  from  Eq. 
F-29  using  assumed  ot  values.  ot  can  then  be  determined  for 
the  given  COV(x)  by  interpolating  the  results  in  the  table  of 
COV(x)  versus  ot  values. 

5.  The  constant  in  Eq.  F-25  is  determined  from  Eq. 

F-30. 
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The  above  general  procedure  can  be  used  to  estimate  the 
EIFSD  parameters  for  the  Weibull  compatible  EIFSD  function  as 
follows.  Expressions  for  the  probabiity  of  crack  exceedance, 
p(i,T  ),  and  the  cumulative  distribution  of  TTCI,  FT  (t) ,  are 
given  in  Eqs.  F-31  and  F-32,  respectively.  In  Eqs.  f-31  and 
F-32 ,  Ji  -  scaling  factor,  xy  »  EIFS  upper  bound  limit,  ot  and£ 
are  empirical  constants.  In  Eq.  F-31,  x  -  EIFS  variable  and 
in  Eq.  F-32,  x1  -  reference  crack  size  for  TTCI,t.  Equations 
F-31  and  F-32  both  account  for  scaling.  In  these  equations 
note  that:  F  }  -  Fft(0)  -  Pa(0)(x)  and  PT  (t)  -  PT(t)  when 
tm  1.0  (no  scaling). 
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Equation  F-31  can  ba  rawrittan  in  an  aquivalant  form  as 
Eq.  F-25  as  shown  in  Eq.  F-33.  Equation  F-33  can  ba  usad 
whan  EZFSs  ara  usad  to  astimata  tha  EZFSD  parameters  oc  and  ^ 
for  a  givan  xu.  In  Eq.  F-33,  sinca  tha  term  je^inc  xu/x)  is 

dimensionless  the  mean  of  this  term  can  ba  determined  for  a 
givan  xu  using  all  tha  EIFS  data  sets  combined.  An  iterative 
procedure  and  method  of  moments  approach  ara  described  later 
for  estimating  <x  and  P  for  a  givan  xu> 


p  (ijt)  =  1  -  exp 


lx) 

U 


a 


<p 


(F-33) 


Equation  F-32  for  the  cumulative  distribution  of  TTCI 
can  also  be  reduced  to  the  same  form  as  Eq.  F-25  by  imposing 
the  condition  that  x1  «  xu.  This  condition  requires  that  the 
reference  crack  size  for  TTCIs  be  equal  to  the  EIFS  upper 
bound  limit  (xu)  used.  The  resulting  expression  for  F_  (t) 
is  given  in  Eq.  F-34.  In  Eq.  F-34,  since  the  term  i"6t  is 
dimensionless/  the  TTCIs  for  one  or  more  data  sets  can  be 
used  to  estimate  oc  and  P  for  a  given  xu  using  the  method  of 
moments . 
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Tha  method  of  momenta  approach  for  estimating  tha  EIFSD 
parameters  for  tha  Waibull  compatible  EIFSD  function  la 
schematically  illuatratad  in  Fig.  F.l.  Varioua  aapacta  and 
dataila  of  thia  approach  will  now  ba  diacuaaad  for  both  EIFS 
and  TTCI  fita.  Eithar  an  EIFS  or  a  TTCI  fit  can  ba  uaad  to 
estimate  tha  EIFSD  parameters  for  tha  Waibull  compatible 
EIFSD  function. 

Let  the  terms  in  the  numerator  of  Eqa.  F-33  and  F-34  be 
defined  as  shown  in  Eqs.  F-35  and  F-36,  respectively. 


WL 

A  =  ((  ^)( nix  lx  ) 

tj  I  U  J 


(F-35) 


V. 


(F-36) 


A  subscripted  notation  is  used  for  fractographic  data  pooling 
and  tracking  purposes.  The  subscripts  i  and  j  denote  tha  ith 
fractographic  data  set  and  jth  value  from  the  ith  data  set. 
In  Eqs.  F-35  and  F-36,  -  scaling  factor  for  ith  data  set, 

trial  value  of  oi  in  iterative  process,  Xj  -  jth  value  of 
EIFS  in  the  ith  data  set,  ■  empirical  crack  growth  con¬ 
stant  ("pooled  Q"  value)  for  the  ith  data  set  and  tj  -  jth 
TTCI  value  for  a  given  reference  crack  size  x1  for  the  ith 
data  set. 

Expressions  for  the  mean,  standard  deviation  and  coef¬ 
ficient  of  variation  (COV)  for  7  are  given  in  Eqs.  F-37, 
F-38  and  F-39 ,  respectively.  These  equations  are  written  in 
a  subscripted  form  for  data  pooling  purposes  and  computer 
applications.  The  same  general  expressions  for  the  mean, 
standard  deviation  and  COV  for  7  can  be  used  for  either  an 
EIFS  or  TTCI  fit.  For  example,  if  a  EIFS  fit  is  used  let 

*ij  "  \j  and  a  TTCI  used  In  E,3®* 
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F-37  -  F-39,  the  subscripts  i  and  j  denote  the  ith  data  set 
and  jth  value  from  the  ith  data  set,  M  ■  number  of  data  sets 
used  in  the  pooling  procedure  and  -  number  of  samples  for 
the  ith  data  set. 
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Since  one  of  the  constants  to  be  determined  (i.e.,  oc  ) 
is  included  in  Eqs.  F-35  and  F-36,  an  iterative  technique  is 
used.  The  general  procedure  is  schematically  described  in 
Fig.  F.l.  Computer  software  has  been  developed  to  implement 
the  method  of  moments  approach  described  in  this  section 
[6]. 

F . 4  HOMOGENEOUS  EIFS  APPROACH 

With  the  homogeneous  EIFS  (HEIFS)  approach  EIFSs  are 
computed  for  each  fractographic  data  set  and  the  resulting 
EIFSs  for  all  data  sets  are  pooled  and  ranked  into  one  homo- 
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geneous  EIFS  population.  The  EIFSD  parameters  can  then  be 
estimated  using  the  applicable  EIFS  population  and  a  least 
squares  fit  procedure. 

The  HEIFS  approach  presented  in  this  section  has  two 
major  limitations:  (l)  it  doesn't  account  for  scaling  (i.e., 
lm  1})  of  different  fractographic  data  sets  used;  and  (2) 
since  the  EIFSs  for  all  data  sets  are  combined  into  one  homo¬ 
geneous  population,  the  inherent  scatter  for  each  data  set  is 
obscured  in  the  process.  Scaling  of  different  fractographic 
data  sets  can  be  accounted  for  using  the  HEIFS  approach  and 
an  iterative  procedure.  However,  details  are  not  presented 
herein.  If  the  HEIFS  approach  is  used  to  estimate  the  EIFSD 
parameters  for  multiple  fractographic  data  sets^  scaling 
should  be  accounted  for. 

The  HEIFS  approach  will  be  discussed  assuming  the  EIFSD 
is  represented  by  the  Weibull  compatible  function.  The 
general  procedure  is  conceptually  described  in  Fig.  F.2  and 
in  the  following. 

1.  Compute  the  EIFSs  using  the  DCGA  for  each  fracto¬ 
graphic  data  set  to  be  used  to  estimate  the  EIFSD  parameters 
in  a  global  sense. 

2.  Pool  and  rank  the  EIFSs  for  all  data  sets  into  one 
homogeneous  population.  Let  -  number  of  EIFSs  in  the  ith 
EIFS  data  set  and  M  ■  total  number  of  EIFS  data  sets  to  be 
used.  Then,  the  total  number  of  EIFSs  of  all  data  sets  com¬ 
bined  is  denoted  by  n  as  shown  in  Eq.  F-40. 
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Figure  F.2. 


Conceptual  Illustration  of  the  Homogeneous  EIFS 
Approach  for  Defining  Initial  Fatigue'TSality 


The  cumulative  distribution  of  EIFS,  (x) ,  cein  b«  deter¬ 

mined  using  Eg.  F-41. 
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3.  The  Weibull  compatible  EIFSO  function,  Fa(0)^' 
given  in  Eg.  F-42,  can  be  transformed  into  a  linear  least 
sguare  fit  form  as  shown  in  Eg.  F-43.  In  Eg.  A-7,  x  «  EIFS 
random  variable;  xu  »  EIFS  upper  bound  limit;  ot  and  are 
empirical  constants. 
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4.  The  terms  in  Eg.  F-43  are  modified  into  a  sub¬ 
scripted  notation  to  facilitate  data  pooling  and  implementa¬ 
tion  using  ti  computer.  Expressions  are  shown  in  Eg.  F-44. 
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5.  The  empirical  constants  oi  and  ^  can  be  estimated  for 
a  given  xu  using  the  least  squares  fit  procedure  as  follows. 
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6.  The  EIFSD  parameters  oc  and  ♦  can  be  optimized  for  a 
given  x  by  minimizing  the  total  standard  error,  TSE,  within 
the  constraints  of  xu  (i.e.,  largest  EIFS  in  the  population  < 
xu  1  0.05").  An  iterative  procedure  can  be  used  to  minimize 
the  TSE. 


F. 5  UPPER  TAIL  FIT 

For  functional  impairment  due  to  fuel  leakage  and  liga- 
m  nt  breakage,  the  main  interest  is  in  the  larger  EIFSs  in 
the  upper  tail  of  the  EIFSD.  A  method  was  investigated  for 
estimating  the  EIFSD  parameters  using  only  the  largest  EIFSs. 
The  interest  here  was  to  determine  the  EIFSD  based  on  an 
"upper  tail  fit"  and  then  to  determine  the  initial  flaw  size 
corresponding  to  a  selected  probability  of  crack  exceedance. 
The  actual  study  is  documented  in  Volume  II  [3]  and  the 
technical  details  are  given  in  the  following. 

Basic  assumptions  of  the  upper  tail  fit  approach  are: 
(1)  the  EIFSD  parameters  can  be  estimated  using  the  largest 
EIFSs  in  the  EIFS  population;  and  (2)  the  initial  flaw  size 
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for  a  given  probability  can  be  obtained  from  the  derived 
EIFSD  and  the  resulting  initial  flaw  size  should  be  more 
realistic  than  if  the  entire  EIFS  population  is  used  to 
estimate  the  EIFSD  parameters. 

Procedure 


1.  Compute  EIFSs  for  each  EIFS  data  set  to  be 

considered  in  the  global  fit.  » 

2.  Once  the  EIFSs  have  been  determined,  the  CLSSA  or 
the  HEIFS  approach  (if  JL  *  1  applies)  can  be  used  to  estimate 
the  EIFSD  parameters.  If  the  CLSSA  is  used,  the  EIFSs  for 
each  data  set  are  ranked  in  ascending  order.  A  probability 
of  crack  exceedance  is  selected  (e.g.,  p(i,T)  -  0.30)  and 
ranked  EIFSs  <  0.30  are  used  to  estimate  the  EIFSD  parameters 
using  the  CLSSA.  If  the  HEIFS  approach  is  used,  all  the 
EIFSs  are  pooled  and  ranked  into  one  homogeneous  population. 

A  probability  of  crack  exceedance,  p(i,7*)  is  selected  and  the 
ranked  EIFSs  in  the  upper  tail  of  the  EIFSD  are  used  to 
estimate  the  EIFSD  parameters.  In  this  case  a  least  squares 
fit  procedure  is  used. 

3.  The  EIFSD  parameters  can  be  optimized  for  a  given 

EIFS  upper  bound  limit,  xu,  by  minimizing  the  total  standard 
error  for  the  fit.  This  can  be  accomplished  using  a  trail 
and  error  procedure  as  follows:  (1)  assume  xQ  within  con¬ 
straint  (i.e.,  largest  EIFS  in  population  <  xu  <  0.05");  (2) 

select  p(i,T);  (3)  compute  0(,  and  total  standard  error, 

TSE;  and  (4)  repeat  (1)  to  (3)  until  the  TSE  is  a  minimum. 

4 

F , 6  NON-LINEAR  APPROACH 

i 

A  non-linear  approach  for  estimating  the  EIFSD  para¬ 
meters  for  the  Weibull  compatible  EIFSD  function  was  investi¬ 
gated.  To  estimate  the  EIFSD  parameters  for  the  Weibull 
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compatible  function  using  tha  least  squares  fit  procedure , 
the  appropriate  equation  must  first  be  transformed  into  a 
linear  least  squares  fit  form.  A  In  In  transformation  is  re¬ 
quired.  Since  the  least  squares  error  is  measured  in  terms 
of  the  transformed  equation,  rather  than  the  original  equa¬ 
tion,  there  is  a  tendency  to  give  more  weight  to  the  extremal 
input  data  than  the  intermediate  data  when  the  linear  least 
squares  fit  method  is  used.  For  this  reason,  a  non-linear 
approach  was  investigated.  Results  are  presented  in  Volume 
II  [3].  The  approach  is  described  in  this  section  and  the 
equations  are  derived  for  implementing  the  method. 

A  general  non-linear  approach  for  estimating  model  con¬ 
stants  has  been  developed  [62].  This  approach  can  be  used  to 
estimate  the  EIFSD  parameters  for  the  Weibull  compatible 
EIFSD  function  in  a  global  sense  (non-linear  fit) .  Details 
are  as  follows. 

The  cumulative  distribution  of  TTCI  based  on  the  Weibull 
compatible  EIFSD  function  and  the  deterministic  crack  growth 
rate  model  is  given  in  Eq.  F-9. 

The  "EIFSD  parameters"  are  xu,  OC  and  <p  •  In  Eq.  F-9, 
note  that:  F^ ( t )  ■  FT(t)  when  1  (no  scaling). 

The  residual  equation  for  Eq.  F-9  is  given  in  Eq.  F-47, 
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»  predicted  cumulative  distribution 
of  TTCI  based  on  a  selected  xu  and 
initial  values  of  and  ^ . 
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ot  and  <p>  are  the  EIFSD  parameters;  and  <pQ  are  initial 
values  for  OL  and  <p  ,  respectively;  Ao/q  and  Afy  are  increments 
of  U  and  to  be  determined.  The  partial  derivative, 

(3f/doc)0  and  (3f/d^)0#  are  given  in  Egs.  F-48  and  P-49, 
respectively,  where  /  is  defined  in  Eq.  F-50. 
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In  Eq.  F-47 ,  Aot^  and  A^q  can  b®  determined  using  the 
linear  least  squares  fit  procedure.  The  sum  squared  error, 
E2,  for  Eq.  F-47  is  given  in  Eq.  F-51. 
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Expression*  for  Ao^  and  A4>0  can  be  obtained  by  minimizing  E2 
(i.e.,  cJE'&x  -0  and  £E2£^  -  0  and  then  solving  two  linear 
simultaneous  equations.  The  resulting  expressions  for 
and  Ao^  are  given  in  Eqs.  F-52  and  F-53,  respectively.  A 
subscripted  notation  is  used  to  facilitate  data  pooling. 
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The  ,  Z^j,  and  terms  in  Eqs.  F-52  and  F-53  are 

defined  in  Eqs.  F-54,  F-55,  F-56,  and  F-57,  respectively. 
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In  Eqs.  F-54  through  F-57,  tha  subscripts  i  and  j  denote  tha 
ith  TTCI  data  aat  and  tha  jth  valua  from  tha  ith  data  sat, 

^ ,  -  scaling  factor  for  tha  ith  TTCI  data  aat,  ■  poolad  Q 
valua  for  tha  ith  fractographic  data  aat,  tj  *  jth  TTCI  valua 
in  tha  ith  data  aat;  othar  terms  hava  alraady  baan  defined. 

A  schematic  of  the  non-linaar  procedure  for  estimating/ 
optimizing  the  EIFSD  parameters  for  one  or  more  TTCI  data 
sets  is  shown  in  Fig.  F.3.  This  general  procedure  can  be 
extended  to  an  "EIFS  fit”  but  appropriate  equations  would 
have  to  be  developed. 
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APPENDIX  G 


FRACTOGRAPH I C  DATA  PROCESSING  METHODS 

G . 1  INTRODUCTION 

Data  processing  methods  are  needed  to  prepare  raw  frac- 
tographic  data  for  durability  analyses  because  fractographic 
data  for  each  fatigue  crack  in  the  data  set  may  not  be  homo¬ 
geneous*  For  example,  the  fractographic  data  for  a  given 
data  set  may  not  cover  the  same  crack  size  range.  The  mini¬ 
mum  and  maximum  crack  size  readings  typically  vary  from 
crack-to-crack  (e.g.,  see  Fig.  G.l).  Fractographic  data  pro¬ 
cessing  methods  are  needed  to:  (1)  obtain  a  homogeneous 

fractographic  data  set  and  (2)  extrapolate  and/or  interpolate 
fractographic  results  to  define  TTCIs  for  a  selected  refer¬ 
ence  crack  size  (x1) .  Methods  are  described  and  discussed  in 
this  section  for  processing  raw  fractographic  data  for  dura¬ 
bility  analyses. 


G.2  RAW  FRACTOGRAPHIC  DATA 

Raw  fractographic  data  situations  are  described  in  Fig. 
G.l.  Various  aspects  are  discussed  in  this  section. 

In  Fig.  G.l  the  fractographic  results  for  six  different 
fatigue  cracks  are  depicted.  The  fractographic  crack  size 
range  for  determining  the  "pooled  Q"  value  is  denoted  by 
AL-AU.  The  reference  crack  size  for  TTCIs  is  denoted  by  x^ 
The  following  observations  are  based  on  Fig.  G.l:  (1)  frac¬ 
tographic  results  for  crack  (J)  cover  the  AL-AU  range,  (2) 
cracks  (2)  through  ©  have  fractographic  data  in  the  AL-AU 
range  but  the  minimum  and/or  maximum  crack  size  reading  do 
not  cover  the  AL-AU  range,  (3)  there  is  no  fractographic 


data  in  the  AL-AU  range  for  crack*  ®  and  ©  ,  (4)  the 
TTCIs  for  cracks  and  (2)  can  be  obtained  by  interpola¬ 
tion,  and  (5)  the  fractographic  results  for  cracks  CD 
through  (§)  must  be  extrapolated  to  determine  the  TTCIs  at 
x^  In  Fig.  G.l,  dashed  lines  indicate  fractographic  data 
extrapolations  to  the  AL  and/or  AU  boundary. 

The  “pooled  Q"  value  for  a  given  data  set  and  the  EIFSD 
parameters  are  identified  for  a  given  AL-AU  crack  size  range 
(e.g.,  0 . 01"-0 . 05" ) .  In  some  cases  there  may  be  no  fracto¬ 
graphic  data  in  the  AL-AU  range.  Such  cases  are  illustrated 
by  cracks  (D  and  ©  in  Fig.  G.l.  To  utilize  the  fracto¬ 
graphic  data  for  cracks  ©  and  ©  ,  a  default  crack  size 
range  is  selected,  DL-DU,  so  that  the  fractographic  results 
within  the  DL-DU  range  can  be  extrapolated  for  applications 
in  the  desired  AL-AU  range. 

G. 3  EXTRAPOLATION  METHODS 

A  reasonable  method  is  described  in  this  section  for 
extrapolating  fractographic  data.  This  method  is  based  on 
the  crack  size-time  relationship  given  in  Eq.  C-3  and  the 
following  least  squares  fit  procedure.  Other  extrapolation 
methods  could  also  be  used.  For  example,  refer  to  extrapola¬ 
tion  methods  reviewed  by  Penny  and  Marriott  [60]. 

In  Fig.  G.2  the  fractographic  data  (i.e.,  a(t)  vs.  t 
results)  for  a  given  fatigue  crack  are  denoted  by  solid  cir¬ 
cles.  Suppose  the  TTCI  value,  t1#  is  beyond  the  range  of  the 
fractographic  data  an  extrapolation  is  required  to  determine 
tx  at  x1< 

The  generalized  crack  size-time  relationship  given  in 
Eq.  c-3  can  be  transformed  into  a  linear  least  square  fit 
form  as  shown  in  Eq.  C-5.  In  this  ase,  Q^,  is  the  deter¬ 
ministic  crack  growth  rate  parameter  for  a  given  fatigue 
crack. 
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The  other  terms  in  Eg.  C-5  have  already  been  defined.  A 
least  squares  fit  procedure  can  be  used  to  determine  the  con¬ 
stants  and  a(0)  in  Eg.  C-5.  Expressions  for  and  a(0) 
are  given  in  Eqs.  C-6  and  G-l,  respectively.  The  sums  in  Eg. 
G-l  can  be  determined  using  the  fractographic  data. 
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Once  the  constants  Qi  and  a(0)  have  been  determined,  the  TTCI, 
t^,  can  be  determined  from  tne  generalized  crack  size-time 
relationship  expressed  in  terms  of  t^.  Such  an  expression 
for  tx  is  given  in  Eg.  G-2. 


(G-2) 


Equation  G-2  can  be  expressed  in  terms  of  x^  as  shown  in 
Eg.  G-3 .  This  equation  can  be  used  to  extrapolate  the 
fractographic  results  to  a  selected  t1  beyond  the  range  of 
the  fractographic  data. 

Xx  =  a(0)ap(Qtlj  (G-3) 


The  Q1  value  in  Eq.  C-5  can  also  be  used  in  the  crack 
growth  rate  model,  given  in  Eq.  1  where  b  -  1)  ,  to  estimate 
the  crack  growth  rates  at  a(t),  t  points  beyond  the  range  of 
the  fractographic  data. 


G-4  DISCUSSION 

Fractographic  results  may  be  limited  or  non-existent  for 
the  desired  durability  analysis  application.  To  effectively 
use  available  fractographic  data,  reasonable  data  processing 
procedures  should  be  used  to  obtain  a  homogeneous  data  set. 
The  fractographic  data  should  be  processed  in  a  consistent 
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and  reasonable  manner.  In  any  case,  beware  of  extrapolating 
far  beyond  the  limits  of  the  actual  fractographic  data.  Com¬ 
puter  software  is  available  for  screening  and  processing  the 
raw  fractographic  data  [ 6 ] . 
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APPENDIX  H 


INITIAL  FLAW  SIZE 

Tha  following  are  discussed  in  this  section:  (1)  ini¬ 
tial  flaw  size  needs  for  durability  analysis,  (2)  the  meaning 
and  limitations  of  EIFS  in  a  durability  analysis,  (3)  EIFSD 
functions  with  an  "open0  and  "closed"  upper  tail,  (4)  how  in¬ 
itial  flaw  sizes  can  be  determined  from  the  EIFSD  function, 
and  (5)  expressions  for  initial  flaw  sizes  for  four  EIFSD 
functions.  The  size  of  the  initial  flaw  depends  on  the  EIFSD 
function  used,  the  desired  probability  of  crack  exceedance 
and  the  EIFS  upper  bound  limit,  xu. 

H.l  DURABILITY  ANALYSIS  NEEDS 

Functional  impairment  design  requirements  for  excessive 
cracking  and  for  fuel  leakage/ ligament  breakage  can  be  evalu¬ 
ated  using  the  entire  EIFSD  for  a  population  of  structural 
details.  In  this  case,  the  entire  EIFSD  is  grown  forward  to 
determine  the  probability  of  crack  exceedance  and/or  func¬ 
tional  impairment. 

Functional  impairment  due  to  fuel  leakage  and  ligament 
breakage  can  also  be  evaluated  using  an  initial  flaw  size  ob¬ 
tained  from  the  EIFSD.  In  this  case,  the  durability  analysis 
is  similar  to  the  typical  damage  tolerance  analysis.  Ex¬ 
cessive  cracking  is  concerned  with  relatively  small  fatigue 
cracks  for  the  entire  population  of  structural  details. 
Whereas,  fuel  leakage  and  ligament  breakage  are  generally 
concerned  with  selected  critical  control  points  within  a 
structure  where  this  type  of  functional  impairment  may  occur. 
In  this  case,  large  through-the-thickness  type  cracks  (e.g., 
0.50"  -  0.75")  may  cause  functional  impairment. 
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There  are  subtle  differences  between  using  a  single  ini¬ 
tial  flaw  size  and  the  entire  EIFSD  for  durability  analysis. 
When  a  single  initial  flaw  size,  x^,  is  used  the  initial  flaw 
is  obtained  from  the  EIFSD  and  it  corresponds  to  a  selected 
crack  exceedance  probability  (e.g.,  0.001).  Suppose  the  in¬ 
itial  flaw  size,  x^  ,  is  grown  forward  to  a  crack  size  x1  at 
time  t.  Then,  it  can  be  said  in  this  case  that  the  proba¬ 
bility  of  exceeding  the  crack  size  x1  is  p(i,T)  »  0.0001, 
which  is  the  same  probability  associated  with  the  selected 
initial  flaw  size.  On  the  other  hand,  suppose  the  entire 
EIFSD  is  to  be  grown  forward  to  time  t  the  same  way  the 
single  initial  flaw  size  x^  was  grown.  In  this  case,  the 
probability  of  exceeding  a  crack  size  x1  at  time  t  would  be 
the  same  as  that  based  on  the  single  initial  flaw  growth. 

The  entire  EIFSD  should  be  used  to  evaluate  functional 
impairment  due  to  excessive  cracking.  In  this  case,  the  pro¬ 
bability  of  exceeding  a  specified  crack  size  xx  at  time  t  is 
the  question.  Since  p(i,T*)  varies  depending  on  the  x^^ 
value,  a  single  initial  flaw  size  x^  cannot  be  used  to  make 
crack  exceedance  predictions  for  different  x^^  values. 

H.2  MEANING  AND  LIMITATIONS  OF  EIFSs 

The  meaning  and  limitations  of  an  initial  flaw  obtained 
from  the  EIFSD  are  as  follows: 

1.  The  initial  flaw  is  artificial  -  not  a  real  flaw, 
such  as  a  scratch,  burr  or  void  in  a  fastener  hole.  There¬ 
fore,  the  size  of  artificial  initial  flaws  are  not  neces¬ 
sarily  comparable  with  minimum  fractographic  crack  sizes  ob¬ 
tained  from  tear-down  inspection  results.  Furthermore,  an 
artificial  initial  flaw  cannot  be  verified  by  NDI. 
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2.  An  artificial  flaw  defines  a  oharactariatic  crack 
aiza  that  appliaa  to  diffarant  crack  ahapaa  and  geometries. 

3.  An  artificial  initial  flaw  dapanda  on  tha:  (l) 
fr act ©graphic  crack  8ize  ranga  uaad,  AL-AU,  (2)  crack  growth 
rate  modal  uaad  for  back-extrapolations  (DCGA  or  SCGA) ,  (3) 
CIFSO  function  used  and  (4)  upper  bounda  limits  (xu) .  These 
factors  define  the  "basis"  for  the  EIFS. 

4.  All  artificial  Initial  flawa  or  EIFSs  must  be  con¬ 
sidered  in  the  same  context  that  such  flaws  ware  dafined. 
For  example,  an  EIFS  must  be  grown  forward  for  crack  growth 
predictions  in  the  same  manner  the  EIFS  was  determined.  In 
other  words,  when  an  EIFS  is  grown  forward  in  a  crack  growth 
analysis  the  EIFS  is  subject  to  the  same  limitations  and  con¬ 
ditions  used  to  define  the  EIFSs  in  the  first  place.  This  is 
a  very  important  point.  For  example,  once  an  artificial  ini¬ 
tial  flaw  has  been  determined  for  a  given  crack  exceedance 
probability,  this  initial  flaw  cannot  be  indiscriminately 
grown  forward  to  any  crack  size  uaing  a  LEFM  crack  growth 
program  [e.g.,  38,  39]  without  considering  the  "basis"  for 
the  EIFSs. 


5.  EIFSs  and  EIFSDs  are  defined  using  an  "engineering 
approach."  This  approach  does  not  reflect  a  mechanistic-based 
characterization  of  the  complex  microcracking  process.  Dura¬ 
bility  "analysis  tools"  are  needed  for  design  purposes.  The 
"engineering  approach"  satisfies  the  needs  for  durability 
analysis  during  the  design  stage. 

H.3  COMPARISON  OF  EIFSD  FUNCTIONS 

Two  types  of  EIFSD  functions  will  be  discussed:  (1) 

distribution  function  with  closed  upper  tail  (e.g.,  Weibull 
compatible  and  lognormal  compatible)  and  (2)  distribution 
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function  with  open  upper  tail  (e.g,,  two-parameter  Waibull 
and  lognormal) .  Tha  two  typaa  of  EZFSD  functiona  are  ahown 
in  Fig.  H.l. 

H.3.1  EIFSD  Function  with  Closed  Upper  Tail 

A  compatible  type  EIFSD  function  is  derived  using  a  re¬ 
asonable  cumulative  distribution  function  for  TTCI  and  an  ap¬ 
propriate  deterministic  crack  growth  law.  In  this  case,  the 
resulting  EIFSD  is  said  to  be  "compatible  with  the  TTCID. " 
As  shown  in  Fig.  H.l (a),  the  compatible  type  EIFSD  function 
has  a  closed  upper  tail.  This  means  that  the  largest  initial 
flaw  size  in  the  EIFSD  is  equal  to  the  EIFS  upper  bound 
limit,  xu<  In  most  cases,  the  initial  flaw  size  of  most  in¬ 
terest  lies  in  the  extreme  upper  tail  of  the  EIFSD  function. 
For  example,  suppose  the  initial  flaw  size,  x^,  corresponding 
to  a  probability  of  crack  exceedance  p(i,T),  of  1  out  of 
10000  (0.0001)  is  desired.  In  Fig.  H.l(a),  p(i,T)  is  equal 
to  the  cross-hatched  area  under  the  EIFSD.  It  is  seen  in 
Fig.  H.l (a)  that  the  initial  flaw  size,  x^  approaches  xu  as 
a  limit  at  the  pCi,^)  approaches  0.  Thus,  xu  has  a  strong 
influence  on  tha  extreme  initial  flaw  size  desired  for  durab¬ 
ility  analysis. 

H.3.2  EIFSD  Function  with  Open  Upper  Tail 

In  some  cases,  the  user  may  want  to  assume  an  EIFSD 
function  for  representing  the  IFQ  which  has  an  open  upper 
tail  as  shown  in  Fig.  H.l(b).  For  example,  the  two-parameter 
Weibull  and  the  lognormal  distribution  functions  have  this 
characteristic . 

A  compatible-type  EIFSD  function  is  preferred  because  it 
has  a  sounder  basis  than  the  non-compatible  type.  For  exam¬ 
ple,  the  compatible-type  EIFSD  function  is  based  on  a  physi- 
cally-meaningful  TTCID  function  which  models  the  fatigue 
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wearout  process  and  tranda.  On  tha  othar  hand,  if  tha  EIFSD 
function  ia  a imply  assumed  without  regard  to  tha  fatlgua 
waarout  process,  tha  raaulting  SIFSD  may  not  ba  physically 
meaningful . 

In  Fig.  K.l  tha  "closed  upper  tail"  EIFSD  is  conceptu¬ 
ally  compared  with  tha  "open  upper  tail"  EIFSD.  Tha  "open 
upper  tail"  type  EIFSD  has  three  general  characteristics: 
(1)  since  there  is  no  upper  bound  limit  on  the  EIFSs,  infi¬ 
nite  initial  flaw  sizes  are  possible  and  (2)  the  shape  of  the 
EIFSD  is  typically  skewed  downward  whereas  the  compatible- 
type  EIFSD  is  skewed  upward  and  (3)  when  the  EIFSD  is  grown 
forward  to  predict  the  cumulation  distribution  of  TTCI, 
FT(t),  a  truncated  TTCID  is  required  to  assure  that  all  TTCls 
will  be  positive  (e.g.,  see  Fig.  H.2). 

H.3.3  Initial  Flaw  Size  Expressions 

Initial  flaw  size  expressions  for  four  different  func¬ 
tions  are  summarized  in  Table  H.l.  Once  the  EIFSD  parameters 
have  been  determined  for  a  given  EIFSD  function,  an  initial 
flaw  size  can  be  determined  for  a  given  probability  of  crack 
exceedance.  For  example,  an  initial  flaw  size,  x^,  corre¬ 
sponding  to  a  probability  of  crack  exceedance,  p(i,7*)  ,  for  1 
out  of  10000  (or  0.0001)  can  be  determined.  Further  consid¬ 
erations  are  required  to  define  the  initial  flaw  size  for  a 
given  probability  and  confidence  level. 
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Figure  H.2.  Xrgnaform»tl°n  °f  Open  Upper  Tall  Type  EIFSD 
Obtain  Cumulative  distribution  of  VTSTl - 


TABLE  H.l.  Summary  of  Expressions  for  Initial  Flaw  Size  for 
Four  Different  EIFSD  Functions. 
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